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I.  SUMMARY 


The  high-temperature  deformation  behavior  of  TiB2/NiAl 
composites  were  investigated  through  the  measurement  and  analyses 
of  the  thermal  activation  parameters  of  the  deformation  process, 
such  as  the  activation  volume,  the  apparent  activation  energy, 
and  the  temperature  dependence  of  yield  stress.  The  stress 
exponent,  n,  and  the  apparent  activation  energies  for  deformation 
of  TiB2/NiAl  composites  are  functions  of  stress,  and  volume 
fractions  of  the  reinforcements.  These  apparent  activation 
energies  are  considerably  higher  than  that  for  self-diffusion  in 
NiAl.  Also,  it  strongly  appears  that  the  deformation  process  is 
not  controlled  by  diffusion.  However  it  may  be  assisted  by 
diffusion.  The  yield  stress  strongly  depends  on  the  temperature 
for  all  0,  10,  and  20  V%  TiB2/NiAl  composite,  in  the  temperature 
range  of  703-1273  K.  A  20  V%  TiB2/NiAl  composites  is  about  two 
and  a  half  times  stronger  than  a  0  V%  TiB2/NiAl.  Activation 
volumes  of  deformation  in  the  temperature  range  of  813-1273  K  for 
TiB2/NiAl  composites  are  similar  to  that  of  monolithic  NiAl  and 
resembles  the  low-temperature  behavior  of  b.c.c.  metals. 
Deformation  of  TiB2/NiAl  composites  at  765  to  1273  K  exhibit 
low-temperature  characteristics,  and  grain  size  refinement  in 
these  composites  could  account  for  the  strengthening  of  the 
composites. 

Interfaces  of  TiB2“NiAl  and  a-Al203-NiAl  in  TiB2/NiAl 
composites  have  been  investigated  by  analytical  electron 


microscopy.  Although  no  consistent  crystallographic  orientation 
relationships  have  been  found  between  NiAl  and  TiB2  or  AI2O3, 
semicoherent  interfaces  between  a-Al203  and  NiAl  have  been 
observed  by  HREM  in  areas  where  the  low  indexed  crystallographic 
planes  of  a-Al203  aligned  with  that  of  NIAl.  No  semicoherent 
interfaces  between  NiAl  and  TiB2  have  been  observed.  Silicon 
segregation  was  consistently  detected  by  EDS  at  the  TiB2/NiAl 
interface  region.  Segregation  has  not  been  detected  in  the 
a-Al203-NiAl  interface  region.  The  segregation  layer  observed  at 
the  TiB2-NiAl  interface  is  too  thin  to  absorb  any  of  the  thermal 
residual  stress. 
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Abstract 

I  he  uddiikm  of  2(1  vol.  "h  dtaniiini  diboride  in 
particulate  form  (1-2  pm)  to  nickel  aluminide 
lTiH,-.\iAl)  results  in  a  twofold  increase  in  the 
high  temperature  strength  of  NiAl.  There  are  at 
least  two  theories  that  have  been  proposed  to 
account  for  the  high  temperature  strength  of  dis¬ 
continuous  reinforced  metal  matrix  composites. 
However,  they  cannot  be  adequately  used  as  a 
basis  to  explain  the  observed  strengthening. 

.■\n  investigation  was  undertaken  of  NiAl. 
l()vol."i:TiB--NiAl  and  2(hvl.%TiB;.-NiAl  in  the 
annealed  condition  and  after  deformation,  allowed 
to  cool  slowly.  There  is  a  low  dislocation  density  in 
the  annealed  samples  and  the  dislocation  density 
did  increase  slightly  as  a  result  of  deformation. 
However,  deformation  did  produce  some  intriguing 
dislocation  arrangements;  for  example,  it  was 
found  that  there  was  a  high  dislocation  density 
within  the  TiR,  in  the  deformed  higher  volume 
fraction  composites  and  the  dislocation  density 
within  the  NiAl  matrix  iva5  not  uniform. 

I.  Introduction 

l.I.  Background 

The  intcre.st  in  intermeiallic  conipound.s  for 
high  temperature  applications  extends  back  in 
time  for  a  number  of  years  j  1  ]  and  there  has  been 
a  continued  interest  since  then  [2-26].  There  are 
numerous  obvious  advantages  for  the  considera- 
ti  'n  of  aluminides  and  silicides.  e.g.  they  have 
good  oxidation  resistance  at  high  temperatures, 
ffowever.  they  are  not  widely  used  in  structural 
applications,  but  they  are  ( e  g.  NiAl  t  used  as  coat¬ 
ings  on  superalloys  for  oxidation  resistance.  The 
reason  that  NiAl  is  not  used  in  high  temperature 
structural  applications  is  that  its  strength  is  at 
least  a  factor  of  10  weaker  than  conventional 
superalloys.  For  example,  at  1 300  K  and  a  strain 
rate  of  10" '  s‘ '  the  stress  required  for  UDIMET 
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I  1.^  127i  is  about  20(1  MPa  whereas  for  Ni.4l  it  is 
less  than  20  MPa  [26j. 

.Another  deterrent  to  the  use  of  NiAl  in  struc¬ 
tural  applications  is  its  lack  of  ductility,  especially 
in  the  large  grain  size  condition  [13].  The  fracture 
toughness  as  measured  by  Charpy  impact  testing 
is  about  0.08  ft  lb  (0.1  J)  I28j.  It  was  assumed  that 
the  fracture  toughness  of  SiC  measured  by 
Charpy  impact  testing  was  low.  but  in  that  case 
the  fracture  toughness  is  1  ft  lb  [29j. 

The  lack  of  high  temperature  strength  and  the 
low  temperature  ductility  can  prtibabb  be 
improved  by  the  addition  of  TiB,.  i.e.  to  produce 
a  discontinuous  composite  of  TiB--NiAl.  The 
major  advantage  of  TiB;  is  that  it  is  inert  with 
respect  to  the  NiAl  matrix  [30],  and  in  high  tem¬ 
perature  applications  the  matrix  should  not  react 
with  the  reinforcement.  If  a  reaction  takes  place, 
it  is  possible  that  reinforcement  will  be  consumed 
during  the  reaction. 

1.2.  .Mechanism  or  mechanisms  of  strengthening 

The  high  temperature  creep  of  composites  as 
an  indication  of  their  high  temperature  strength 
has  been  analyzed  in  terms  of  a  classical  contin¬ 
uum  composite-.strengthening  mechanism  |31j. 
Simply,  this  means  that  there  is  a  reduction  of  the 
stress  in  the  matiix  owing  to  the  load  transfer  to 
the  reinforcement.  However,  if  these  mechanisms 
are  applied  to  a  discontinuous  spherical  reinforce¬ 
ment.  the  load  transfer  is  very  small,  i.e.  5%. 
which  would  mean  that  the  predicted  decrease  in 
the  creep  rate  of  the  spherical  reinforced  com- 
(K)site  would  be  very  small  compared  with  that  (>t 
the  matrix  material. 

If  one  attempted  to  perform  a  more  rigorous 
analysis  by  the  use  of  the  Eshelby  method,  there 
is  a  restriction  which  is  presently  placed  on  the 
Eshelby  method  which  requires  that  the  matrix 
has  to  be  treated  as  a  newtonian  fluid,  i.e.  the 
stress  exponent  n  would  have  to  be  unity.  The 
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observed  stress  exponents  are  much  greater  than 
unity,  i.e.  in  the  range  5-20(32]. 

Another  approach  would  be  to  consider  the 
dislocation  mechanisms  that  have  been  proposed 
to  predict  creep  behavior.  However,  very  few- 
detailed  mechanisms  have  been  proposed  which 
take  into  account  the  presence  of  large  second- 
phase  particles  [32, 33]. 

From  a  consideration  of  a  dislocation  climb 
over  second-phase  particles,  it  is  possible  to 
obtain  the  following  equation  [32]  for  the  creep 
rate; 

i  -  '  ^  ’  /''  ’/>  ^Ep„,  I ’ 

D~  2.83,-r  rkT  [e) 


where  D  is  the  diffusion  coefficient  for  the  matrix, 
i'  is  Poisson's  ratio,  Z  is  the  interparticle  spacing, 
E  is  the  Young  modulus  of  the  matrix,  is  the 
density  of  mobile  dislocations,  is  the  effective 
climb  stress,  b  is  the  Burgers  vector,  r  is  the 
radius  of  the  particle,  k  is  Boltzmann’s  constant 
and  T  is  the  temperature.  If  it  is  assumed  that  the 
second  phase  is  approximately  spherical,  then 
there  is  a  simple  relationship  [33]  between  /. 
and  r: 


where  I '.  is  the  volume  fraction  of  the  reinforce¬ 
ment,  As  one  would  guess  intuitively,  for  a  con¬ 
stant  Vf.  e/D  is  greater  for  a  larger  size 
reinforcement  and  decreases  with  an  increase  in 
the  volume  fraction  of  reinforcement.  If  the  actual 
f  vs.  is  determined  for  a  composite  of  20  vol.% 
and  particles  of  radius  0.5  pm  in  the  NiAl  matrix, 
then  eqn.  { I )  predicts  a  stress  of  less  than 
6.9  X  10  '  MPa  for  an  f  of  10  '  s’ '.  This  stress 
is  at  least  three  orders  of  magnitude  smaller  than 
the  strength  of  the  matrix.  There  are  other  pos¬ 
sible  dislocation  models  which  have  been  pro¬ 
posed  by  Blum  and  Reppich  133).  but  the  result  is 
that  if  the  Orowan  stress  is  low  for  the  composite, 
the  effect  of  the  reinforcement  on  the  creep  rate  is 
small. 

There  have  been  a  few  transmission  electron 
microscopy  (TEM)  investigations  of  dislocations 
overcoming  large  particles  at  high  temperatures, 
but  details  of  the  rate-controlling  mechanisms  are 
not  given  [34], 


It  can  be  concluded  that  there  is  a  large  differ¬ 
ence  in  the  predicted  high  temperature  strengths, 
Le.  as  manifested  by  the  creep  rates  of  discontinu¬ 
ous  composites  as  predicted  by  the  continuum 
mechanics  and  dislocation  models.  If  we  consider 
the  preliminary  data  obtained,  there  is  approxi¬ 
mately  a  factor-of-2  increase  of  the  strength  of 
the  composite  (20  vol.%j  ix  the  matrix. 

The  purpose  of  this  investigation  was  to  deter¬ 
mine  the  microstructures  of  TiB.-NiAl  samples 
which  were  deformed  at  high  temperatures,  and 
from  these  determinations  to  obtain  some  insight 
into  the  possible  strengthening  mechanism. 

2.  Materials 

The  materials  used  in  this  investigation  were 
purchased  from  the  Martin  Marietta  Corporation 
and  consisted  of  0.  10  and  20  vol.%  TiB-  in  NiAl. 
The  materials  were  ground  into  compression 
.samples  with  a  height-to-diameter  {ft/d)  ratio  of 
1.5-2,  which  were  then  compression  tested  at 
1140  and  1300  K  to  various  strain  levels.  The 
samples  were  also  annealed  at  the  test  tempera¬ 
ture  for  periods  of  time  from  2  to  1 2  h. 

3.  Testing  methods 

The  samples  were  tested  in  compression  at 
various  strain  rates  ranging  from  lO"*  to  10'- 
s  ■ '  at  11 40  and  1 300  K.  A  decision  w-as  made  to 
deform  the  samples  well  past  initial  yielding,  since 
if  there  was  a  substantial  increase  in  dislocation 
density,  the  assumption  was  made  that  there 
would  be  a  greater  chance  of  retaining  this  higher 
density,  because  at  the  completion  of  testing  the 


Fig,  1.  The  yield  stress  at  H)00°C  1 1273  K)  of  TiB,-NiAI 
composites  as  a  function  of  strain  rate  and  volume  fraction  of 
TiB.-NiAl  sample. 
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samples  were  not  cooled  rapidly;  i.e.  they  were 
cooled  over  a  15  min  period  and  not  in  a  time 
frame  of  less  than  2  min  [35]. 

After  compression  testing  aryi/or  after  anneal¬ 
ing.  the  samples  were  cut  into  3  mm  discs  with  a 
thickness  of  0.5  mm  by  electrical  discharge 
machining.  After  lightly  grinding  to  remove  the 
debris  on  the  surface,  the  discs  were  dimple 
ground  until  the  thickness  reached  0.07-0.1  mm 
in  the  center.  1  he  final  electropolishing  was  per¬ 
formed  using  5%  perchloric  acid-ethanol  solu¬ 
tion  at  253  K  or  below  until  a  small  hole 


Fij:.  2.  SH  M  micrograph  of  a  20vol.%TiB;-NiAI  sample. 


appeared  at  the  center  of  the  sample.  The  applied 
voltage  was  50-70  V. 

An  investigation  was  also  carried  out  both 
optically  and  by  scanning  electron  microscopy 
(SEM)  to  determine  the  panicle  size  and  distribu¬ 
tion.  This  was  then  followed  by  an  X-ray  powder 
diffraction  investigation  to  determine  the  content 
ofTiB,. 

The  TEM  was  carried  out  at  the  high  voltage 
electron  microscopy  (HVEM)  facility  at  .Argonne 
National  Laboratory. 

4.  Experimental  results  and  discussion 

The  yield  stress  a^.  which  is  only  slightly  less 
than  the  ultimate  stress,  was  determined  as  a 
function  of  volume  fraction  of  TiB.  and  stram 
rate.  As  expected,  o,  increases  with  increasing 
volume  fraction  (Fig.  1 )  and  increasing  strain  rate. 
However,  the  increase  in  with  increasing 
volume  fraction  of  TiB;  is  greater  at  faster  strain 
rates. 

Figure  2  is  an  SEM  micrograph  of  a  2()vol.‘/o- 
TiB;-NiAl  sample.  There  are  particles  other  than 
TiB,  and  these  have  been  tentatively  identified  as 
a-Al;0-,  from  an  examination  of  TEM  .selected 
area  diffraction  patterns.  The  dark  particles  in 
Fig.  2  are  TiB;  and  the  light  particles  .  a-Al.O, 


Fig  3.  X-ray  powder  diffraction  pattern  obtained  from  a  30vol.“oTiB;-NiAI  powder  sample. 
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can  constitute  as  high  as  10  vol.%.  Figure  3  is  an 
X-ray  powder  diffraction  pattern  obtained  from  a 
20%  TiB;  powder  sample.  The  purpose  of  the 
X-ray  investigation  was  to  obtain  some  insight  as 
to  what  were  in  the  SEM  micrographs,  i.e.  an 
attempt  to  determine  the  nature  of  the  other  par¬ 
ticles. 

Figures  4  and  5  are  TEM  micrographs  of  0 
vol.%  samples  deformed  30%  at  1 143  K  (870  °C/. 
It  is  evident  that  the  second-phase  panicles  are 
acting  as  barriers  to  dislocation  motion.  Figures 
6  and  7  are  again  of  ()vol.%TiB,-NiAI  matrix 
samples  and  show  ihe  formation  of  a  network 
which  is  prevalent  in  these  alloys.  It  is  also  evident 
that  the  dislocation  structure  is  not  planar,  i.e. 
there  is  no  well-defined  slip  plane,  and  the  dislo¬ 
cations  have  convoluted  shapes  as  shown  in  Fig. 


t-ij;  4  TtM  micrograph  of  a  (Ivol  "..TiB;-NiAI  sample, 
.til",,  deformed,  mosi  of  ihe  dislocaiion  having  a  Burgers 
seeior  (>f  a  i  IflOi 


Fig.  5.  TbM  micrograph  of  a  (lvol,”..TiB.-NiAI  sample. 
.t(l%  deformed,  most  of  the  dislocation  having  a  Burgers 
vector  of  a  |  UMI|. 


8.  This  result  agrees  with  the  observations  on 
crept  .single-crystal  NiAl  [36). 

Another  phenomenon  which  is  only  observed 
in  0  vol.%  NiAl  matrix  samples,  especialb  those 
which  were  deformed  at  higher  stram  rates,  are 
bands  or  regions  of  very  high  dislocation  density 


Fig.  6.  TEM  micrograph  of  a  Ovol.%TiB.-NiAI  sample. 

deformed.  Networks  formed  bv  reactions  of 
!l(tOj  +  iOI(l|  =  |l  I0)and  |OU)l  +  !OO)l  =  (0lij. 


Fig  7.  TEM  micrograph  of  a  livol.''.,TiB.-NiAI  sample, 
.'ll".,  delormed.  Networks  tormed  b\  reactions  of 
,1110' s- it) 111!  =  i  1  lOiandlOlOi-t- 1110)  1  =  101  i 


V 


(.K^IN  )t«H  MUM\ 

Vivi.  5^.  Schemaiic  illuMTauon  lor  convoluied  shape  disU'ca- 
lums. 
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which  we  define  as  shear  bands.  They  could  be 
the  same  type  of  bands  which  were  observed  [36] 
in  crept  single-crystal  NiAl  by  using  surface  rep¬ 
licas.  Figure  9  is  an  example  of  these  shear  bands 
in  a  sample  which  was  deformed  at  1140  K.  The 
shear  bands  appear  to  go  straight  across  the  grain 
boundary  without  any  apparent  shifting  of  the 
grain  boundary  or  of  the  shear  band. 

If  we  now  consider  the  samples  containing 
TiB;  particulate,  several  things  become  evident. 
Figure  10  is  a  low  magnihcation  TEM  micro¬ 
graph  of  20vol.7oTiB2-NiAl  which  was  deformed 
30“/o.  There  is  a  very  high  dislocation  density  in 
the  area  with  higher  than  average  volume  fraction 
of  TiB,.  whereas  in  regions  of  tower  than  average 
TiB;  concentration  there  is  a  very  low  dislocation 
density.  Also,  it  is  evident  in  these  pictures  that 
within  the  TiB,  there  are  a  large  number  of 
dislocations. 

This  difference  in  dislocation  density  which 
depends  upon  the  local  volume  fraction  of  1  iB,  is 
evident  in  Figs.  11  and  12  where  there  are  two 
different  localized  volume  fractions  of  TiB,. 


Figure  1 1  contains  a  very  high  dislocation  density 
in  the  matrix  whereas  in  Fig.  12  there  is  a  very 
low  dislocation  density  in  the  matrix. 

It  is  generally  observed  that  there  are  local¬ 
ized  regions  of  high  dislocations  in  the  deformed 
(10-20)vol.7oTiB2-NiAl  composites.  However, 
if  the  average  dislocation  density  is  considered, 
then  for  the  composite  the  dislocation  density  is 
less  than  the  average  dislocation  density  in  the 
deformed  Ovol.%TiB;-NiAl  matrix  material. 

Returning  to  a  consideration  of  the  plastic 
deformation  which  occurs  in  the  TiB^,  there  are 
several  unique  observations  as  shown  in  Figs.  13 
and  1 4.  First,  the  TiB;  particles  are  not  uniformly 
deformed  in  all  cases.  Secondly,  in  the  TiB^  par¬ 
ticles  which  contain  a  variable  dislocation  density 
there  is  a  variable  dislocation  density  in  the  sur¬ 
rounding  matrix.  There  is  an  inverse  relationship; 
if  a  portion  of  a  particle  has  a  high  dislocation 
density,  then  the  dislocation  density  in  the  matrix 
immediately  surrounding  this  portion  of  the  par¬ 
ticle  has  a  low  dislocation  density;  if  a  portion  of 
the  particle  has  a  low  dislocation  density,  then  the 


Fig.  9.  TEM  micrograph  of  a  Ovol.'/oTiBj-NiAl  sample.  30%  deformed.  Shear  bond  structure. 


t-  ig.  1 0  TEM  micrograph  of  a  20vol.%TiB;-NiAI  sample.  .''0%  deformed. 


Fig  1  I  TEM  micrograph  of  a  20vol.‘’..TiB;-NiAI  sample.  Fig  12.  TEM  micrograph  of  a  20vol.” oTiB.-NiAl  sample. 
I".,  deformed.  Highly  conceniraied  TiB;  area.  30”» deformed.  Low  concemration  of  TiB.. 


matrix  immediately  surrounding  this  portion  of 
the  particle  ha.s  a  high  dislocation  density. 

As  stated  in  the  introduction,  the  purpose  of 
this  investigation  was  to  determine  the  changes 
in  microstructure  of  TiB;-NiAI  composites 


deformed  at  high  temperatures  for  the  purpose  of 
obtaining  insight  into  possible  strengthening 
mechanisms  in  this  composite.  The  results 
obtained  do  not  appear  to  support  any  of  the 
presently  proposed  mechanisms.  The  theories  of 
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Fii!.  ]3.  TEM  micrograph  of  a  2()vol.%TiB;-NiAI  sample. 

deformed.  Non-uniformly  deformed  TiB,  panicles. 
Deformed  matrix  always  in  contact  with  dislocation-free  end 
of  TiB;  panicle. 


Fig.  14.  TEM  micrograph  of  the  same  area  as  in  Fig.  13, 
except  deformed  end  of  TiB-  panicle  with  a  practically  dislo¬ 
cation-free  area  surrounding  matrix. 


Taya  and  Lilholl  |3 1 )  which  are  ba.sed  on  the  load 
transfer  concept  require  that  there  should  be  zero 
or  very  little  plastic  deformation  of  the  reinforce¬ 
ment.  This  requirement  does  not  .seem  to  be  satis¬ 
fied  owing  to  high  dislocations  within  the  TiB- 
particles.  The  theories  of  dislocation  climb  based 
on  climb  over  second-phase  particles  are  inade¬ 
quate,  since  they  predict  .strengthening  which  is 
three  orders  of  magnitude  too  .small  [32.  33). 
Finally,  the  theories  of  composite  strengthening 
proposed  by  Arsenault  and  coworkers  [37-45) 
which  are  based  on  an  increased  dislocation 
density  and  reduced  subgrain  size  cannot  be 
applied.  The  average  dislocation  density  in  the 
1 0-20  vol.7o  composites  is  less  than  the  average 
dislocation  density  in  the  0  vol.%  matrix  material. 


Therefore  an  entirely  new  approach  to  causes  of 
strengthening  will  have  to  be  considered. 

5.  Conclusions 

From  a  consideration  of  the  data,  there  are 
several  conclusions  at  which  we  can  arrive. 

First,  in  regard  to  the  NiAl  matroc 

i  1 !  there  are  other  particles  present  within  the 
matrix,  possibly  constituting  as  much  as  1 0  \  ol.‘’n; 

i2i  the  Burgers  vector  is  a  J 100!  type  with  no 
defined  slip  plane; 

l3i  deformation  at  low  strain  rates  results  in 
no  large  change  in  dislocation  densities,  while  at 
higher  strain  rates  shear  bonds  are  observed  in 
the  Ovol.%TiB;-NlAl. 

Secondly,  in  regard  to  the  TiB. -NiAl  com¬ 
posites; 

i  1 1  there  is  heterogeneous  distribution  of  dis¬ 
locations  in  the  TiB, -NiAl  deformed  samples; 

i2)  the  TiB,  is  deformed; 

',3 1  at  pre.seni.  a  strengthening  mechanism  is 
not  obvious. 
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Abstract 

In  this  investigation,  it  was  shown  that  titanium  diboride  would  plastically 
deform  at  a  temperature  of  1273  K  and  at  a  stress  level  much  lower  than  predicted 
by  previous  experimental  data.  The  crystal  structure  of  titanium  diboride  is  a 
topologically  hexagonal  structure  which  allows  for  a  number  of  possible  slip 
systems.  Experimentally,  it  was  observed  that  three  slip  directions  were  present: 
[0001],  ]<2TT0>  and  j<2TT3>.  The  probable  slip  planes  are  (0001),  |lofo|  and 

;ioTi|. 


()l.  Introduction 

Titanium  diboride  (TiBj)  is  an  interesting  material.  It  has  a  high-temperature 
strength  which  is  greater  than  other  ceramics  such  as  silicon  carbide  (SiC)  and 
aluminium  oxide  (AljOiKRamberg  and  Williams  1987).  Also,  TiBj  is  a  reinforcement 
candidate  in  high-temperature  composite  for  it  is  inert  with  respect  to  potential  matrix 
materials.  One  of  the  composite  systems  which  has  received  some  attention  is  TiB,  in 
nickel  aluminide  (NiAl)  (Wang  and  Arsenault  1990).  It  was  found  that  the  addition  of 
30  vol.",,  particulate  TiBj  resulted  in  a  threefold  increase  in  the  yield  stress  at  1273  K. 
The  reason  why  this  increase  occurs  is  not  evident. 

A  microstructural  investigation  was  undertaken  of  these  composites  to  determine 
the  effect  of  the  addition  ofTiBj  on  the  dislocation  arrangement  and  density  within  the 
NiAl  matrix.  It  was  assumed  that  the  TiBj  did  not  deform  at  1 273  K.  for  the  yield  stress 
at  1273K  which  is  obtained  by  a  lower  limit  extrapolation  would  be  2000 MPa 
(Ramberg  and  Williams  1987).  The  stress  concentration  in  the  TiBj  that  could  occur 
owing  to  load  transfer  would  not  be  large  enough  to  exceed  2000  MPa.  However,  it 
should  be  pointed  out  that  the  TiBj  used  by  Ramburg  and  Williams  (1987)  contained 
titanium  carbide  (TiC)  precipitates. 

In  the  investigation  of  dislocation  configuration  and  density  (Wang  and  Arsenault 
1 990)  in  the  NiAl  matrix  ofTiBj  NiAl  composites  deformed  at  1273  K,  it  was  observed 
that  there  was  a  high  density  of  dislocations  within  TiBj  particles.  Therefore,  at  least 
two  questions  arose. 

( 1 )  Did  the  dislocations  within  the  TiBj  arise  from  deformation  of  the  compiosites? 

(2)  What  types  of  dislocation  are  within  the  TiBj  particles? 

There  appears  to  be  a  large  number  of  possibilities. 

Nakano,  Imura  and  Takenchi  (1973)  and  Nakano,  Matsubar  and  Imura  (1974)  in 
an  investigation  of  the  high-temperature  hardness  ofTiBj  concluded  that  the  primary 
slip  system  was  { I010}<T2T0>.  However,  Ramberg  and  Williams  proposed  two  other 
systems:  {10T0}[{)001]  and  (000l)<I0T0>. 
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L.  Wang  and  R.  J.  Arsenault 
Fig.  1 


Ball  model  of  TtB2.  Boron  atoms  form  a  net  connection  only.  Titanium  atoms  sit  In  the  open 
centre  of  the  net  (Westbrook  1967). 


Table  I.  Yield  stress  of  10vol.%TiB2-NiAI  deformed  at  1273  K. 


Strain  rate 
(s') 

Yield  stress  (MPa) 

Annealed  at  1673  K 

Not  annealed 

11  X  10  “ 

809 

706 

2'2x  10  “ 

84-6 

81  6 

2'2x  10  ■’ 

1361 

1113 

The  reason  that  so  many  possibilities  exist  is  related  to  the  crystal  structure  of  TiBj 
which  is  schematically  shown  in  fig.  I.  The  crystal  structure  could  be  defined  as  a 
topologically  hexagonal  structure  (Parthe  1967),  in  which  the  boron  atoms  form  a  net 
connection  only.  For  the  case  of  packing  of  touching  spheres  (Laves  1956),  the  ideal 
c,a=  107,  and  c/a=  1063  for  TiBj.  From  an  examination  of  fig.  1,  it  is  evident  that 
there  is  no  well  defined  close-packed  plane,  but  the  closest-packed  plane  is  (0001)  and 
the  shortest  repeat  distance  in  this  plane  is  ^[1120]. 

The  purpose  of  the  present  investigation  was  to  determine  whether  an  increase  in 
the  dislocation  density  did  indeed  occur  upon  deformation  of  the  TiB2  NiAl 
composite,  and  to  characterize  the  various  types  of  dislocation  Burgers  vector  that  were 
assumed  to  be  present  in  the  TiBj  particles  within  the  deformed  TiBj-NiAl  composite. 

§2.  Materials  and  experimental  procedures 
The  TiBj  NiAl  composite  used  in  this  investigation  was  purchased  from  Martin 
Marietta  Corporation.  The  composite  contained  10  or  20vol.%TiB2  and  the  TiB2 
particles  were  of  various  shapes  and  size  but  on  average  the  largest  dimension  was  3  pm. 

Prior  to  testing,  the  compression  samples,  with  a  length-to-diameter  ratio  of  two. 
were  annealed  at  1273  and  1673  K  for  periods  of  time  of  2  12  h.  Variation  in  annealing 
times  had  no  effect  on  the  strength  of  the  composite.  Surprisingly,  annealing  at  1673  K 
resulted  in  a  slight  increase  in  strength.  The  samples  were  tested  at  1273  K  at  various 
strain  rates  and  a  summary  of  the  results  is  given  in  table  I  (details  of  the  testing  have 
been  given  elsewhere  (Wang  and  Arsenault  1990)).  The  transmission  electron 
microscopy  (TEM)  foils  were  obtained  from  the  deformed  composite  samples  and 
examined  at  both  l(X)kV  and  1  MeV.  Further  details  of  the  foil  preparation  have  been 
given  elsewhere  (Wang  and  Arsenault  1990).  Since  it  is  assumed  that  there  will  be  a  large 
number  of  possible  Burgers  vectors,  it  is  necessary  to  use  a  number  of  operating 
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diffraction  vectors.  However,  if  the  relative  intensity  of  the  operating  diffraction  vector 
is  low,  then  the  dislocations  would  be  in  very  low  contrast.  Therefore  all  micrographs  in 
this  paper  were  taken  under  operating  conditions  where  the  relative  intensity  is  large 
enough  to  show  the  difference  between  the  residue  contrast  of  invisible  dislocations  and 
the  image  of  visible  dislocations. 

§  3.  Results  and  discussion 

The  initial  portion  of  the  investigation  was  to  determine  whether  the  dislocation 
density  within  the  TiBj  increased  as  a  result  of  a  deformation  of  the  composite.  In 
annealed  samples,  there  is  a  low  dislocation  density  in  both  the  matrix  and  the  TiB^,  as 
shown  in  fig.  2  (a),  except  in  the  case  where  some  of  the  particles  are  in  contact  with  each 
other.  A  comparison  of  TEM  micrographs  of  annealed  samples  and  deformed  samples 
clearly  shows  that  the  density  of  dislocations  within  the  TiB2  is  significantly  higher  in 
the  deformed  composite  sample,  as  shown  in  figs.  2  and  3. 

It  was  found  that  the  planar  sides  of  a  TiBj  particle  were  almost  always  parallel  to 
the  (0001 )  plane.  The  dislocations  labelled  A  and  B  do  not  lie  on  the  plane  of  the  foil  but 
lie  on  a  plane  whose  trace  is  parallel  to  the  bottom  edge  of  the  particle.  In  other  words, 
dislocations  A  and  B  lie  on  the  ((XIOl)  plane  which  is  perpendicular  to  the  plane  of  the 
foil.  Under  conditions  of  (X)0l  and  202T  operating  diffractions,  dislocations  labelled  A 
in  fig.  3  (a)  are  almost  completely  invisible.  Therefore,  the  dislocations  labelled  A  have 
Burgers  vectors  of  |[1 2 10].  The  dislocation  labelled  Bin  figs.  3(a)  and  (c)  which  is  only 
invisible  under  the  operation  diffraction  0001  has  a  Burgers  vector  of  ^[2110]  or 
K1I30]. 

If  the  electron  beam  is  exactly  parallel  to  the  basal  plane,  then  the  projection  of 
dislocations  on  the  basal  plane  should  be  straight  lines.  However,  there  are  exceptions  if 
climb  occurs  as  indicated  by  the  arrow  in  fig.  3(a).  Also,  it  has  been  found  that  stacking 
faults  exist  on  the  basal  plane  in  TiB^  (Mersol,  Lynch  and  Vahldiek  1968).  Again,  if  the 
electron  beam  is  parallel  to  the  basal  plane  as  in  fig.  4(a),  the  dislocation  on  the  basal 
plane  should  project  as  a  straight  line  labelled  A.  If  the  sample  is  tilted  as  shown  in 
fig.  4(6),  then  stacking-fault  ribbons  are  evident.  The  presence  of  stacking  faults  is  not 
universal;  only  a  few  TiBj  particles  contained  stacking  faults. 

It  is  suspected  that  slip  on  the  prism  plane  also  occurs  and  that  prismatic  punching 
into  the  TiBj  particles  from  the  interface  takes  place,  that  is  the  dislocations  labelled  C 
in  figs.  5(a)  and  (6).  If  different  operation  diffractions  arc  used  (TOlO  and  3l  10),  then 
prismatically  punched  dislocations  become  invisible  (figs.  5  (c)  and  [d ))  in  areas  that  the 
Burgers  vector  of  prismatically  punched  dislocations  is  [0001]. 

There  are  dislocations  with  other  Burgers  vectors  within  this  particle.  The 
dislocations  labelled  E  in  figs.  5(a)  and  (c)  are  out  of  contrast  under  the  diffraction 
vector  0(X)1  (fig.  4(d)).  The  dislocations  which  have  Burgers  vectors  of  3<2I10> 
obviously  do  not  lie  on  the  basal  plane.  It  is  therefore  possible  that  slip  directions  of  the 
]<2l  10>-type  can  occur  on  other  slip  planes. 

Besides  dislocations  with  Burgers  vectors  of  ]<2ll0>  and  [(X)0l],  there  are  still 
other  dislocations  with  different  Burgers  vectors,  for  example  the  dislocations  labelled 
D  in  fig.  5.  It  was  difficult  experimentally  to  find  two  independent  extinction  conditions 
for  these  dislocations  because  of  the  tilting  limit  of  the  sample  stage.  Every  operating 
diffraction  condition  was  employed  which  had  an  intensity  high  enough  to  show  the 
difference  between  the  residual  contrast  of  invisible  dislocations  and  the  contrast  of 
visible  dislocations  in  the  (I2l0)*,  (OTIO)*,  (Olll)*  and  (I2l3)*  reciprocal  planes. 
However,  the  Burgers  vector  of  dislocations  labelled  D  could  not  be  determined. 


(6) 

TEM  micrographs  of  (a)  annealed  20vol.%TiB2~NiAI  composite  (low  dislocation  density  both 
in  matrix  and  TiBj)  and  (h)  a  TiBj  particle  in  annealed  20  vol.^TiBj-NiAl  composite. 
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TEM  micrographs  of  a  TiBj  particle  in  deformed  20vol.%TiB2-NiAl  composite  with  different 
operating  diffractions:  (a)  dislocation  with  Burgers  vector  of  ^<2l  10>  type  on  the  basal 
plane  are  visible;  (b)  all  i<2l  10>-type  dislocations  are  invisible;  (c)  only  dislocations  with 
Burgers  vector  of  j[l2l0]  are  invisible. 


It  is  possible  to  obtain  the  Burgers  vector  of  dislocations  labelled  D  by  considering 
the  dislocation  reaction  which  occurs  at  the  arrow  in  figs.  5(b)-(d).  Dislocations 
labelled  A  which  have  a  Burgers  vector  of  |[2l  10]  reacted  with  the  dislocation  labelled 
D  which  results  in  the  short  segment  of  dislocation  indicated  by  the  arrow  in  these 
pictures.  This  short  segment  of  dislocations  has  a  Burgers  vector  of  [0001],  for  it 
followed  the  same  extinction  conditions  as  the  dislocation  labelled  C  which  has  a 
Burgers  vector  of  [0001].  Therefore  the  reaction  can  be  written  as 

l>A  +  l>D~^C- 

The  subscripts  denote  the  dislocations  labelled  in  fig.  5.  Since  is  10]  and  be  is 
[0001],  then  the  dislocations  labelled  D  must  have  a  Burgers  vector  of  the  type 
K3l  13).  The  reason  why  it  is  so  difficult  to  find  two  independent  extinction  conditions 
to  determine  the  Burgers  vector  of  dislocation  labelled  D  directly  is  that  only  the  (TOT I ) 
diffraction  condition  satisfies  the  extinction  condition.  This  one  extinction  condition  is 
not  enough  to  determine  the  Burgers  vector  of  these  dislocations.  It  was  impossible  to 
obtain  extinction  conditions  owing  to  the  limitations  of  the  sample  tilting  stage. 

The  possible  combinations  of  slip  directions  and  slip  planes  are  listed  in  table  2.  The 
three  directions  were  determined  in  this  investigation,  and  the  planes  listed  are  the  three 
of  highest  atomic  packing. 


(b) 

TEM  micrographs  of  a  TiBj  particle  in  deformed  20vol.%TiB2-NiAl  composite  with  stacking 
faults  on  the  basal  plane;  (a)  basal  plane  is  parallel  to  the  electron  beam;  (b)  basal  plane  is 
not  parallel  to  the  electron  beam. 


id) 

TEM  micrographs  of  a  TiBj  particle  in  deformed  20  vol.%TiB2-NiAI  composite;  (o)many-beam 
condition,  all  the  dislocations  are  visible;  (ft)  all  the  }<3ll0>-type  dislocations  are 
invisible:  (c),  id]  prismatically  punched  dislocation  out  of  contrast. 
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Table  2.  Slip  systems  in  titanium  diboride. 


Slip 

direct 


Presence  of  slip  system  in 
the  following  slip  planes 


(0001)  {lOTO}  {1011} 


[0001]  Yes 

K2110>  Yes  Yes  Yes 

K2113>  Yes 


§4.  Conclusions 

From  the  data  obtained  in  this  investigation,  it  is  possible  to  arrive  at  the  following 
conclusions. 

(1)  The  TiBj  particles  within  a  deformed  TiB^  NiAl  composite  can  be  plastically 
deformed  as  shown  by  the  increased  dislocation  density  within  the  TiBj 
particles. 

(2)  Three  slip  directions  were  observed  within  the  TiBj  particles;  they  are  ^<2l  10>, 
[0001]  and  ]<2113>. 

(3)  Dislocations  with  Burgers  vectors  of  the  type  ^<2 1 10>  on  the  basal  plane  could 
split  into  two  partials  with  a  stacking  fault  between  the  partials. 

(4)  Most  of  the  dislocations  in  TiB2  had  Burgers  vectors  of  [0001]  or  13>. 
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Abstract 

Tensile  tests  of  SiC-6061  At  composites  con¬ 
taining  various  volume  fractions  of  whiskers  or 
particles  (20,  5  and  0  vol.%)  showed  that  for 
samples  containing  a  high  volume  fraction  (20 
vol.  %)  the  fracture  process  was  very  localized,  i.e.  a 
very  narrow  neck.  As  the  volume  fraction  of 
whiskers  or  particles  decreased,  the  deformed 
region  spread  out.  One  might  expect  that  the 
microstructure  should  correspond  to  the  macro¬ 
scale  changes.  In  the  highly  deformed  region  the 
dislocation  density  is  expected  to  be  higher,  in  the 
less  deformed  regions  the  dislocation  density 
should  be  lower,  and  if  the  deformation  is  very 
localized,  then  the  high  dislocation  density  should 
also  be  limited  to  a  very  narrow  region. 

Overall,  there  is  good  agreement  between  the 
microstructure  (dislocation  density)  change  and 
the  macroscale  deformation  of  SiC-Al  composite 
tensile  samples.  The  mechanism  proposed  to 
account  for  this  change  in  deformation  behavior  as 
a  function  of  volume  fraction  of  SiC  in  aluminum 
is  related  to  the  expansion  of  the  plastic  zone  (due 
to  differences  in  thermal  coefficients  of  expansion 
between  SiC  and  aluminum)  when  the  external 
stress  is  applied.  Also,  the  localized  deformation  is 
related  to  localized  clusters  of  SiC  particles.  There 
is  a  cooperative  effect  which  leads  to  a  region  of 
very  localized  plastic  deformation. 


“necking”  is  observed  (in  high  volume  fraction 
composites),  but  the  fracture  surface  has  both 
brittle  and  ductile  characteristics.  There  is  no 
indication  of  fracture  of  SiC  when  the  size  of  SiC 
is  less  than  10  /rm.  The  fracture  may  initiate  at 
large  intermetallic  inclusions,  clusters  of  SiC 
particulates  or  whiskers  and  voids  associated  with 
the  SiC  particulates  or  whiskers  [1], 

The  results  [2]  of  tensile  tests  of  SiC-6061  A1 
composites  with  various  volume  fractions  of 
whiskers  or  particulates  (20,  5  and  0  vol.%)  indi¬ 
cated  that  for  high  volume  fraction  samples  (20 
vol.%)  the  fracture  process  was  very  localized.  As 
the  volume  fraction  of  whiskers  or  particulates 
decreased,  the  deformed  region  “spreads  out”. 
Figure  1  demonstrates  how  “necking”  changes 

Omm 
3  mm 
5mm 


0.0% 

0.5% 

3% 


(a) 


(b) 


1.  Introduction 

There  are  many  unusual  deformation  charac¬ 
teristics  associated  with  SiC-Al  discontinuous 
metal  matrix  composites.  One  of  them  is  that  the 
ductility  and  fracture  toughness  of  SiC-Al  com¬ 
posites  are  very  low  [1],  Generally,  examinations 
of  the  fracture  surface  do  not  reveal  the  presence 
of  SiC  particulates  or  whiskers  in  the  quantities 
corresponding  to  their  volume  fraction  [1].  The 
fracture  process  is  very  localized,  very  little 
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Fig.  1 .  Reduction  in  area  at  various  distances  from  fracture 
surface  for  SiC-Al  composite  tensile  samples:  (a)  un¬ 
deformed;  (b)  20  vol.%  SiC;  (c)  5  vol.%  SiC;  (d)  0  vol.%  SiC. 
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with  a  decrease  in  volune  fraction.  In  the  case  of 
the  20  vol.%  composite  the  reduction  in  area  at 
the  fracture  surface  was  only  3%,  and  5  mm  away 
from  the  fracture  surface  the  reduction  in  area 
was  zero.  However,  in  the  case  of  the  5  vol.% 
composite,  at  the  fracture  surface  the  reduction  in 
area  was  14%,  and  5  mm  away  from  the  fracture 
surface  the  reduction  in  area  was  4%. 

One  might  expect  in  a  highly  deformed  region 
that  the  dislocation  density  is  very  high  and  in  the 
less  deformed  region  the  dislocation  density 
should  be  lower.  In  other  words,  the  high  disloca¬ 
tion  density  should  be  limited  to  a  very  narrow 
region.  It  is  also  of  interest  to  know  that  previous 
work  [1]  has  indicated  that  along  the  fracture  path 
the  local  density  of  SiC  particles  is  higher  than  the 
average  volume  fraction  of  the  SiC.  In  other 
words,  the  fracture  path  is  attracted  to  the  par¬ 
ticle  clusters.  However,  it  is  still  unclear  whether 
such  a  “particle  friendliness”  is  only  restricted  to 
the  fracture  process  or  if  it  also  governs  the  onset 
of  plastic  deformation. 

If  there  are  high  dislocation  densities  in  regions 
containing  a  larger-than-average  volume  fraction 
of  SiC,  this  indicates  that  is  where  deformation  is 
taking  place.  This  is  intuitively  difficult  to  under¬ 
stand.  It  has  been  clearly  shown  that  if  the  volume 
fraction  of  SiC  is  increased  in  an  aluminum 
matrix,  the  strength  increases.  Therefore  the 
question  arises;  why  should  deformation  begin 
and  be  restricted  to  regions  of  higher-than- 
average  volume  fraction  of  SiC?  The  logical 
assumption  would  suggest  that  deformation 
should  begin  in  regions  of  less-than-average 
volume  fraction  of  SiC. 

There  have  been  numerous  investigations  and 
several  theories  of  the  ductile  fracture  of  two- 
phase  alloys  where  the  second  phase  is  a  discrete, 
non-deforming  precipitate  or  particle  [3-19].  All 
the  theories  are  based  on  the  assumption  that  a 
void  nucleates  at  the  interface  of  a  non-deforming 
particle  and  then  grows,  leading  to  an  ultimate 
fracture.  An  examination  of  transverse  sections 
by  optical  microscopy  and  scanning  electron 
microscopy  (SEM)  methods  indicated  that  there 
are  very  few  voids  near  the  fracture  surface.  It  is 
possible  that  the  voids  nucleated  are  very  small 
and  can  only  be  observed  by  transmission 
electron  microscopy  (TEM)  techniques. 

To  answer  the  questions  raised  above,  an 
investigation  was  proposed  to  determine  the 
dislocation  density  along  the  tensile  axis  of  the 
gauge  section  by  using  transmission  electron 


microscopy  to  relate  macro-  and  microstructural 
changes  to  the  plastic  deformation  and  to  deter¬ 
mine  if  a  correlation  exists  between  slip  line 
generation  and  density  and  the  localized  volume 
fraction  of  SiC  in  the  SiC-Al  composite.  Also,  the 
purpose  of  this  investigation  was  to  determine  if 
there  was  a  higher  density  of  voids  near  the  frac¬ 
ture  surface  in  SiC-Al  composites. 

To  investigate  the  underlying  mechanism  for 
such  a  phenomenon,  a  set  of  two-dimensional 
finite  element  method  (FEM)  meshes  were  gener¬ 
ated  to  simulate  the  effect  of  inhomogeneous 
distribution  of  the  SiC  reinforcement  whiskers, 

1. e.  the  effect  of  clustering.  Furthermore,  it  is  also 
of  importance  to  know  the  plastic  deformation 
process  around  each  SiC  whisker.  Therefore  a 
three-dimensional  FEM  mesh  was  generated  to 
monitor  the  behavior  of  plastic  zones  around 
whiskers. 

2.  Experimental  procedures 

2.1.  Materials 

SiC-whisker-  or  particulate-reinforced  6061 
AI  alloy  composites  which  were  obtained  com¬ 
mercially  were  used  for  this  investigation. 

The  composites  with  20,  5  and  0  vol.%  SiC 
whiskers  (SiC,,-Al)  were  purchased  from  ARCO 
SILAG  (presently.  Advanced  Composite  Materi¬ 
als)  in  the  form  of  extruded  rods  15.5  mm  in 
diameter.  The  composites  with  20  and  5  vol.% 
SiC  particulates  (SiCp-Al)  were  purchased  from 
DWA  Composite  Specialties  in  the  form  of  plates 
15.9  mm  thick.  All  materials  were  fabricated  by 
the  powder  metallurgy  method.  The  whiskers  [20] 
were  small  in  diameter,  ranging  from  0.1  to  1  /im, 
and  the  average  aspect  ratio  was  2.  The  particu¬ 
lates  were  platelet  shaped,  5-7  /<m  long  and  had 
an  aspect  ratio  of  2-3. 

Sample  blanks  were  cut  from  the  rods  and 
plates,  then  machined  into  tensile  samples.  The 
doubly  reduced  samples  were  required  to  ensure 
that  the  fracture  occurred  in  the  gauge  section 
and  not  in  the  transition  region  [20].  The  dimen¬ 
sions  of  the  samples  are  given  in  Fig.  2.  Those 
that  were  to  be  heat  treated  were  given  a  modified 
T6  heat  treatment,  which  means  a  solutionizing 
step  of  3  h  (modified)  as  opposed  to  the  standard 
1  h.  The  modified  heat  treatment  was  used  in 
case  a  divorced  eutectic  formed  during  the  pro¬ 
cessing.  The  others  were  annealed  for  1 2  h  at  a 
solutionizing  temperature  of  810  K  and  furnace 
cooled.  The  method  of  producing  the  “melted” 
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Fig.  2.  Dimensions  of  a  round  tensile  sample  (in  millimeters). 


250  /<m  SiCp-1 100  A1  alloy  composites  has  been 
described  elsewhere  [21]. 

2.2.  Tensile  testing 

All  samples  were  tested  in  an  Instron  testing 
machine  at  room  temperature.  Most  of  the  sam¬ 
ples  were  tested  to  failure;  however,  some  were 
tested  to  a  strain  where  a  measurable  "neck” 
occurred  in  the  sample. 

2.3.  TEM  thin  foil  preparation  and  examination 

The  method  of  TEM  foil  preparation  has  been 

given  elsewhere  [22].  In  a  recent  publication  by 
Arsenault  et  al.  [23]  it  has  been  clearly  shown  that 
in  order  to  obtain  representative  dislocation  den¬ 
sities  it  is  necessary  to  examine  the  foils  at  a  high 
operating  voltage.  All  the  foils  were  observed  in 
the  AEI-EM7  high  voltage  electron  microscope 
(HVEM)  operated  at  1  MeV  at  the  Argonne 
National  Laboratory. 

The  TEM  foil  thickness  was  determined  by 
counting  the  extinction  fringes  or  from  stereo- 
pairs.  All  the  micrographs  were  taken  under  two 
beam  conditions  with  the  diffraction  vectors 
being  g  =  [lll]  or  [200],  but  in  most  cases  the 
diffraction  vector  was  g  =  { 1 1 1 ) .  If  a  higher  order 
of  diffraction  vectors  such  as  g  =  ]220)  is  used, 
the  contrast  is  reduced,  which  results  in  difficul¬ 
ties  in  determining  the  dislocation  density.  Since 
the  Burgers'  vector  of  aluminum  is  of  H 1 10)  type 
and  there  are  six  specific  Burgers'  vectors,  any 
given  1111}  reflection  will  result  in  three  of  the  six 
being  visible.  Therefore  in  calculating  the  disloca¬ 
tion  density  the  measured  density  was  multiplied 
by  a  factor  of  2. 

For  each  TEM  foil  examined  in  the  HVEM  a 
minimum  of  10  micrographs  were  taken  and 
there  could  be  as  many  as  100.  The  procedure 
was  to  examine  the  micrographs  and  make  a 
visual  selection,  removing  the  micrographs  which 
had  extremely  high  or  low  dislocation  densities. 
This  selection  procedure  greatly  reduced  the  time 


necessary  to  determine  the  dislocation  density  of 
a  given  foil. 

2.4.  Dislocation  density  calculation 
The  dislocation  densities  from  the  selected 
micrographs  were  determined  by  the  method 
proposed  by  Keh  [24],  In  this  method  a  grid  of 
two  sets  of  lines,  normal  to  each  other  but  with 
different  spacing,  is  placed  over  the  TEM  micro¬ 
graph  and  then  the  intersection  of  dislocations  is 
given  by 


NJL.+N.JL, 

P  = - - -  (1) 

where  L,  and  L,  are  the  lengths  of  the  grid  lines 
normal  to  each  other,  N,  and  N,  are  the  inter¬ 
sections  along  the  lines  and  /  is  the  thickness  of 
the  foil  in  the  area  where  the  micrograph  was 
taken. 

The  average  density  was  determined  for  a 
given  foil  by  simply  taking  the  numerical  average 
of  the  examined  foils.  This  .same  procedure  was 
followed  for  all  foils  of  a  given  condition,  e.g.  the 
same  distance  from  the  fracture  surface.  The 
error  bar  associated  with  data  points  of  disloca¬ 
tion  vx  distance  from  the  fracture  surface  simply 
covers  the  range  of  the  averages. 

The  upper  limit  [25]  of  dislocation  density 
calculated  by  using  eqn.  (1)  is  about  10'-'-10‘^ 
m"-.  In  .some  micrographs  the  dislocations  were 
so  den,se  that  it  was  impossible  to  see  the  indi¬ 
vidual  dislocation.  The  dislocation  density  was 
estimated  by  using  the  width  of  the  micrograph 
divided  by  the  width  of  the  dislocation  image. 
Clearly,  the  value  of  dislocation  density  was 
underestimated. 

In  the  case  of  very  high  dislocation  densities  it 
may  be  argued  that  the  weak  beam  technique 
should  be  employed.  However,  the  weak  beam 
technique  fails  because  it  is  extremely  sensitive  to 
the  diffraction  condition,  i.e.  the  s  vector.  For  a 
sample  with  a  very  high  dislocation  density  the  s 
vector  does  change  from  point  to  point  within  the 
sample.  This  is  especially  true  in  the  thicker  foils 
which  are  used  in  the  HVEM. 

2.5.  Slip  line  analysis 

While  composite  samples  with  a  small  reinforce¬ 
ment  size  are  more  suitable  for  TEM  analyses, 
composites  with  larger  reinforcement  are  pre¬ 
ferred  for  optical  observation  of  slip  lines.  As  a 
result,  melted  250  pm  particulate-reinforced 
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composite  material  was  chosen  for  the  analysis. 
From  this  melted  particulate-reinforced  com¬ 
posite,  rectangular  samples  8  mm  wide  by  50  mm 
long  and  1.5  mm  thick  were  electrical  discharged 
machined  and  polished.  The  samples  were  glued 
to  an  E.  Fullam  tensile  substage  model  18202 
JEOL  JSM-35.  The  substage  was  mounted  on  a 
specially  constructed  holder  to  fit  on  to  a  Zeiss 
ICM-405.  The  entire  effective  gauge  was  photo¬ 
graphed.  Then  the  sample  was  deformed  by  a 
smdl  amount  and  the  area  which  contained  slip 
lines  was  re-photographed,  and  again  the  sample 
was  deformed  and  re-photographed  where  slip 
lines  occurred.  This  process  was  repeated  until 
the  sample  fractured. 

2.6.  FEM procedures 

2.6. 1.  Two-dimensional  FEM  procedures 
In  the  present  study  the  composites  are  repre¬ 
sented  by  a  periodic  array  of  transversely  aligned 
whiskers  as  shown  in  Fig.  3,  i.e.  an  infinite  two- 
dimensional  array  of  SiC  whiskers  embedded  in 
the  aluminum  matrix  in  the  plane  strain  condi¬ 
tion.  In  order  to  consider  the  effect  of  clustering, 
the  array  is  transformed  in  such  a  manner  that  an 
array  of  periodic  geometric  cluster  centers  was 
selected  and  the  neighboring  whiskers  (in  this 
case,  four  whiskers  (Fig.  4))  were  attracted  to¬ 
ward  their  cluster  centers  so  that  a  periodic  array 
of  clusters  were  formed  as  shown  in  Fig.  5.  The 
degree  of  clustering  was  defmed  by  the  values  of 
SJ  V  and  SJH  in  Fig.  5. 

In  applying  the  above  geometric  transforma¬ 
tion  scheme,  a  commercially  available  ADINA 
FEM  code  was  employed.  To  represent  the 
reality,  the  material  is  assumed  to  be  stress  free  at 
the  annealing  temperature,  i.e.  773  K,  and  is 
subsequently  cooled  to  room  temperature. 
Therefore  the  effect  of  thermal  residual  stresses  is 


implicitly  considered.  After  thermal  treatment  the 
composites,  clustered  and  unclustered,  were 
uniaxially  loaded  in  the  whisker  directions. 
Owing  to  the  periodic  arrangement,  the  condition 
of  symmetry  can  be  so  utilized  that  a  “unit  cell” 
was  selected  as  a  basic  building  block  such  that 
only  a  quadrant  of  whisker  (cluster)  is  contained, 
as  shown  in  Figs.  6(a)  and  6(b),  and  multicon¬ 
straint  boundary  conditions  were  employed  to 
represent  the  symmetry  conditions  at  the 
boundary. 

2. 6. 2.  Three-dimensional  FEM  procedure 

In  this  case  the  reinforcement  particles  are 
hexagonally  distributed  through  three-dimen¬ 
sional  space.  Because  of  the  symmetry,  the 
response  of  the  composite  can  be  fully  repre¬ 
sented  simply  by  considering  a  building  block 
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Fig.  4.  This  figure  shows  how  SiC  whiskers  move  to  their 
geometric  cluster  center  to  form  a  periodic  array  of  clusters. 
Whiskers  with  dotted  boundaries  represent  the  final  loca¬ 
tions  of  the  whiskers  after  the  movement  indicated  by  the 
arrows  in  the  figure,  and  5.  and  5^  show  the  degree  by  ver¬ 
tical  and  horizontal  clustering. 


Fig.  3.  Periodic  hexagonal  arrangement  of  the  SiC  reinforce¬ 
ment  of  the  whisker.  The  shaded  areas  are  SiC  whiskers. 


I . I 


Fig.  5.  Configuration  of  clustered  composites. 
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Fig.  6.  (a)  Unit  cell  for  the  uniform  distribution  of  whiskers 
employed  by  the  FEM  analysis,  where  “h"  and  “v"  stand  for 
the  horizontal  and  vertical  half-spacing,  (b)  Unit  ceil  for  the 
clustered  whisker  configuration  employed  by  the  FEM 
analysis,  where  5,  and  5^  are  defined  in  the  same  way  as  in 
Fig.  4. 


which  only  contains  one-eighth  of  the  reinforce¬ 
ment  particle,  as  shown  in  Fig.  7,  and,  in  consid¬ 
eration  of  symmetry,  multipoint  constraint 
equations  were  employed  in  the  boundary  condi¬ 
tions.  As  in  the  case  of  two-dimensional  analysis, 
the  composites  were  assumed  stress  free  at 
773  K,  which  is  the  annealing  temperature  for 


Fig.  7.  Unit  cell  employed  for  the  three-dimensional 
analysis. 


SiC-Al  composites,  and  the  sample  was  subse¬ 
quently  cooled  to  room  temperature.  Owing  to  a 
large  difference  in  thermal  expansion  coefficient 
between  SiC  and  aluminum,  severe  thermal  resid¬ 
ual  stress  would  be  generated  around  particles. 
As  a  result,  it  is  apparent  that  a  plastic  zone  will 
form  around  each  particle.  Upon  further  uniaxial 
loading,  such  a  plastic  zone  will  certainly  respond 
to  the  external  load.  Depending  on  how  such  a 
response  would  be,  it  will  be  helpful  to  under¬ 
stand  the  deformation  mechanism  of  the  com¬ 
posites.  Therefore  it  is  our  intention  to  monitor 
the  geometric  change  of  the  plastic  zone  as  a 
result  of  external  load.  This  can  be  done  simply 
by  imaginarily  slicing  the  sample  longitudinally 
along  the  loading  direction  and  monitoring  the 
two-dimensional  change  of  the  plastic  zone  in 
each  plane. 

3.  Experimental  results 

The  discussion  of  the  experimental  results  will 
be  divided  into  four  parts:  (a)  dislocation  density 
vs.  distance,  (b)  slip  line  analysis,  (c)  FEM  results 
and  (d)  voids. 

3.  J.  Dislocation  density  vs.  distance 

3. 1. 1.  20  vol.  %  SiC-Al  composite 

The  plot  of  dislocation  densities  {p)  vs.  dis¬ 
tance  from  the  fracture  surface  (x)  is  shown  in 
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Fig.  8.  Dislocation  density  i'5.  distance  from  fracture  surface 
for  20  vol.%  SiC-Al  composite. 


Fig.  9.  Dislocation  density  vs.  distance  from  fracture  surface 
for  5  vol.%  SiC-AI  composite. 


Fig.  8.  The  value  of  p  is  higher  near  the  fracture 
surface  region  and  decreases  as  x  increases. 
When  the  distance  is  3  mm  or  more,  p  is  constant 
and  the  value  of  p  is  approximately  equal  to  the 
value  of  p  for  an  undeformed  sample.  This  value 
is  about  one  order  of  magnitude  lower  than  that 
in  the  0.3-0. 5  mm  region,  i.e.  near  the  fracture 
surface  region.  In  general,  the  values  of  p  are 
higher  in  the  whisker-reinforced  heat-treated 
samples. 

The  difference  between  the  di.sk>cation  den¬ 
sities  of  each  micrograph  in  the  same  foil 
increases  as  x  increases,  which  means  the  disloca¬ 
tion  densities  are  more  uniformly  distributed  near 
the  fracture  surface  region. 

In  some  foils  the  difference  in  the  dislocation 
densities  is  quite  large  (it  may  exceed  two  orders 
of  magnitude).  To  avoid  confusion,  the  reported 
data  are  the  average  values  of  micrographs  taken 
from  each  location  rather  than  those  of  a  single 
micrograph. 

Several  attempts  were  made  to  deform  20  vol.% 
SiC-Al  composites  until  a  “neck"  occurred,  then 
stopping  the  test  prior  to  fracturing  the  sample, 
and  then  cutting  TEM  foils  from  the  sample  in 
the  "neck"  region.  All  attempts  were  failures; 
either  the  sample  fractured  before  the  test  was 
stopped  or  upon  examination  of  the  TEM  foils. 
There  was  no  difference  in  density  between  the 
perceived  “neck"  region  and  the  region  away 
from  the  “neck". 

3. 1.2.  5  vol.  %  SiC-At  composite 

The  plot  of  dislocation  densities  (p)  vs.  dis¬ 
tance  from  the  fracture  surface  lx)  is  shown  in 
Fig.  9.  Similar  to  the  20  vol.%  composite  (Fig.  8). 


at  0.3-0. 5  mm,  i.e.  near  the  fracture  surface 
region,  p  is  higher  and  approximately  equals  that 
of  the  20  vol.%  samples  in  the  same  region;  then 
p  decreases  as  x  increases.  However,  at  3  mm 
and  beyond,  p  remains  constant  and  the  values  of 
p  are  3-5  times  higher  than  that  of  undeformed 
samples,  a  result  which  is  different  from  the 
20  vol.%  samples.  Also,  in  agreement  with  the 
20  vol.%  samples,  p  is  higher  in  the  whisker- 
reinforced  heat-treated  samples. 

As  in  the  case  of  the  20  vol.%  SiC-Al  com¬ 
posite,  several  samples  were  tested  until  a  "neck" 
had  formed.  Then  TEM  foils  w'ere  cut  from  the 
“neck"  region  and  the  ends  of  the  samples.  Figure 
10(a)  is  a  TEM  micrograph  of  the  end  region  of 
the  sample  and  Fig.  10(b)  is  a  TEM  micrograph 
showing  the  high  dislocation  density  in  the  “neck" 
region.  The  variation  of  the  dislocation  density 
from  the  "neck"  region  to  the  ends  of  the  sample 
is  shown  in  Fig.  11. 

3. 1.3.  0  vol.  %  SiC-AI  composite 

The  plot  of  dislocation  densities  (p)  vs.  dis¬ 
tance  from  the  fracture  surface  (.r)  is  shown  in 
Fig.  12.  Near  the  fracture  surface  region  (0.3-0. 5 
mm),  p  is  higher,  but  the  value  of  p  is  about  one- 
third  of  that  for  the  20  vol.%  samples  in  the  same 
region  and.  unlike  the  20  and  5  vol.%  samples,  as 
X  increases,  p  decreases  at  a  much  slower  rate. 
The  values  of  p  in  this  plot  are  about  two  orders 
of  magnitude  higher  than  for  an  undeformed  sam¬ 
ple.  Similar  to  the  20  tmd  5  vol.%  samples,  p  is 
higher  in  the  heat-treated  samples. 

In  summary,  from  the  aforementioned  experi¬ 
mental  observations,  a  fact  can  be  clearly  stated 
that,  irrespective  of  the  reinforcement  morphol- 


61 


Fig.  10.  la  TEM  micrograph  of  foil  taken  from  the  end 
section  of  5  vol."  .  whisker  SiC-AI  tensile  sample  b  TEM 
micrograph  of  foil  taken  from  the  neck  section  of  5  vol.‘’o 
whisker  SiC'-.-\l  tensile  sample 
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Fig.  11.  Plot  of  dislocation  density  in  the  neck  region  and 
the  end  sections  of  a  tensile  sample  which  has  been  deformed 
to  produce  a  neck. 
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Fig.  1 2.  Dislocation  density  is.  distance  from  fracture 
surface  for  0  vol.%  SiC-AI  composite. 


ocy.  a  high  dislocation  density  is  restricted  to  the 
vicinity  of  the  fracture  surface,  which  represents 
the  highest  plastic  deformaton  site,  i.e.  plastic 
deformation  is  very  hicalized. 

3.2.  Slip  line  aniilysis 

The  initial  slip  lines  began  in  a  region  which 
contained  a  higher-than-average  density  of  SiC 
particles  i.e.  a  cluster .  The  very  first  slip  lines  in 
this  region  occur  at  the  corners  of  SiC  particles, 
as  shown  in  Fig.  13.  Slip  lines  were  observed  at 
the  corners  of  other  SiC  particles  through  the 
gauge  length  of  the  sample.  By  successive  incre¬ 
mental  loadings,  the  gauge  section  of  the  sample 
was  examined  to  determine  the  local  density  of 
slip  lines.  The  region  with  the  initial  higher-than- 
average  density  of  SiC  particles  contains  almost 
all  the  subsequent  increase  in  slip  line  density. 


Figures  14-16  are  micrographs  taken  at  succes¬ 
sive  loadings.  In  other  areas  of  the  sample  there 
are  no  or  very  few  slip  lines. 

3.3.  FEM 

3.3. 1.  Two-dimensional  FEM 

The  global  stress-strain  curves  of  the  com¬ 
posites  produced  by  the  calculation  are  shown  in 
Fig.  17.  It  is  clear  that  the  more  the  whiskers  are 
clustered  together,  i.e.  the  smaller  5,  and  5,,  in 
Fig.  6(b),  the  lower  is  the  yield  stress.  This  is 
consistent  with  the  experimental  data  described 
in  the  previous  subsection,  i.e.  plastic  deforma¬ 
tion  is  initiated  within  the  region  where  there  is  a 
higher  localized  particle  concentration  (cluster) 
and  this  clustering  results  in  a  reduction  of  com¬ 
posite  flow  stress.  This  leads  to  an  intuitively 
obvious  conclusion:  the  more  uniform  the 


Fig.  1 3.  Sequence  of  slip  line  formation  as  a  result  of  incremental  loading  along  the  direction  shown  by  the  arrow  s.  As  the  load 
increases  (from  left  to  right),  slip  lines  first  form  around  the  corners  of  the  SiC  particles;  then,  as  the  load  funher  increases,  slip 
lines  formed  around  each  particle  w  ill  propagate  and  interconnect  and  finally  form  a  set  of  densely  packed  slip  bands. 


arrangement  of  discontinuous  reinforcement,  the 
stronger  is  the  composite. 

It  should  be  pointed  out  that  while  localization 
of  plastic  deformation  is  evident  throughout  the 
experimental  observations,  as  we  described 
previously,  FEM  modeling  on  a  parallelpiped 
array  is  best  suited  for  a  whisker  composite.  The 
results  only  give  a  qualitative  approximation  for 
the  particulate  composite. 

3.3.2.  Three-dimensional  FEM 

It  was  found  that  the  plastic  zone  extended 
about  one  particle  diameter  from  the  edge  of  the 
particle,  which  is  in  agreement  with  the  experi¬ 
mental  results  [26],  and  the  nature  of  the  arrange¬ 
ment  of  the  residual  stress  along  the  longitudinal 
direction  is  such  that  at  the  whisker  end  the 
matrix  is  under  compression  while  the  rest  of  the 
matrix  is  under  tension.  Also,  it  was  found  that 
the  plastic  zone  expanded  transversely  upon 
loading  the  sample  in  tension,  which  is  in  agree¬ 
ment  with  the  slip  line  analysis,  whereas  longi¬ 
tudinally  it  contracts  near  the  whisker  ends  owing 
to  the  compressive  nature  of  the  plastic  zone  in 
that  part  of  the  matrix.  However,  this  process 


generally  means  the  generation  of  dislocations  of 
opposite  sign;  annihilation  of  dislocations  does 
not  make  a  significant  contribution  to  the  process. 

3.4.  Voids 

The  TEM  micrographs  were  examined  to 
determine  the  density  of  voids  at  the  interface 
between  the  matrix  and  the  SiC  whiskers  or  par¬ 
ticles.  Voids  were  found;  however,  the  density  of 
voids  near  the  fracture  surface  was  greater  than 
the  density  of  voids  in  the  non-deformed  regions 
of  the  samples. 

4.  Discussion 

The  discussion  of  the  experimental  results  will 
be  related  almost  exclusively  to  explaining  how 
and  why  the  deformation  is  so  localized  in  the 
SiC-Al  composites.  There  will  be  a  brief  discus¬ 
sion  of  the  lack  of  voids. 

The  evidence  obtained  from  the  dislocation 
density  and  slip  line  investigations  is  overwhelm¬ 
ing:  deformation  in  the  higher  volume  fraction 
composites  is  very  localized.  From  an  analysis  of 
the  configuration  of  SiC  particles  in  the  samples 


Fig.  1 4.  1  he  >.amplo  is  slighiK  delormod.  It  can  he  seen  that  slip  lines  begin  to  emerge  from  the  particle-matrix  interface  and  the 
slip  line  densits  is  higher  where  there  is  higher  localized  particle  densits. 


which  were  used  for  the  slip  line  analysis,  it  is 
evident  that  the  slip  lines  began  in  regions  which 
contained  a  clustering  of  SiC  particles.  The  slip 
propagated  from  one  cluster  region  to  another, 
and  deformation  in  the  remainder  of  the  samples 
is  almost  zero.  The  data  obtained  from  the  TEM 
investigation  also  indicated  that  the  deformation 
is  very  localized  in  the  higher  volume  fraction 
discontinuous  composites.  In  the  case  of  the  20 
vol.‘’'o  SiC-Al  alloy  composite,  only  in  a  small 
region  less  than  3  mm  from  the  fracture  surface) 
is  there  an  increase  in  dislocation  density  over 
that  of  the  annealed  sample  in  a  sample  which  has 
been  tested  in  tension  to  failure.  The  phenom¬ 
enon  of  localized  deformation  was  duplicated  by 
an  FE.Vl  analysis. 

.As  the  volume  fraction  of  SiC  is  decreased,  the 
degree  of  localization  decreases,  i.e.  the  5  vol.'’<. 
SiC-.Al  compissite  has  a  dislocation  density  p' 


v.v.  distance  from  the  fracture  surface  (.v)  rela¬ 
tionship  which  is  approaching  that  of  the  0  vol.% 
composite.  This  is  a  change  which  one  would 
intuitively  guess  at.  and  the  reason  that  this  p  vs.  x 
relationship  is  different  from  that  of  the  20  vol.% 
SiC-Al  composite  is  due  to  the  statistical  nature 
of  the  clustering.  The  higher  the  volume  fraction 
of  the  particles,  the  greater  is  the  probability  of 
obtaining  a  cluster  containing  a  given  number  of 
panicles  and  the  greater  is  the  probability  of  find¬ 
ing  clusters  within  close  proximity. 

Although  the  FEM  analysis  can  account  for 
the  localized  deformation,  the  following  question 
still  remains  to  be  addressed:  why  should  pla.stic 
deformation  initiate  in  the  clustered  region  where 
there  is  a  local  high  volume  fraction  of  reinforce¬ 
ment  when  it  is  obvious  that  the  higher  the 
volume  fraction  of  reinforcement,  the  stronger  is 
the  material  and  the  less  prone  it  should  be  to 


Fig.  15.  Upon  further  straining,  the  dcnsiis  ul  ihc  'lip  line'  incrca'C'.  However,  they  are  limited  to  the  localized  region  where 
there  is  a  higher  particle  volume  traction.  The  arrows  indicate  ihe  direction  of  loading.i 


65 


Fig.  16.  As  the  external  load  increases,  plastic  deformation 
becomes  more  localized  and  finally  a  crack  is  formed  within 
such  a  highly  deformed  region.  iThe  arrows  indicate  the 
direction  of  loading.) 

plastic  deformation?  To  explain  this  phenomenon 
fully,  it  is  necessary  to  consider  the  fact  that 
deformation  compatibility  has  to  be  satisfied 
when  there  is  an  inhomogeneous  distribution  of 
reinforcement  particles,  i.e.  clustering.  Since  the 
clusters  are  stiffen  upon  satisfying  deformation 
compatibility,  the  stress  would  be  distributed  in 
such  a  way  that  the  clusters  would  bear  more  load 
than  the  rest  of  the  matrix.  Secondly,  since  the 
local  volume  fraction  is  higher  within  clusters,  for 
the  same  amount  of  deformation,  the  cluster  must 
activate  more  slip  systems  in  the  matrix  to  accom¬ 
modate  the  same  amount  of  deformation.  Figure 
18(a)  shows  effective  plastic  strain  as  a  function 
of  global  strain  at  various  locations  within  the 
matrix  as  shown  in  Fig.  1  S(  b;.  From  the  plot  it  can 
easily  be  seen  that  there  exists  a  strong  deforma- 


Fig.  1 7.  Stress-strain  curves  produced  as  a  result  of  cluster¬ 
ing:  curve  1,  SJ2V=  1,  SJ2H=  1;  curve  2,  S,/2V'=0.75,  SJ 
277=0.429;  curve  3,  5,/2F=0.5,  5;/27/  =  0.l43;  curve  4, 
SJ2  1^=  0.25. 5, /277  =  0.071. 


tion  gradient  and  the  plastic  strain  as  a  function  of 
external  applied  load  between  the  fiber  ends  in 
the  cluster  is  much  higher  than  elsewhere  in  the 
matrix.  The  divergence  of  effective  plastic  strain 
at  various  locations  in  the  matrix  as  the  global 
applied  strain  approaches  zero  indicates  that 
there  is  a  large  gradient  on  plastic  deformation  in 
the  matrix  owing  to  initial  thermal  residual 
stresses,  and  the  initial  flat  region  on  some  of  the 
curves  is  an  indication  that  the  initial  deformation 
in  these  regions  is  purely  elastic,  which  may 
include  possible  unloading  as  in  the  case  of  that  at 
the  whisker  tip.  Figure  18(c)  shows  a  comparison 
of  effective  plastic  strain  at  the  whisker  ends  in  a 
regular  and  a  clustered  whisker  distribution. 
From  the  figure  it  can  be  clearly  seen  that 
although  whiskers  within  a  regular  distribution  of 
whiskers  have  a  higher  plastic  deformation  rate  at 
the  whisker  end  as  compared  to  the  rest  of  the 
matrix  [27],  the  clustered  distribution  has  both  a 
higher  effective  plastic  strain  rate  and  a  higher 
effective  plastic  strain  value.  This  signifies  that 
the  matrix  material  within  the  cluster  is  prone  to 
plastic  deformation,  which  is  what  is  observed 
experimentally,  as  discussed  in  the  previous 
section. 

In  order  to  characterize  the  response  of  the 
effective  plastic  strain  to  the  external  load  at 
different  locations  in  the  matrix,  the  response  of 
an  imaginary  homogeneous  material  was  incor- 
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Fig.  18.  (a)  Localized  effective  plastic  strain  at  various  locations  vs.  global  strain.  “Homogeneous  material"  applies  to  a 
homogeneous  material  subject  to  a  tensile  load,  (b)  Locations  where  effective  plastic  strain  is  taken,  ci  Effective  plastic  strain  vs. 
global  strain  at  whisker  ends  for  clustered  and  regular  arrays  (at  location  I)  for  clustered  array  as  shown  in  (b).  It  is  clearly  seen 
that  the  clustered  array  has  a  higher  effective  plastic  strain  rate  for  the  same  global  strain  increment,  which  simply  means  the 
clustered  arrangement  of  the  reinforcement  particles  is  more  susceptible  to  external  loading  in  terms  of  plastic  deformation, 
'd)  The  deformation  rate  is  higher  within  the  reinforcement  cluster  except  at  locations  near  the  longitudinal  matrix-whisker 
interface. 


porated  into  the  plot.  Because  of  the  nature  of 
homogeneity,  the  internal  effective  plastic  strain  is 
always  equal  to  the  external  applied  global  strain. 
Furthermore,  if  localized  deformation  exists  in 
the  cluster,  the  ratio  of  the  increase  of  the  internal 
effective  plastic  strain  to  the  increase  of  the  exter¬ 


nal  applied  global  strain,  which  may  be  defined  as 
the  internal  effective  plastic  ratio,  should  be 
larger  than  that  of  a  homogeneous  material, 
which  has  a  value  of  unity.  Figure  18(d)  shows 
internal  effective  plastic  strain  rate  vs.  applied 
global  strain  at  various  locations  in  the  matrix.  It 
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is  clearly  sef  i  that  except  in  the  immediate  region 
near  the  vertical  interface  of  the  reinforcement 
and  matrix,  the  effective  plastic  strain  ratio  within 
the  cluster  is  greater  than  that  of  a  homogeneous 
material,  while  the  ratio  approaches  unity  at 
locations  away  from  the  cluster  in  the  matrix.  The 
effective  plastic  strain  ratio  of  a  homogeneous 
material  is  always  unity. 

The  lack  of  voids  at  the  matrix-whisker  inter¬ 
face  is  probably  related  to  several  phenomena. 
First,  it  has  been  shown  [28]  that  the  bond 
between  SiC  and  aluminum  is  very  good.  This 
bond  is  several  times  stronger  than  the  bond 
between  FejC  and  iron.  Secondly,  there  is  a 
region  of  cold  work,  i.e.  a  high  dislocation  density 
region,  around  the  whisker.  As  a  result  of  this 
hardened  region  around  the  whisker,  a  triaxial 
stress  develops  in  the  matrix  some  distance  away 
from  the  whisker.  It  is  into  this  region  that  the 
crack  propagates.  The  end  result  is  that  very  few 
voids  are  found  at  the  interface  between  SiC  and 
aluminum. 


5.  Conclusions 

From  the  experimental  results  and  the  FEM 
calculations  it  is  possible  to  arrive  at  several  con¬ 
clusions. 

( 1 )  In  discontinuous  metal  matrix  composites 
containing  a  high  volume  fraction  of  particle.s  and 
whiskers,  the  deformation  is  highly  localized. 

(2)  There  are  manifestations  of  three  distinct 
localized  deformations:  (a)  there  are  very  narrow 
“necks”  in  fracture  samples;  (b)  the  close-in  dis¬ 
location  density  is  confined  to  a  small  region  very 
close  to  the  fracture  surface;  (c)  the  slip  line 
generation  is  confined  to  a  localized  region. 

(3)  The  FEM  study  has  shown  that  the  effec¬ 
tiveness  of  strengthening  from  the  reinforcement 
in  the  composites  is  reduced  if  there  is  clustering 
of  the  reinforcement.  Internal  deformation  maps 
generated  by  FEM  have  shown  that  deformation 
is  more  localized  in  comparison  to  composites 
with  uniformly  distributed  reinforcement,  and 
deformation  in  the  cluster  is  more  severe  in  the 
reinforcement  cluster  than  in  the  remainder  of 
the  matrix. 
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STRENGTHENING  OF  COMPOSITES  DUE  TO 
MICROSTRUCTURAL  CHANGES  IN  THE  MATRIX 
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{RcirifcJ  23  December  1989;  in  revised  form  22  June  1990) 

.Abstract  -  The  addition  of  discontinuous  silicon  carbide  (SiC)  to  aluminum  (Al)  alloys  can  result  in  a 
five-fold  increase  in  the  yield  stress.  The  magnitude  of  the  increase  is  obviously  a  function  of  the  volume 
fraction  and  the  particle  size  of  the  SiC.  Previously,  it  was  proposed  that  the  strength  increase  due  to  SiC 
addition  to  Al  alloys  was  the  result  of  change  in  the  matri.x  strength,  i.c.  an  increase  in  dislocation  density 
and  a  reduction  of  subgr:"’'  size.  The  data  obtained  from  a  series  of  experiments  indicate  that  dislocation 
density  increases  w  ith  an  n. crease  in  volume  fraction  of  SiC  and  decreases  with  an  increase  in  particle  size. 
The  suhgrain  size  decreases  as  the  volume  fraction  increases  and  increases  as  the  particle  size  increases 
There  is  a  good  correlation  between  the  microstructural  changes  in  the  matrix  and  the  changes  in  the  yield 
stress  of  the  composites. 


R«ume  L'addition  de  carbure  de  silicium  (SiC)  di.scontinu  a  des  alliages  d'aluminium  ( Al)  pent  conduirc 
a  une  elevation  d  un  facteur  5  de  la  limite  elaslique.  La  grandeur  de  cette  augmentation  est  evidemment 
fonction  de  la  fraction  volumiquecl  de  la  tailledcs  particules  de  SiC.  On  a  deja  suggerc  que  faugmentalion 
de  la  resistanee  me’canique  due  ;i  l'addition  de  SiC  dans  les  alliages  d'aluminium  est  le  rcsultat  d'un 
changemeni  de  la  resistance  mecamquc  de  la  matricc.  e'est  a  dire  d  une  augmentation  de  la  densitc  de 
dislocations  et  d  une  reduction  de  la  taille  des  sous-grains.  Les  donnees  obtenues  pour  une  serie 
d'experiences  indiqucnt  que  la  densitc  de  dislocations  croit  lorsque  la  fraction  volumique  de  SiC  augmente 
et  devroit  lorsque  la  taille  des  particules  croit.  La  taille  des  sous-grains  decroit  lorsque  la  fraction 
volumique  augmente  et  croit  lorsque  la  taille  des  particules  diminue.  II  existe  une  bonne  correlation  entre 
les  changements  microstructuraux  dans  la  matrice  el  les  variations  de  la  limite  daslique  des  composites. 

/.usammenfas.sung  Zugabe  von  diskoniinuierlichem  Siliziumkarbid  (SiC)  zu  Aluminium  (.Al)-Legierun- 
gen  kann  die  Flietispannung  verfunfi'achen.  Diese  Erhohung  hiingt  offenkundig  vom  Volumanliel  und  der 
TeilchengrdBe  des  SiC'  ab  Fruher  wurde  dargclegl.  daB  diese  Festigkeitserhohung  von  einer  Anderung 
in  der  Matrixfestigkeit.  d  h.  eineni  Anstieg  in  der  Versetzungsdichte  und  einer  Verringerung  der 
SubkorngrdBe,  herriihrl  Fine  Rcihe  von  crhaltenen  cxperimentallen  F.rgebnis,sen  zeigt.  daB  die  Verset¬ 
zungsdichte  mil  zunehmenden  Volumanliel  von  SiC  ansteigi  und  mit  zunehmendcr  Teilchengr6/(e 
abnimmt.  Die  SubkvirngrdUe  ninimi  mit  zunehmendem  Volumanicil  ab  und  mit  zunehmendcr  Teilchen 
groBc  zu  Zwischen  den  Anderungen  in  der  Mikrosiruklur  der  Matrix  und  den  Anderungen  ier 
FlieBspannung  des  Materials  besleht  gulc  Korrelation. 


1.  INTRODI  C  TION 

The  concept  of  strengthening  due  to  an  increase  in 
dislocation  density  has  been  in  existence  for  several 
decades,  probably  since  shortly  after  the  original 
proposal  of  edge  dislocations  by  Orwan  et  at.  [I|. 
Over  the  subsequent  years,  there  have  been  numerous 
investigations  whieh  have  correlated  a  strength  in¬ 
crease  with  an  increase  in  dislocation  density.  The 
strengthening  due  to  a  reduction  in  subgrain  size 
is  not  as  old  as  the  concept  of  strengthening  due  to 
an  increa.se  in  the  dislocation  density  (2).  However, 
the  correlation  has  been  demonstrated  quite  clearly 
by  several  individuals;  in  particular,  by  McQueen 
et  al.  [3]. 

It  had  been  proposed  by  Arsenault  [4]  a  few  years 
ago  that  the  strengthening  of  Al  and  Al  alloys  due  to 
the  addition  of  discontinuous  SiC  reinforcements 


occurred  as  a  result  of  generation  of  dislocations 
which  produced  an  mcrea.se  in  dislocation  density 
and  a  reduced  subgrain  size  in  the  matrix.  The 
generation  is  due  to  a  difference  in  thermal  coefficient 
of  expansion  between  the  matrix  and  the  reinforce¬ 
ment  (ACTE).  The  presence  of  the  high  dislocation 
density  within  the  matri.x  has  been  demonstrated 
along  with  a  very  .small  subgrain  size.  .An  in  situ 
investigation  has  demonstrated  the  validity  of  the 
ACTE  mechanisms  as  a  means  by  which  the  high 
dislocation  density  is  produced  within  the  matrix  [5]. 

Further,  it  has  been  shown  that  a  very  simple 
model  based  on  prismatic  punching  is  capable  of 
producing  an  adequate  dislocation  density  [6].  How¬ 
ever,  it  should  be  clearly  pointed  out  that  this  model 
will  only  predict  the  lower  limit  of  the  dislocation 
density.  Attempts  to  calculate  the  upper  limit  of  the 
dislocation  density  due  to  the  ACTE  effect  have 
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proven  unsuccessful,  for  all  of  the  calculations  pre¬ 
dicted  a  very  high  dislocation  density,  i.e.  approach¬ 
ing  infinity,  which  obviously  is  not  correct. 

However,  there  remain  several  questions  pertaining 
to  the  details  of  the  strengthening  mechanism  of 
discontinuously  reinforced  metal  matrix  composites 
based  on  the  strengthening  of  the  matrix  due  to  a 
change  in  the  microstructure  of  the  matrix. 

The  first  question  is  related  to  transmission  elec¬ 
tron  microscopy  (TEM)  procedures  for  determining 
dislocation  densities;  (a)  does  the  sample  preparation 
technique,  i.e.  dimpling  and  ion  milling,  introduce 
damage  into  the  TEM  foil?  and  (b)  if  a  thin  region 
of  the  foil  is  examined  (i.e.  the  viewing  region)  at 
100-200  kV,  will  a  lower  dislocation  density  be  ob¬ 
served?  In  regard  to  (b)  of  the  above  question,  there 
have  been  a  few  investigations  which  have  determined 
the  dislocation  density  as  a  function  of  foil  thickness 
and  the  most  recent  are  by  Fajita  et  al.  [7-9].  In  the 
case  of  pure  Al,  it  was  found  that  the  foil  thickness 
had  to  be  greater  than  1  /tm  to  obtain  a  bulk 
dislocation  density.  The  question  arises,  since  SiC/Al 
is  not  high  purity  Al,  does  this  requirement  of  foil 
thickness  still  apply? 

The  second  question;  is  there  an  overall  relation¬ 
ship  between  dislocation  density  (p),  subgrain  size 
(/),  reinforcement  particle  size  (D)  and  volume  frac¬ 
tion  (F),  and  the  increase  in  the  yield  stress  of  the 
composite?  There  have  been  several  investigations 
where  the  dislocation  density  was  determined  for  a 
specific  particle  size  and  volume  fraction,  but  there 
have  been  no  systematic  investigations  of  these  par¬ 
ameters.  Also,  there  are  no  systematic  investigations 
of  the  subgrain  size  as  a  function  of  particle  size  and 
volume  fraction. 

Another  question:  is  it  possible  to  increase  the 
strength  of  the  matrix  sufficiently  to  account  for  the 
strengthening?  In  other  words,  if  the  Al  matrix  alone 
(with  no  Sic  present)  had  the  same  dislocation 
density  as  in  the  composite,  would  the  Al  matrix  have 
the  same  strength  as  the  composite?  Is  it  possible  to 
increase  the  dislocation  density  and  reduce  the  sub¬ 
grain  size  within  the  AI  alloy  matrix  by  cold  working 
to  approach  the  yield  strength  of  a  composite  con¬ 
taining  20  F%  reinforcement  in  the  same  Al  alloy 
matrix. 

An  investigation  was  undertaken  to  address  the 
questions  raised  above. 

2.  MATERIALS  AND  TESTING 

A  discussion  of  materials  and  testing  will  be  di¬ 
vided  into  three  portions  corresponding  to: 

2.1.  Composite  materials; 

2.2.  Cold  rolled  0  F%  1100  Al  alloy; 

2.3.  High  purity  Al. 


tNow,  Advanced  Composite  Materials  of  Greer,  South 
Carolina. 

{Composite  Specialties  of  Chatsworth,  California. 


2.1.  Composite  materials 

All  of  the  composite  materials  were  produced  by  a 
powder  metallurgical  procedure.  The  Al  alloy  pow¬ 
ders  were  mixed  with  particulate  or  whisker  SiC, 
hot-compacted  and  hot-extruded  into  12.7  mm  diam¬ 
eter  rods  or  hot  pressed  into  a  plate.  The  0,  5  and 
20  V%  SiC  whiskers/6061  Al  alloy  composites  plus 
the  0  F%/1100A1  alloy  composites  were  purchased 
from  ARCO  Silagt.  The  20  F%  SiC-250  pm  particu¬ 
late/1 100  Al  alloy  composite  was  purchased  from 
DWAt  in  the  form  of  a  plate.  The  remainder  of  the 
material  was  produced  in  the  Metallurgical  Materials 
Laboratory  of  the  University  of  Maryland  in  the 
form  of  an  extruded  rod. 

From  the  extruded  rod  and  plate,  tensile  samples 
were  machined  into  a  configuration  as  described 
elsewhere  [4],  Prior  to  tensile  testing,  the  samples  were 
annealed  for  12  h  at  803  K  and  furnace  cooled.  All  of 
the  tensile  testing  was  conducted  at  room  tempera¬ 
ture  at  a  strain  rate  of  10  ‘‘s"'. 

Short  sections  of  extruded  rod  (12  mm)  and  por¬ 
tions  of  as-pressed  plate  were  annealed  at  803  K  for 
12  h  and  furnace  cooled.  From  these  sections  and 
portions,  0.5  mm  thick  disks  were  electrically  dis¬ 
charge  machined.  The  procedure  for  thinning  the 
TEM  foils  and  the  method  of  determining  the  dislo¬ 
cation  densities  is  given  elsewhere  [5].  The  electron 
microscopy  investigation  involved  the  use  of  two 
microscopes,  a  JEOL  100  CX  operating  at  100  kV 
and  the  Argonne  National  Laboratory  modified 
Kratos/AEI  EM7  HVEM  operated  at  I  MV. 

2.2.  Cold  rolled  0  V%  1100  Al  alloy 

The  cold  rolled  strips  were  produced  by  taking 
0  F%  1 100  Al  alloy  material  purchased  from  ARCO 
Silag  and  rolling  without  any  intermediate  annealing 
until  a  reduction  in  thickness  of  ~  90?'o  was  obtained. 
From  the  cold  rolled  strips,  “dog  bone"  shaped 
tensile  samples  were  machined  with  a  24.5  mm  gage 
length  and  a  width  of  2.6  mm  and  a  thickness  of 
1.42  mm. 

The  tensile  sample  preparation  and  testing  were  the 
same  as  described  in  Section  2. 1 .  This  also  applies  to 
the  TEM  foil  preparation  and  examination. 

2.3.  High  purity  A I 

Relatively  high  purity  Al  (99.99%)  was  obtained  in 
the  form  of  a  rod.  12.7  mm  dia  from  ALCOA  as  a 
gift.  This  material  was  annealed  for  12  h  at  803  K  and 
furnace  cooled.  The  TEM  foils  of  high  purity  Al  were 
produced  by  two  different  procedures.  Some  of  the 
foils  were  produced  using  the  electrical  discharge 
machining,  dimpling  and  ion  milling  procedure,  using 
the  same  operating  parameters  as  used  for  the  com¬ 
posite  TEM  foils  (Section  2.1).  The  other  samples, 
after  being  cut  into  0.5  mm  thick  disks  by  electrical 
discharge  machining,  were  electrochemically  thinned 
using  a  jet  polishing  device.  The  foils  were  examined 
at  100  kV  and  1  MV. 
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3.  RESULTS 

The  discussion  of  the  results  will  be  divided  into 
three  portions: 

3.1.  Yield  strength  vs  particle  size  and  yield 
strength  vs  volume  fraction; 

3.2.  Dislocation  density  determination; 

3.2.1.  Dislocation  density  due  to  foil  prep¬ 
aration  procedures  and  TEM  operating 
voltage; 

3.2.2.  Dislocation  density  and  subgrain  size  vs 
particle  size  and  volume  fracture; 

3.3.  Changes  in  dislocation  density  and  sub¬ 
grain  size  due  to  cold  rolling  and  the 
determination  of  the  yield  strength. 

Yield  strength 

Of  interest  is  the  incremental  increase  in  the  yield 
stress  (which  is  defined  as  the  stress  at  a  plastic  strain 
of  0.2%)  due  to  the  addition  of  the  reinforcement. 
Therefore,  the  yield  stress  of  the  0  K%  matrix  alloy 
is  subtracted  from  the  yield  stress  of  the  composite. 
The  incremental  yield  stress  (Aff,)  is  what  is  plotted 
and  listed  in  all  cases.  The  variations  in  of  the  SiC 
whisker  6061  A1  alloy  composites  were  considerably 
less  than  that  in  the  SiC  1100  A1  alloy  composites. 
Th'S  is  simply  due  to  the  greater  homogeneity  of  the 
reinforcement  within  the  matrix  in  the  6061  A1  alloy 
composites. 

As  expected,  as  the  size  of  the  SiC  particulate 
increased,  the  strength  decreased,  as  shown  in  Fig.  1. 
The  squares  represent  the  average  of  at  least  three 
tests  and  there  could  be  as  many  as  ten  tests  for  a 
given  particular  size.  The  arrow  bars  define  the  range 
of  scatter.  The  ultimate  strength  is  only  slightly  larger 
than  the  yield  strength  indicating  that  the  work 
hardening  rate  past  yielding  is  small.  The  uniform 
strain  was  fairly  constant,  between  3  and  4%  for 
the  0.5-9  um  composite  materials,  then  increasing 
to  6-7%  uniform  strain  for  the  20  and  70 /rm 
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Fig.  I .  The  change  in  the  yield  stress  as  a  function  of  SiC 
particulate  size.  The  matrix  is  a  IIOOAI  alloy.  The  com¬ 
posite  is  in  the  annealed  condition. 
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Fig.  2.  The  change  in  the  yield  stress  as  a  function  of  the 
volume  fraction  of  SiC  (SiC  whisker  is  ^0.5  /jm  in  dia  (and 
on  average  ~ 2  |im  long).  The  matrix  is  a  6061  A1  alloy.  The 
composite  is  in  the  annealed  condition. 

particulate  and  slightly  larger  on  average  for  the 
250 /tm  particulate  size  composites. 

Increasing  the  volume  fraction  of  SiC  resulted  in  an 
increase  in  strength  (Fig.  2).  Experimentally,  it  was 
observed  that  Aoj  increases  to  100  MPa  due  to  the 
addition  of  20  V%  SiC  whisker.  The  classical-contin¬ 
uum-load  transfer-shear  lag  model  [10]  would  predict 
~  15  MPa  increase  in  strength,  whieh  is  indicated  by 
the  dashed  line  in  Fig.  2. 

3.2.  Dislocation  density  determinations 

3.2.1.  Dislocation  density  due  to  foil  preparation 
procedures  and  TEM  operating  voltage.  The  two 
different  series  of  investigations  were  undertaken  to 
determine  if  the  TEM  foil  preparation  procedure  was 
introducing  damage  into  the  composite.  The  first 
investigation  involved  preparing  TEM  foils  from 
annealed  and  furnace  cooled  high  purity  A1  rods  by 
two  different  routes.  The  first  route  involved  electrical 
discharge  machining,  followed  by  dimpling  using 
0.25 /im  diamond  paste  and  then  followed  by  ion 
milling,  i.e.  the  same  procedure  which  is  used  for  the 
composites.  The  second  technique  involved  electrical 
discharge  machining,  followed  by  electro-chemically 
thinning.  Figure  3(a)  is  a  micrograph  taken  of  an 
ion  milled  thinned  sample  and  Fig.  3(b)  is  taken 
from  a  sample  which  was  electrochemically  thinned. 
As  can  be  seen,  the  dislocation  densities  are  low  in 
both  cases,  and  the  dislocation  configurations  are  the 
same.  The  only  recognizable  difference  between 
the  ion  milled  and  the  electro-chemically  thinning  is 
the  “mottled"  appearance  of  the  ion  milled  micro¬ 
graph.  This  mottled  appearance  is  the  result  of  the 
non-uniform  thickness  of  the  ion  thinned  foil;  the 
reason  for  the  non-uniformity  is  not  known. 
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Fig.  3.  (a)  A  transmission  electron  micrograph  of  high  purity  Al  in  the  annealed  condition  thinned  by 
dimpling  and  ion  milling,  (b)  A  transmission  electron  micrograph  of  high  purity  Al  in  the  annealed 
condition  thinned  by  dimpling  and  jet  polished. 


However,  it  is  safe  to  conclude  that  the  dimpling  The  second  investigation  was  related  to  foil  thick- 
and  ion  milling  do  not  produce  any  detectable  in-  ness  (viewing  thickness)  as  related  to  dislocation 
crease  in  the  dislocation  density.  density.  If  the  operating  voltage  of  the  TEM  is 
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100  kV,  the  thickness  of  the  foi)  in  which  dislocations 
may  be  viewed  is  0.2-0.4/im.  However,  if  the  operat¬ 
ing  voltage  is  I  MV,  the  thickness  of  the  foil  increases 
to  1-2  If  we  consider  the  micrograph  shown  in 
Fig.  4,  we  can  see  the  effect  of  foil  thickness  on  the 
dislocation  densities.  Figures  4(a)  and  (b)  are  micro¬ 
graphs  taken  at  100  kV  and  the  foil  thickness  in 
Fig.  4(b)  in  which  a  few  dislocations  are  visible  is 
0.3 /tm  and  the  dislocation  density  is  1.6  x  10'^  m  ’. 


Figure  4(c)  and  (d)  are  micrographs  taken  at  an 
operating  voltage  of  1  MV.  The  dislocation  density  in 
Fig.  4(d)  is  4.5  x  10”  m  ’  and  the  foil  thickness  is 
0.9  /im.  Figures  4(a),  (b),  (c)  and  (d)  are  micrographs 
taken  from  the  same  area  of  a  foil  under  the  same 
operating  diffraction  condition,  only  a  thinner  area 
around  the  rim  of  the  hole  is  transparent  for  the 
100  kV  microscope.  This  is  a  general  result  for 
all  composites  examined.  If  a  low  operating  voltage 


Fig.  4.(a)  Caption  otvrleaf. 


Fig.  4.(b,c)  Caption  on  facinp  page. 
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Fig.  4.  (a)  A  transmission  electron  micrograph  taken  at  an  operating  voltage  of  100  kV.  (b)  A  transmission 
electron  micrograph  taken  at  an  operating  voltage  of  100  kV.  This  micrograph  is  of  the  enclosed  area  in 
(a),  (cl  A  transmission  electron  micrograph  taken  at  an  operating  voltage  of  I  MV.  (d)  A  transmission 
electron  micrograph  taken  at  an  operating  voltage  of  1  MV.  This  micrograph  is  of  the  enclosed  area  in  (c). 


is  used,  i.e.  100“200kV,  the  apparent  dislocation 
density  is  low  compared  to  the  dislocation  density 
obtained  from  foils  examined  at  1  MV. 

A  factor  which  has  to  be  recognized  when  examin¬ 
ing  foils  at  high  accelerating  voltages  is  the  possibtlity 
of  radiation  damage  resulting  in  the  formation  of 
dislocation  loops.  It  was  observed  that  if  a  given  area 
was  observed  for  ^  i  h  at  1  MV  and  the  dislocation 
density  was  low  (<1  x  10"  m  ■).  then  dislocation 
loops  became  visible  from  the  condensation  of  the 
radiation  damage.  These  dislcKation  loops  were  very 
small  (<0.1  /I ml  and  could  be  readily  distinguished 
from  the  dislocation  produced  by  the  ACTE  and 
were  never  counted  when  determining  the  dislocation 
density. 

Therefore,  in  order  to  obtain  realistic  measures  of 
dislocation  density  in  SiC  A1  composites,  it  is  necess¬ 
ary  to  examine  thicker  foils  which  needs  a  operating 
voltage  (800 kV  to  1  MV), 

3.2.2.  Di.sloi  alien  dcn.'iily  and  suhgrain  she  vs  par- 
lute  she  and  volume  fravtion.  The  determination  of 
dislocation  density  by  TEM  is  fraught  with  many 
perils.  Two  of  the  many  possibilities  have  already 
been  addressed,  i.c.  introduction  of  dislocations 
during  sample  preparation  and  the  effect  of  the  foil 
thickness  on  the  observed  density.  There  arc  certainly 
several  others,  especially  if  numerous  slip  systems  are 
operative  as  they  certainly  are  in  composites.  In  a 
two-beam  condition,  this  means  that  several  operat¬ 
ive  -vectors  have  to  be  considered,  but  when  operat¬ 
ing  at  I  MV,  a  many-beamed  condition  can  easily  be 
used  and  a  larger  fraction  of  the  dislocations  present 
become  visible.  However,  it  is  still  necessary  to  tilt  the 
sample  and  take  several  different  pictures.  Also,  there 
is  an  inherent  variation  of  dislocation  density  from 
place  to  place  within  the  matrix.  For  example,  the 
dislocation  density  is  higher  near  the  particles  and 
decreases  with  distance  away  from  the  particle.  This 
variation  is  especially  critical  in  large  particle  size  and 
low  volume  fraction  samples.  In  order  to  overcome 


this  difficulty,  it  is  necessary  to  make  measurements 
on  many  micrographs.  The  procedure  for  measuring 
the  dislocation  density  is  given  elsewhere  [5].  The  data 
given  for  each  condition  (i.e,  particle  size  or  volume 
fraction)  in  Figs  5  and  6  represent  tens  to  hundreds 


Fig  5.  The  change  in  dislocation  density  as  a  function  of 
volume  fraction  of  SiC  (SiC  whiskers  have  a  length  to 
diameter  on  average  of  -2  and  the  dia  '0.5 /im). 


Fig.  6.  The  change  in  dislocation  density  as  a  function  of 
particulate  size. 
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Table  1.  Comparison  of  9 particle  size  2(1  T"©  SiC  Al  compt>Mte 
with  cold  roiled  0  I  ®©  I1(X)A1  compt^iie 


Dislocation 

Subgrain 

A  Experiment 

densitv 

size 

yield  stress 

Material 

( fi  X  10''* m') 

(A  X  ^m) 

(Aa,,,  ...  MPa) 

20  y%  SiCp 

I  I 

2 

H8 

0  r-*©  CR 

0.15 

2 

116 

A  Experimental  yield  stress  =  the  yield  stress  of  the  composite  minus 
the  yield  stress  of  0  r**©  matrix. 


of  micrographs  of  a  given  foil  and  5-6  different  foils 
were  examined  for  each  condition.  The  dislocation 
densities  reported  are  the  average  dislocation  den¬ 
sities  in  the  matrix  of  annealed  samples.  As  stated 
above,  the  dislocation  density  is  not  uniform  in  the 
matrix. 

Upon  increasing  the  volume  fraction,  the  dislo¬ 
cation  density  increases,  but  upon  increasing  the 
volume  fraction,  the  subgrain  size  decreases  as  shown 
in  Fig.  5.  and  upon  increasing  the  particle  size,  the 
opposite  is  the  case,  the  dislocation  density  decreases 
(and  the  subgrain  size  increases)  as  the  particle  size 
increases  (Fig.  6). 

3.3.  Effect  of  cold  rolling 

The  yield  stiength  and  the  dislocation  density  were 
determined  of  a  0  F“-o  1 100  Al  alloy  which  had  been 
cold  rolled  ~90%.  The  cold  rolling  reduced  the 
subgrain  size  slightly  from  5  lO/im  to  2/(m.  but 
there  was  a  significant  increase  in  dislocation  from 
less  than  10'"  m  'to  1.5  x  10"  m  ^  The  cold  rolling 
increased  the  yield  strength  from  .^0  to  146  MPa.  If  a 
comparison  is  made  with  the  9  ftm  SiC  particle  size 
in  20  I'  %  llOOAl  alloy  matrix  composite,  it  can  be 
seen  that  the  Aa,  is  greater  for  the  cold  rolling 
material  than  that  of  the  SiCp  llOOAl  alloy  com¬ 
posite  (Table  1).  Flowever.  the  interesting  point  is 
that  the  dislocation  density  in  the  composite  is  much 
greater  than  that  of  the  0  r”o  ll(K)AI  alloy  cold 
rolled  material. 

4.  DISCUSSION 

The  discussion  will  follow  along  the  lines  of  the 
questions  raised  in  the  Introduction 

4.1.  Dislocation  and  suhgrain  size  determinations 

From  the  investigation  of  high  purity  Al.  it  was 
shown  that  the  dislocation  density  in  the  foils,  pro¬ 
duced  by  dimpling  and  ion  milling  was  the  same  as 
that  in  the  electrochemically  produced  foils.  A  main 
concern  of  ion  milling  is  that  the  heat  generation 
occurs  and  that  this  heat  generation  would  result  in 
a  decrease  in  the  dislocation  density.  In  terms  of 
dimpling,  which  involves  grinding  the  foil  with 
0.25  ftm  diamond  paste,  this  could  plastically  deform 
the  sample,  i.e.  introduce  dislocations.  These  dislo¬ 
cations  would  have  an  entirely  different  appearance 
as  compared  to  the  remnant  dislocations  which  re¬ 
main  from  the  long  term  high  temperature  annealing. 
The  outcome  of  this  investigation  is  that  the  dimpling 


and  ion  milling  do  not  introduce  any  dislocations  into 
the  TEM  foil. 

In  regard  to  the  question  raised  concerning  the 
viewing  thickness  of  the  TEM  foil,  it  has  been 
previously  shown  for  the  case  of  high  purity  Al  that 
the  viewing  thickness  of  the  TEM  foil  has  to  be 
greater  than  I  /im  in  order  to  obtain  a  dislocation 
density  representative  of  bulk  samples  [7-9].  High 
purity  Al  is  a  homogeneous  material  and  the  main 
factor  contributing  to  lowering  dislocation  density  is 
the  image  force  acting  on  the  dislocations  within  the 
foil.  In  thinner  foils,  the  presence  of  the  image  force 
should  remove  a  proportionately  larger  number  of 
dislocations.  In  the  case  of  SiC.  Al  alloy  composites, 
we  arc  not  dealing  with  a  homogeneous  material. 
Besides  the  image  force  which  is  still  present,  there  are 
other  factors  which  contribute  to  the  loss  of  dislo¬ 
cations  from  the  thinned  foil.  There  are  thermal 
residual  stresses  within  the  SiC  Al  alloy  composite 
prior  to  thinning.  If  we  consider  the  points  labeled  1 
and  2  in  Fig.  7.  there  is  a  tensile  residual  stress  acting 
in  the  direction  .4  .4  '.  B -B\  C-C  and  D  -D  '.  If  the 
sample  is  cut  along  the  lines  A -A'  and  B-B'.  the 
stress  at  the  surlaces  will  have  to  relax  to  zero.  A 
simple  manner  of  stress  relaxation  is  by  dislocation 
motion,  i.e.  dislocation  motion  out  of  the  foil.  For  a 
thinner  foil,  there  will  be  a  larger  proportional  dislo¬ 
cation  loss  from  the  sample.  The  image  force  is  still 
present  and  would  lead  to  a  further  reduction  of  the 
dislocation  density  in  the  thinner  area  of  the  foil,  just 
as  in  the  homogeneous  material.  Therefore,  in  order 
to  obtain  reasonable  measures  of  the  dislocation 
density  of  a  SiC  Al  alloy  composite,  it  is  necessary  to 
examine  thicker  portions  of  the  TEM  foil,  which 
means  examining  the  foil  at  1  MEV.  The  determi¬ 
nation  of  the  subgrain  size  is  not  as  dependent  on  foil 
thickness  or  the  foil  preparation  procedure. 

4  2.  Correlation  between  microstructurai  changes  and 
yield  stress 

The  increase  in  the  yield  stress  with  an  increase  in 
the  volume  fraction  or  the  decrease  of  yield  stress 
with  an  increase  in  the  particle  size  is  to  be  expected. 
Correspondingly,  the  dislocation  density  and  the 
subgrain  change  in  the  exfiected  manner.  However,  it 
is  necessary  to  consider  if  the  changes  in  the  yield 
stress  can  he  related  to  the  changes  in  dislocation 
density  and  the  .suhgrain  size. 

Before  considering  these  correlations,  it  is  necess¬ 
ary  to  point  out  again  that  the  dislocation  density 
within  the  matrix  is  not  uniform,  and  secondly,  there 
can  be  a  significant  difference  in  stress  at  the  pro¬ 
portional  limit  and  the  yield  stress  defined  at  a  plastic 
strain  of  0.2%.  This  would  indicate  that  there  is 
significant  work  hardening  on  a  local  scale  occurring 
within  the  matrix.  The  reason  for  stating  that  it  is  on 
a  local  scale  is  due  to  the  fact  that  very  little  plastic 
strain  is  occurring  between  the  proportional  limit  and 
the  defined  yield  stress.  It  has  been  observed  that 
the  initial,  and  in  some  cases,  the  entire  plastic 
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Fig  7,  A  schcmalic  arrangemcni  of  SiC  whiskers  and  plaieJets  with  arbitrary  thinning  to  a  TEM  fojJ. 
Tahk  2  Predicted  increase  in  yield  strength  of  20  f'®'o  SiC  IlOO  matrix  and 


experimentally  measured  values 


Panicle  si/c 

0  (;im) 

Aff. 

(MPa) 

Affs<', 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

n  5 

126.6 

131 

17 

240,6 

15.^ 

9 

99 

69 

17 

ISI 

88 

70 

79 

20 

J7 

82 

69 

250 

59  7 

Mt.S 

17 

42.7-50.5 

26 

‘  stress  defined  at  0.2%  plastic  strain  minus  the  0.2%  yield  stress 

of  the  0  r^o  I  lOO  matrix,  =  the  predicted  increase  due  to  the  increased 
dislocation  density;  A<?st.  =  the  predicted  increase  due  to  the  reduc*ed  sub- 
gram  si/c;  Aff^.,  =  the  predicted  average  tensile  residual  stress;  =  the 

sum  of  Aff^,  +  A<7,i^,  —  A(Tj^. 

deformalion  occurs  in  regions  containing  higher  than  where  Act,  is  the  increase  in  the  yield  stress  of  the 

average  volume  fraction  of  particulate  or  whisker  composite  over  that  of  the  0  y%  matrix  material, 

[1 1).  The  dislocation  density  is  the  largest  on  average  a  is  a  constant  and  Hansen  [12]  has  shown  that 

in  the  regions  of  higher  than  average  volume  fraction  for  Al,  this  is  1 .25.  n  is  the  shear  modulus  ot  the 

of  particulate  or  whisker,  but,  there  is  not  an  overall  matrix  (2.64  x  10*  MPa),  b  is  the  Burgers  vector 

change  in  the  density  of  dislocations  in  the  sample.  (2.86  x  10  '“m)  and  p  is  the  dislocation  density. 
Therefore,  an  assumption  is  made  that  the  average  Then  we  can  calculate  the  incremental  increase  in 
increase  in  the  dislocation  density  is  small  due  to  the  strength  due  to  the  presence  of  the  increase  in 

the  initial  loading  to  a  plastic  strain  of  0.2%,  and  that  dislocation  density.  This  is  the  increase  which  should 

the  stress  field  which  opposes  dislocation  motion  be  added  to  the  annealed  0  1^%  Al  alloy  matrix  yield 

had  a  strain  of  0.2%  which  is  due  to  the  average  stress.  In  Table  2  are  listed  the  incremental  increases 

dislocation  density  of  the  annealed  undeformed  in  the  strength  due  to  an  increase  in  dislocation 

sample.  Therefore,  we  can  use  the  well-established  density  as  a  function  of  changes  in  particle  size  at  a 

relationship  constant  volume  fraction  of  20%.  The  incremental 

increase  in  strength  is  large  for  the  case  0.5  pm 
particle  size,  and  decreases  as  the  particle  size 


Act,  =  upb  y/p 
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Table  3.  Predicted  increase  in  yield  stress  of  6061  AI  alloy  as  a 
function  of  volume  f^raction  of  SiC  whisker 


Volume 

(y%) 

A<t^ 

(MPaJ 

Ae^ 

(MPa> 

(MPa) 

(Mrik) 

(MPa; 

1 

37.7 

0 

L7 

36.0 

6.9 

5 

56.6 

7.0 

8.6 

55.0 

52 

20 

70 

48.3 

34.5 

83.8 

100 

=  the  yield  stress  defined  at  0.2*o  plastic  strain  minus  the  0.2®o 
yield  stress  of  the  0  f'/o  1100  matrix;  =  the  predicted 
increase  due  to  the  increased  dislocation  density;  ~  the 
predicted  increase  due  to  the  reduced  subgrain  size;  >=  the 
predicted  average  tensile  residual  stress;  Ao.^^  =  the  sum  of 

A<r^  +  AiXsti  —  . 

increases.  In  the  larger  particle  size  range,  the  pre¬ 
dicted  ijicremental  increase  in  strength  due  to  the 
increased  dislocation  density  alone  exceeds  the 
increase  in  the  experimental  yield  stress. 

If  we  consider  the  increase  in  dislocation  density 
with  the  increase  in  volume  fraction,  it  is  obvious  that 
the  strengthening  due  to  the  increase  in  dislocation 
density  alone  is  very  large  (Table  3).  In  this  case,  it 
is  necessary  to  subtract  out  the  dislocation  density  of 
the  0  V%  material,  for  this  is  appreciable  when 
considering  smaller  volume  fractions.  In  a  case  of 
particle  size  increase,  this  subtraction  is  not  necessary 
for  the  dislocation  density  of  the  0  V%  annealed 
llOOAl  alloy  matrix  was  very  low. 

In  terms  of  strengthening  due  to  a  reduction  in 
subgrain  size,  it  is  necessary  to  rely  on  the  data  of 
McQueen  [3].  In  Table  2  are  listed  the  predicted 
increases  in  strength  due  to  a  reduction  in  the  sub- 
grain  size  Aosc-  The  subgrain  size  of  the  IIOOAl 
matrix  is  much  larger  than  the  subgrain  size  of  the 
composite  so  a  subtraction  is  not  necessary.  The 
amount  of  strengthening  due  to  a  reduction  in  sub¬ 
grain  size  in  the  small  particle  size  range  is  greater 
than  that  in  the  large  particle  size  range.  If  we 
consider  the  effect  of  volume  fraction  on  the  change 
in  the  subgrain  size,  in  this  case,  it  is  necessary  to 
subtract  the  increment  of  0  T%.  The  increase  in 
strength  due  to  the  reduction  in  subgrain  size  be¬ 
comes  predominant  at  20  !■'%. 

Although  the  increase  in  the  dislocation  density 
and  the  decrease  in  subgrain  size  are  the  major  factors 
affecting  the  strength  change  of  the  matrix,  there  are 
other  possible  contributions  to  the  strength  change  of 
the  matrix  [13].  These  include  residual  stresses,  differ¬ 
ences  in  texture  between  the  composite  matrix  and 
the  matrix  material  without  reinforcement,  classical 
composite  strengthening  (i.e.  load  transfer)  and  dis¬ 
persion  strengthening. 

The  average  residual  stress  is  the  only  factor  which 
has  to  be  considered;  the  other  factors  are  very  small 
or  zero  [13].  The  residual  stress  in  the  matrix  is  in 
tension  [14]  and  the  magnitude  depends  upon  the 
volume  fraction,  morphology  and  size  of  the  re¬ 
inforcement.  the  type  of  alloy  of  the  matrix,  and  the 
effective  temperature  change.  In  the  calculations  of 
Taya  and  Arsenault  [14]  it  was  assumed  that  the 
thermal  stresses  were  completely  relieved  in  the  tem¬ 
perature  range  from  the  annealing  temperature  down 


to  473  K.  From  473  K  to  room  temperature,  it  was 
assumed  that  there  was  no  plastic  relaxation  of  the 
thermal  stress  and  the  thermal  stress  which  developed 
between  473  K  and  room  temperature  became  the 
residual  stress.  The  residual  stress  increases  almost 
linearly  with  an  increase  in  volume  fraction,  and 
whiskers  result  in  a  larger  residual  stress  than  aligned 
platelets.  In  the  case  of  spherical  particles,  the  average 
residua)  stress  in  the  matrix  is  pure  hydrostatic 
tension  which  would  have  no  effect  on  the  long  range 
dislocation  motion  land  hence,  yielding).  The  analyti¬ 
cal  analysis  of  Arsenault  and  Taya  [14]  indicates  that 
the  size  of  the  reinforcement  would  have  no  effect  on 
the  average  residual  stress,  but  in  actuality,  the 
average  residua)  stress  probably  decreases  as  the  size 
of  the  reinforcement  increases.  This  reduction  is  due 
to  greater  than  predicted  dislocation  generation  be¬ 
cause  there  would  be  less  interference  from  dislo¬ 
cations  generated  by  neighboring  particles.  To  recap, 
the  greater  than  predicted  generation  for  the  large 
panicle  size  does  not  mean  that  '.he  dislocation 
density  is  greater  in  larger  particle  size.  It  means  that 
the  average  dislocation  density,  for  example,  for  the 
250 //m  particulate  SiC  1 100  AI  composite,  would  be 
greater  than  that  predicted  by  the  simple  model  of 
Arsenault  and  Shi  [6].  The  experimental  average 
dislocation  density  is  3xl0‘'m  -.  whereas  the 
predicted  value  is  3  x  10"  m  v 

Since  the  average  residual  stress  is  in  tension,  it 
would  be  a  negative  contribution  to  strengthening 
of  the  composite.  The  values  in  Tables  2  and  3  arc 
from  the  analytical  predictions  of  .Arsenault  and 
Taya  [14]. 

In  general,  the  predicted  increase  in  the  yield 
strength  is  much  larger  than  the  experimentally 
determined  yield  strength.  There  are  at  least  two 
possible  explanations  for  this  discrepancy.  The  first 
possibility  is  that  the  value  taken  for  the  constant  y 
in  equation  ( 1 )  is  loo  large.  The  values  reported  in  the 
literature  for  y  range  from  0.5  to  1.25.  The  second 
possibility  is  that  the  assumption  of  using  an  average 
dislocation  density  is  not  valid  for  the  plastic  strain 
at  0.2%  and  may  be  due  to  the  motion  of  dislocations 
in  the  matrix  where  the  dislocation  density  is  much 
lower  than  the  average  density.  This  certainly  could 
account  for  the  discrepancy  in  the  250  p  m  particle 
size  composite  and  the  1  1%  composite  where  the 
variation  in  the  dislocation  density  is  the  greatest, 

4.3.  Slrenglheiiing  of  the  matrix 

As  discussed  in  the  Introduction,  there  is  a  ques¬ 
tion  as  to  whether  the  matrix  can  be  strengthened  to 
the  required  levels  by  changes  in  the  microstructure, 
i.e.  can  there  be  a  sufficient  increase  in  the  dislocation 
density  or  a  sufficient  decrease  in  the  subgrain  size? 

The  cold  rolled  investigation  has  shown  that  the 
0  r%  1 100  AI  alloy  material  can  have  a  yield  stress 
which  is  comparable  in  strength  to  a  20  F%  SiC  At 
alloy  composite. 


ARSENAULT  ei  at.:  STRENGTHENING  DUE  TO  CHANGE  IN  DISLOCATION  DENSITY 


57 


5.  CONCLUSIONS 

This  investigation  had  as  its  purpose  to  address 
several  questions  concerning  the  observed  strength  of 
discontinuous  SiC/Al  alloys,  especially  related  to  the 
mechanism  proposed  which  is  based  on  the  changes 
in  the  microstructure  of  the  matrix  due  to  differences 
in  thermal  coefficient  of  expansion  between  the 
reinforcement  and  the  matrix. 

From  a  consideration  of  the  data,  several 
conclusions  were  obtained: 

•  The  observed  strength  of  the  discontinuous 
Sic  A1  composite  can  be  accounted  for  in  terms  of  a 
mechanism  based  on  a  change  in  the  microstructure 
of  the  matrix, 

•  This  change  in  the  microstructure,  i.e.  an  in¬ 
creased  dislocation  density  and  a  reduced  subgrain 

"s  ccmpaicJ  to  C  f  ‘'o  matrix  alloy,  is  due  to  the 
difference  in  coefficient  of  expansion  between  the 
reinforcement  (SiC)  and  the  matrix  (Al). 

•  The  procedure  of  TEM  foil  preparation  of 
the  composites  which  involves  electrical  discharge 
machining,  dimpling  and  ion  milling  does  not  intro¬ 
duce  dislocations  into  the  matrix. 

•  In  order  to  obtain  realistic  values  of  the  dislo¬ 
cation  density  of  SiC  Al  alloy  composites,  it  is  necess¬ 
ary  to  view  thick  foils,  i.e.  an  operating  voltage  of 
I  MV. 
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Abstract 

At  least  four  distinct  separation  fracture  modes  (mechanisms)  have  been  proposed  for  the  fracture  of 
SiC-Al  composites;  firstly,  cracking  of  the  SiC  particles  or  whiskers  in  front  of  an  advancing  crack  tip 
and  the  linking  of  these  cracks;  secondly,  SiC  particles  or  whiskers  cracking  and  this  crack  initiating  a 
fracture;  thirdly,  void  nucleation  and  growth  at  the  interfaces  of  SiC  particles  or  whiskers  and  the 
matrix;  fourthly,  fracture  of  the  matrix.  There  is  evidence  which  supports  all  these  mechanisms,  and  con¬ 
versely  there  are  data  which  indicate  that  these  mechanisms  have  no  relationship  to  the  actual  fracture 
process.  It  was  found  that  the  mode  of  fracture  is  through  the  matrix. 


1.  Introduction 

Crack  initiation  and  propagation,  Le.  the 
separation  mechanism,  in  SiC-Al  composites 
have  been  investigated  to  some  extent  over  the 
last  few  years.  Essentially,  there  are  three  differ¬ 
ent  questions  that  have  been  raised  concerning 
the  crack  initiation  and  propagation.  These  are  as 
follows. 

( 1 )  Does  the  fracture  of  the  SiC  particles  con¬ 
tribute  to  the  separation  of  the  composite,  i.e.  are 
they  an  integral  part  of  the  fracture  mode? 

(2)  Does  void  nucleation  at  the  SiC  whisker  or 
particle  and  growth  contribute  to  the  fracture 
process? 

(3)  What  is  the  significance  of  the  crack-tip- 
opening  displacement? 

The  fracture  of  SiC  whiskers  or  particles  can 
conceivably  initiate  fracture  of  the  composite.  If 
the  particle  or  whisker  were  large  enough,  then  a 
crack  of  critical  length  could  form  and  if  the 
dynamic  Kj^  of  the  matrix  were  low  enough,  the 
crack  would  continue  to  propagate  into  the 
matrix,  leading  to  failure  of  the  composite  [1].  If 
the  fracture  of  SiC  could  initiate  fracture,  then  the 
fracture  of  SiC  should  aid  in  the  propagation  of 
the  crack.  It  has  been  shown  that  large  intermetal- 
lic  inclusions  can  act  as  a  fracture  initiation  site  in 


'Deceased. 


SiC-Al  composites  [2].  The  most  recent  and 
probably  the  most  thorough  examination  of  crack 
propagation  in  particulate  SiC-Al  which  sup¬ 
ports  the  concept  that  cracking  of  SiC  particles  in 
front  of  the  end  crack  is  the  significant  mode  of 
separation  in  SiC-Al  composite  was  carried  out 
by  Shang  and  Ritchie  [3].  There  have  also  been 
some  prior  investigations,  e.g.  by  You  et  al.  [4],  in 
which  it  was  concluded  that  the  fracture  of  SiC 
particles  was  important.  In  these  investigations 
they  examined  what  are  defined  as  post-mortem 
samples.  The  crack  has  been  propagated  and 
opposite  sides  of  the  fracture  surface  are 
examined  for  chunks  of  SiC  particles  and,  indeed, 
they  did  see  matching  pairs.  Also,  in  an  investiga¬ 
tion  by  Flom  and  Arsenault  [5]  of  the  fractme 
SiC-Al  composites,  it  is  shown  that  there  are 
large  numbers  of  particles  on  the  fracture  surface, 
especially  in  composites  containing  larger  (20 
fim)  SiC  particles.  However,  if  one  examines 
these  fracture  surfaces  by  scanning  Auger  mi¬ 
croscopy,  then  the  number  of  places  at  which  it  is 
possible  to  obtain  a  clear  and  coherent  SiC  peak 
is  relatively  rare  but  it  should  be  pointed  out  that 
in  some  cases  a  clear  and  distinct  Auger  spectrum 
of  SiC  can  be  obtained  from  SiC  particles  on  the 
fracture  surfaces,  as  shown  in  Fig.  1 . 

Also,  in  an  experiment  to  measure  the  bond 
strength,  a  tensile  sample  of  aluminum  (purity, 
99.99%)  was  fabricated  which  contained  a  SiC 
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Fig.  1.  A  scanning  electron  micrograph  of  the  fracture  sur¬ 
face  of  a  tension  sample  fractured  in  the  scanning  Auger  mi¬ 
croprobe  (sample  of  5  vol.%  5  /tin  SiC  particulate  in  an 
aluminum  alloy  2124  matrix  1. 


sphere  3  mm  in  diameter  in  the  center  of  the 
gauge  section.  The  gauge  section  was  25  mm  long 
and  had  a  diameter  of  6  mm.  When  the  .sample 
was  tested,  shortly  after  macro-yielding,  the  SiC 
sphere  fractured  perpendicular  to  the  tensile  axis 
with  a  small  drop  in  load.  The  test  was  continued 
and  the  sample  finally  failed,  but  the  sample  did 
not  fail  by  extension  of  the  crack  in  the  SiC  into 
the  aluminum  matrix,  even  though  there  was  a 
25%  reduction  in  the  cross-sectional  area  as  a 
result  of  the  fractured  SiC  sphere.  The  sample 
failed  at  another  place  in  the  gauge  length  (6j. 

Another  mechanism  which  has  been  propctsed 
by  numerous  investigators  is  related  to  void 
nucleation  at  SiC  particles  and  subsequent 
growth  of  these  voids  [7)  as  the  .separation 
mechanism.  If  we  consider  Fig.  2  which  is  a  plot 
of  the  logarithm  of  the  initial  cross-sectional  area 
of  the  sample  divided  by  the  area  of  fracture  (as  a 
measure  of  ductility)  vs.  volume  fraction  of  par¬ 
ticulate,  then  it  can  be  seen  that  as  the  volume 
fraction  increases,  the  ductility  decreases.  The 
area  between  the  two  broken  curves  in  Fig.  2  con¬ 
tain  a  great  deal  of  data  which  has  been  generated 
on  numerous  systems  in  which  there  is  a  hard 
particle  (inclusion  or  precipitate)  in  a  relatively 
ductile  matrix.  SiC  particles  in  an  aluminum 
matrix  would  .satisfy  this  criterion,  i.e.  a  hard  par¬ 
ticle  in  a  ductile  matrix  and.  as  can  be  seen,  for 
various  types  and  conditions  of  particulate  and 
whisker  SiC-AI  composites,  the  results  show  the 
same  general  trend,  i.e.,  as  the  volume  fraction 
increases,  the  ductility  decreases. 

Further,  if  we  can  continue  this  discussion,  the 
work  by  Gerberich  18]  and  .several  other  investi- 


Fig.  2.  A  plot  of  the  ductility  rv.  volume  percentage  of  par¬ 
ticulate  in  ductile  matrix  materials: - . - .the  upper 

and  lower  bounds  respectively  of  a  great  deal  of  data  which 
has  been  collected  on  particulate-in-ductile-matrix  material. 

- ,  particulate  and  whisker  SiC  in  aluminum  alloy  bdb  1 

in  the  annealed-and-heat-treated  T6  condition. 


gators  have  shown  that,  if  the  particle  size 
increases  and  the  volume  fraction  is  held  con¬ 
stant.  then  A.',^  should  increase.  Flowever.  experi¬ 
mental  observations  by  Flom  and  Arsenault  [5] 
have  shown  that  A.',,  in  the  small  particle  range. 
i.e.  20  /rm  and  less  in  size,  is  not  dependent  upon 
the  particle  size  or  the  interparticle  spacing.  For  a 
much  larger  particle  size  (250  /<m),  then  the 
mode  of  fracture  changes  dra.stically;  there  is  a 
significant  amount  of  fracture  of  the  SiC  particles. 

Finally,  a  remaining  topic  is  related  to  crack  tip 
opening  and  displacement.  The  question  is:  what 
is  the  crack-tip-opening  displacement,  if  any,  and 
what  does  it  mean?  It  has  been  argued  by  several 
individuals  (3,  7]  that  crack-tip-opening  displace¬ 
ment  is  a  real  value  which  has  real  meaning.  If  it  is 
assumed  that  void  nucleation  at  SiC  particles  and 
whiskers  and  subsequent  growth  of  the  void  are 
the  dominant  mode  of  fracture,  then  the  crack  tip 
opening  should  be  some  large  fraction  of  the 
whisker  or  particle  spacing.  However,  work  on 
post-mortem  samples  indicates  that  the  crack-tip- 
opening  displacement  is  vanishing  small  (Fig.  3). 
but  a  counter-argument  could  be  made  that,  when 
examining  post-mortem  samples,  the  sample  is 
not  under  load  and  therefore  there  may  be  signifi¬ 
cant  closure  of  the  crack. 

Thus  an  investigation  was  undertaken  to  deter¬ 
mine  whether  cracking  of  SiC  particles  did  or  did 
not  occur  during  crack  propagation  of  SiC-Al 
deposits.  Also,  as  part  of  this  investigation,  the 
possible  role  of  void  nucleation  and  growth  at  SiC 
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Fig.  3.  A  scanning  electron  micrograph  of  a  post-mortem 
compact  tension  sample  showing  the  crack  tip. 
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Fig.  4.  A  schematic  view  of  a  compact  tension  sample  also 
showing  the  plane  on  which  the  sample  was  sliced  in  half  for 
examination  of  the  crack  tip  and  subsequent  deformation  in 
the  mini-tensile  stage. 


particles  was  considered.  Finally,  measurements 
of  crack-tip-opening  displacements  were  under¬ 
taken. 

2.  Materials  and  test  procedures 

Particulate  SiC  ranging  in  size  from  2.4  to  20 
//m  was  mixed  with  aluminum  alloy  1100  powder 
(44  fim  in  diameter)  and  this  mixture  was  then 
hot  compacted  and  extruded  into  a  bar  which  was 
6  mm  thick  by  3 1  mm  wide.  The  actual  details  of 
production  have  been  described  elsewhere  [9]. 
The  250  //m  SiC  particulate  composite  was  pur¬ 
chased  from  DWA  Composites  Specialities, 
Chatsworth,  CA.  Also,  commercial  aluminum 
alloy  6061  heat  treated  to  the  T6  condition  was 
rolled  to  a  75%  reduction  in  thickness. 

The  CL<mr'’site  material  and  the  aluminum 
alloy  6061  were  machined  into  compact  tensile 
samples  as  shown  in  Fig.  4  and  tested  in  an 
Instron  testing  machine  until  a  crack  propagated 
into  the  matrix.  After  crack  propagation,  the  sam¬ 
ple  was  sliced  in  half  by  electrical  discharge 
machining.  Each  half  of  the  sliced  compact 
tension  .sample  was  glued  to  an  individual  sample 
holder,  so  that  the  cut  surface  could  be  polished. 
In  other  words,  the  cutting  plane  (Fig.  4)  was  the 
surface  which  was  polished.  The  samples  were 
then  placed  on  a  Vibramai  with  a  Buchler 
Mastermet''*  polishing  compound  for  approxi¬ 
mately  3-4  days  to  obtain  a  polished  surface. 
Attempts  were  made  to  use  diamond  paste  and 
diamond  sawing  to  reduce  the  time  required  for 
polishing  the  samples,  but  these  proved  to  be 
unsatisfactory  for  they  produced  “chunks”  of  SiC 
in  the  matrix  and  also,  in  some  cases,  the 


diamond  paste  actually  became  embedded  into 
the  surface  of  the  sample.  After  polishing,  the 
samples  were  mounted  in  a  mini-tensile  stage. 
The  method  of  loading  the  sample  was  by  means 
of  pins  through  the  holes  in  the  compact  tension 
sample.  The  stage  was  constructed  such  that 
these  pins  move  either  toward  or  away  from  each 
other  by  the  screw  drive  system.  The  displace¬ 
ment  of  the  pins  was  calibrated,  and  one  revolu¬ 
tion  of  the  main  drive  feed  resulted  in  the  pins 
being  displaced  by  1 1  fim.  A  drive  system  was 
constructed  for  this  tensile  stage  with  feed¬ 
throughs  so  that  it  could  be  mounted  in  the 
scanning  electron  microscope. 

The  crack  tip  region  was  then  photographed, 
producing  a  map  for  the  purpose  of  determining 
where  and  which  of  the  SiC  particles  were  frac¬ 
tured.  Then  the  crack  was  opened  at  the  pins  in 
increments. 

Transmi.ssion  electron  microscopy  (TEM)  foils 
were  produced  from  annealed  2.4,  3.2,  8  and  20 
^m  SiC  particulate-AI  composite  in  the  same 
manner  as  described  elsewhere  [10].  The  foils 
were  examined  in  the  high  voltage  electron  mi¬ 
croscope  at  Argonne  National  Laboratory  in  the 
in  situ  tensile  stage.  A  video  tape  was  produced 
during  the  straining  of  the  foil. 

3.  Results 

Composite  samples  containing  various  sizes  of 
SiC  particulates  were  investigated;  however,  only 
a  few  of  them  will  be  considered. 

The  specific  case  to  be  discussed  in  detail  is  a 
composite  containing  20  vol.%  8  /rm  SiC  par- 
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Kig.  5.  A  scanning  electron  transmission  micrograph  of  the 
crack  tip  in  the  unloaded  condition  of  a  20  vol.%  8  /rm  par¬ 
ticulate  SiC-AI  composite. 


Fig.  6.  A  scanning  electron  micrograph  of  the  propagation 
of  a  crack  of  a  20  vol.%  8  /tm  SiC  parliculaie-Al  alloy  1 100 
matrix  composite. 


tides.  Figure  5  shows  the  crack  tip  in  the  pre¬ 
crack  condition,  i.e.  in  the  unloaded  condition. 
The  sample  was  then  loaded  by  displacing  the 
pins  by  22  fim  and  the  crack  propagated  250  /t/m; 
Fig.  6  is  a  micrograph  of  the  crack  tip  area.  (In  all 
subsequent  discussions  of  the  crack  tip,  the  load 
is  applied.)  The  pins  were  displaced  by  another  6 
fivn.  Figure  7  again  shows  the  crack  tip  area.  The 
crack  has  advanced  about  5  fun.  If  Fig.  7  is  con¬ 
sidered  in  detail,  it  can  be  seen  that  the  crack-tip¬ 
opening  displacement  is  a  small  fraction  of  a 
micron,  and  the  load  is  still  being  applied  to  the 
sample.  The  crack  opening  is  extremely  small; 
even  at  a  distance  of  5-10  fitn  behind  the  crack 
tip,  the  crack-opening  displacement  is  less  than  1 
fitn  and  probably  less  than  0.1  ^m.  With  another 
33  /<m  pin  displacement,  the  crack  has  pro¬ 
gressed,  as  shown  in  Fig.  8.  In  this  particular  case, 
the  crack  has  stopped  in  front  of  a  rather  large 


Fig.  7.  A  scanning  electron  micrograph  of  the  propagation 
of  a  crack  of  a  20  vol.%  8  /<m  SiC  particulate- Al  alloy  11 00 
matrix  composite. 


Fig.  8.  A  scanning  electron  micrograph  of  the  propagation 
of  a  crack  of  a  20  vol.%  8  ^m  SiC  particulate-AI  alloy  1 100 
matrix  composite. 


particle.  Continued  loading  of  the  sample  results 
in  no  apparent  change  in  the  crack,  except  that 
the  crack  tip  has  opened  further.  When  the  pins 
are  opened  by  an  additional  22  ^m,  the  crack  dis¬ 
places  upward,  as  shown  in  Fig.  9  and  the  crack 
has  made  a  significant  advance  of  about  90  // m. 

Upon  further  loading  of  the  sample,  as  shown 
in  Fig.  10,  the  crack  tip  has  stopped  somewhere 
in  front  of  a  particle  which  apparently  is  cracked 
in  three  different  locations.  The  crack  proceeds, 
as  shown  in  Fig.  11,  through  one  of  pre-existing 
cracks  in  the  cracked  sample.  Figures  12  and  13 
are  continuations  of  the  crack.  Figure  14  is  of  the 
same  area  as  Figs.  12  and  1 3,  except  after  further 
displacement  of  the  loading  pins.  In  this  case, 
there  is  a  distinct  crack-tip-opening  displacement. 
This  is  a  typical  result  when  the  crack  has  propa- 


Fig.  10,  A  scanning  electron  micrograph  of  the  propagation 
of  a  crack  of  a  20  vol.%  8  nm  SiC  particulate- Al  alloy  1 100 
matrix  composite. 


Fig.  1 1 .  A  scanning  electron  micrograph  of  the  propagation 
of  a  crack  of  a  20  vol.%  8  Aim  SiC  particulate-Al  alloy  1 100 
matrix  composite. 


gated  a  significant  distance  from  its  initial  pre¬ 
cracked  po.sition  in  the  bulk  compact  temsion 
sample. 

An  observation  which  occurs  infrequently  was 
evident  in  the  stimple  which  contained  3.2  nm 
SiC  particulate.  In  this  particular  case,  there  is 
evidence  of  cracking  around  the  existing  sub¬ 
grains  within  the  material  (Fig.  15).  Also  in  this 
particular  sample,  there  happens  to  be  an  inclu¬ 
sion,  presumably  of  AI2O3,  which  cracks  as 
shown  in  Fig.  16. 


Another  result  which  is  not  typical  but  seems 
to  occur  more  frequently  in  the  smaller  particu¬ 
late  size  (in  this  particular  case.  2.4  /tm  SiC)  is  the 
greater  appearance  of  the  separation  of  the 
matrix  about  SiC  particles  in  front  of  the  crack 
tip.  This  is  evident  if  we  consider  Figs.  17  and  18 
of  a  2.4  ^m  SiC  composite.  Also,  Fig.  19  shows  a 
good  example  where  the  crack  has  propagated  a 
significant  distance,  and  the  crack-tip-opening 
displacement  becomes  much  more  pronounced. 
However,  if  you  consider  the  crack-tip-opening 
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Fig.  12.  A  scanning  electron  micrograph  of  the  propagation  Fig.  15.  A  scanning  electron  micrograph  of  a  20  vol.%  3.2 
of  a  crack  of  a  20  vol.%  8  ^tm  SiC  paniculate- Al  alloy  1 100  fim  particulate  SiC-Al  alloy  1 100  matrix, 
matrix  composite. 


Fig.  13.  A  scanning  electron  micrograph  of  the  propagation  Fig.  16.  A  scanning  electron  micrograph  of  a  20  vol.%  3.2 
of  a  crack  of  a  20  vol.%  8  /tm  SiC  particulate- Al  alloy  1 100  ^m  particulate  SiC-Al  alloy  1 100  matrix, 
matrix  composite. 


Fig.  1 7.  A  scanning  electron  micrograph  of  a  20  vol.%  2.4 
Fig.  1 4.  A  scanning  electron  micrograph  of  a  propagation  of  ^m  particulate  SiC  composite, 
a  crack  in  a  20  vol.%  8  ftm  SiC  particulate- Al  alloy  1100 
matrix  composite. 

Figure  20  is  a  micrograph  of  the  unloaded 
displacement,  these  displacements  are  trivial  crack  in  the  pre-crack  compact  tension  sample  of 

compared  with  the  crack-tip-opening  displace-  cold-rolled  aluminum  alloy  6061.  Figure  21  is 

ments  which  occur  in  the  cold-rolled  aluminum  after  the  sample  has  been  loaded  to  a  pin-opening 
alloy  6061.  displacement  of  200  fim,  the  crack  tip  has  not 
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Fig.  1 8.  A  scanning  electron  micrograph  of  a  20  vol.%  2.4  Fig.  20.  The  crack  lip  propagation  in  cold-rolled  aluminum 

urn  particulate  SiC  composite.  alloy  606 1  heat  treated  to  the  T 6  condition. 


Fig.  19.  A  scanning  electron  micrograph  of  a  2(1  vol.‘’o  2.4 
fivix  paniculate  SiC  composite. 


Fig.  21.  The  crack  tip  propagation  in  cold-rolled  aluminum 
alloy  6061  heal  treated  to  the  T6  condition. 


progressed  at  all,  but  now  there  is  measurable 
( 1-2  /fm)  crack-tip-opening  displacement  right  at 
the  crack  tip  and  at  about  2-10  ^m  behind  the 
crack  tip.  This  crack-tip-opening  displacement  is 
significantly  different  from  that  of  the  composite 
material  which  contains  a  small  SiC  particle  size. 
Also,  we  observed  measurable  voids  associated 
with  the  inclusions  which  occur  in  this  particular 
material  which  is  not  evident  in  the  composite 
material.  There  were  no  similar  observations  in 
the  SiC-AI  composites,  i.e.  there  is  no  evidence 
of  void  formation  at  SiC  particles  in  the  front  of 
the  advancing  crack. 

The  in  situ  TEM  straining  of  the  SiC-AI  com¬ 
posites  resulted  in  the  formation  of  a  crack  at  the 
edge  of  the  hole  and  then  crack  propagation  into 
the  foil.  In  most  cases,  the  propagation  was  very 
similar  to  that  observed  by  scanning  electron 
microscopy  (SEM),  i.e.  the  crack  propagated 
around  the  SiC  particles.  However,  in  some  ca.ses 
an  opening,  i.e.  a  void,  forms  near  an  SiC  particle 


leaving  a  distinct  aluminum  layer  in  the  SiC  par¬ 
ticle.  The  location  of  the  void  was  predictecJ  by 
Wilner[l  1). 

4.  Discussion 

From  a  consideration  of  the  above  experi¬ 
mental  results  it  is  obvious  that  there  are  differ¬ 
ences  from  some  of  the  previously  generated  data 
on  the  fracture  characteristics  of  SiC-Al  com¬ 
posites,  but  these  differences  can  be  readily 
accounted  for.  Also,  a  clearer  explanation  of 
some  of  the  previous  data  is  now  possible. 

The  previous  data  from  matching  fracture  sur¬ 
faces  and  unloaded  propagated  cracks  supported 
the  concept  that  fracture  of  SiC  particles  in  front 
of  the  advancing  crack  played  a  significant  role  in 
the  fracture  process  and  can  be  readily  explained 
or  accounted  for.  The  above-mentioned  previous 
data  are  based  on  the  examination  of  post¬ 
mortem  samples.  The  basic  problem  associated 
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with  examining  post-mortem  samples  is  that  it  is 
impossible  to  know  whether  the  SiC  particle  or 
whisker  was  broken  owing  to  processing  prior  to 
the  actual  fracture  test.  In  the  present  investiga¬ 
tion  there  was  no  case  where  it  could  be  shown 
that  the  advancing  crack  resulted  in  the  fracture 
of  SiC  particles. 

It  was  clearly  observed  that,  as  the  crack 
advances,  it  has  a  distinct  preference  to  proceed 
through  existing  cracked  SiC  particles.  This 
would  be  entirely  logical,  since  with  the  crack 
advancing  through  existing  cracked  SiC  particles 
there  is  no  need  to  create  two  new  surfaces  since 
they  already  exist.  Therefore  matching  fracture 
surfaces  would  have  a  matching  “fractured”  SiC 
particle.  Also,  the  observation  by  Flom  and  Arse¬ 
nault  [5]  that  the  number  of  SiC  particles  along 
the  fracture  path  is  greater  than  a  random  line 
through  the  composites  is  simply  a  consequence 
that  the  crack  would  rather  propagate  through  an 
existing  cracked  (fractured)  SiC  particle. 

In  terms  of  the  sharpness  of  the  crack  tip  in  the 
composites,  i.e.  lack  of  crack  tip  blunting,  this  can 
be  traced  back  to  the  lack  of  plastic  deformation 
which  does  not  occur  in  these  composites.  In 
order  to  have  crack  tip  blunting,  it  is  necessary  to 
have  in-plane  plastic  deformation.  If  plastic 
deformation  in  the  crack  tip  region  (which  can 
occur)  is  very  small,  there  will  be  no  blunting. 
Considering  a  very  simple  case,  if  the  crack  was 
initially  very  sharp  and  then  blunted  to  a  crack- 
tip-opening  displacement  of  2  fim,  as  in  the  ca.se 
of  the  cold-rolled  aluminum  alloy  606 1 ,  the  effec¬ 
tive  plastic  strain  required  would  be  15-30%  as 
obtained  by  McMeeking  [12)  in  a  finite  element 
method  investigation  of  crack  tip  blunting.  It  may 
be  argued  that  the  analysis  of  McMeeking  is  not 
valid  since  his  analysis  is  for  a  stationary  crack. 
The  cracks  in  the  SiC-Al  composite  propagate  in 
a  discontinuous  manner.  The  major  component 
of  this  effective  plastic  strain  is  a  tensile  strain 
perpendicular  to  the  crack  plane.  However, 
Davidson  (13)  believes  that  the  tensile  strain  is 
very  small  in  SiC-Al  composites.  The  tensile 
strain  is  small  because  of  the  elastic  constraint  of 
the  high  modulus  (in  comparison  with  aluminum) 
SiC  particulate  present  in  the  composite,  and  the 
fact  that  the  matrix  within  the  composite  has  a 
very  high  dislocation  density,  i.e.  equivalent  to  a 
iughly  cold-worked  condition. 

Even  though  the  aluminum  alloy  6061  was 
cold  rolled  75%,  it  is  still  possible  to  have  signifi¬ 
cant  plastic  strain  occurring  at  the  crack  tip  in  the 


plane  stress  condition.  Prior  to  slicing  the  sample 
in  half,  the  crack  propagates  in  a  nearly  plane 
strain  condition.  As  a  result  of  the  greater  plastic 
strain,  there  is  significant  crack  tip  opening 
(several  orders  of  magnitude)  in  comparison  with 
the  SiC-Al  composites. 

The  lack  of  void  nucleation  and  growth  at  the 
SiC  matrix  interface  is  more  than  likely  due  to 
very  high  bond  strength  [14]  between  SiC  and 
aluminum. 


5.  Conclusions 

From  a  consideration  of  the  experimental 
results  obtained  in  this  investigation,  the  follow¬ 
ing  conclusions  can  be  arrived  at. 

( 1 )  The  crack  appears  to  go  around  intact  SiC 
particles,  leaving  a  distinct  outer  coating  on  the 
SiC  particles  which  is  in  agreement  with  in  situ 
TEM  investigations  of  crack  propagation  in  this 
particular  material. 

(2)  The  observance  of  cracking  of  the  SiC  par¬ 
ticles  in  front  of  the  crack  tip  due  to  the  advance 
of  the  crack  tip  is  a  rare  event.  If  there  are  pre¬ 
existing  fractured  SiC  particles,  the  crack  pro¬ 
ceeds  throu^  the  pre-cracked  SiC  particles. 

(3)  The  existence  of  void  nucleation  and 
growth  in  front  of  the  crack  at  the  SiC  particles 
was  not  observed  in  the  SEM  investigations;  how¬ 
ever,  voids  some  distance  from  the  SiC  particles 
were  observed  in  the  TEM  investigations. 

(4)  The  crack-tip-opening  displacement  at  the 
crack  tip  and  5-15  ;<m  behind  the  crack  tip  as  it 
initially  propagates  in  the  composites  conttiining 
SiC  particles  of  2.4-20  ^m  size  was  extremely 
small,  i.e.  less  than  0.1  f/m. 

(5)  If  a  crack  is  propagated  under  the  same 
conditions  in  cold-rolled  aluminum  alloy  6061- 
T6,  then  there  is  significant  blunting  at  the  crack 
tip,  i.e.  crack-tip-opening  displacement  becomes 
measurable  (1-2  /rm)  and,  in  this  particular  case, 
void  nucleation  occurs  at  the  oxide  inclusion 
which  occurs  in  this  alloy. 
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Analytical  Evaluation  of  the  Thermal  Residual  Stresses 

in  SiC/Al  Composites* 


Ning  SHI  and  Richard  Joseph  ARSENAULT** 

An  analytical  investigation  of  the  difference  in  the  constitutive  behavior  of  SiC  -  A1  composites  under  uniaxial 
compressive  and  tensile  loading  was  carried  out  and  the  associated  changes  of  the  residual  stresses  were  studied. 
It  is  suggested  that  the  observed  asymmetric  response  of  the  constitutive  behavior  is  primarily  due  to  the 
existence  of  residual  stresses  in  the  composite,  and  the  residual  stresses  are  initially  introduced  due  to  the 
differences  in  the  coefficients  of  thermal  expansion  (JCTE),  and  subsequently  changed  after  the  external  load 
is  applied.  A  two-dimensional  finite  element  analysis  of  a  hexagonal  array  of  SiC  whiskers  in  an  A1  matrix  was 
performed,  the  result  of  which  testifies  to  the  appropriateness  of  the  proposed  explanation.  It  is  concluded  that 
the  2-D  finite  element  analysis  is  an  economical  way  of  adequately  reproducing  the  most  prominent  features  of 
the  material’s  constitutive  behavior.  Based  upon  the  theory  of  mechanics  of  composite  material,  simplified 
analytical  models  were  developed  which  can  be  used  to  investigate  the  influences  of  the  whisker  matrix  interface 
shear  load  transfer  and  the  volume  mismatch  on  the  residua!  stresses.  It  is  concluded  that  the  mechanism  of  load 
transfer  in  terms  of  normal  stress  at  the  whisker  tip.  which  is  governed  by  the  volume  (bulk)  mismatch  between 
the  matrix  and  the  whisker,  is  predominantly  responsible  for  altering  the  specific  pattern  of  the  residual  stresses 
under  the  applied  load.  In  contrast,  the  effect  of  shear  load  transfer  at  the  longitudinal  whisker-matrix  interface 
is  only  restricted  to  a  very  small  region  in  the  vicinity  of  the  whisker  and  does  not  significantly  influence  the 
overall  proposed  asymmetric  constitutive  behavior. 
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I.  Introduction 

A  previous  investigation  of  the  deformation 
behavior  of  short  whisker  reinforced  SiC-Al  compos¬ 
ites'"  has  indicated  that  the  constitutive  behavior  of 
these  composites  is  distinctly  different  under  uniaxial 
compressive  deformation  as  compared  to  tensile 
deformation.  It  was  generally  found  that  the  compres¬ 
sive  yield  strength  of  these  composites  was  larger 
than  that  of  the  tensile  yield  strength,  while  the 
Young’s  modulus  for  compression  was  smaller  than 
that  in  tension.  If  monolithic  aluminum  is  initially 
deformed  plastically  in  tension,  then  reverse  deformed 
in  compression,  a  difference  in  stress  is  required  to 
continue  plastic  deformation.  Usually,  this  stress  is 
less  upon  reverse  loading.  This  reduction  in  stress  (or 
the  strain  necessary  to  reach  the  previous  level  of  flow 
stress)  is  defined  as  the  Bauschinger  Effect.  In  the 
case  of  monolithic  aluminum,  the  direction  of  initial 
deformation  (i.  e..  tension  or  compression)  does  not 
have  any  effect  on  the  Bauschinger  Effect.  The  SiC  Al 
composites  exhibit  a  remarkable  characteristic,  i.  e.  the 
flow  stress  drop  (or  alternatively,  the  Bauschinger 
strain)  is  larger  when  the  composite  is  first  pre-com- 
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pressed  plastically  and  subsequently  deformed  in  ten¬ 
sion  than  when  the  composite  is  first  deformed  in 
tension  and  then  compressed'". 

It  is  the  opinion  of  the  authors  that  the  afore¬ 
mentioned  differences  in  the  constitutive  behavior 
under  tensile  and  compressive  loading  as  well  as  the 
asymmetric  Bauschinger  phenomenon  are  primarily 
due  to  the  influence  of  the  residual  stresses.  It  has 
been  shown  that  thermal  residual  stresses  affect  not 
only  the  elasto- plastic  deformation  behavior  of  the 
composites'",  but  also  the  fracture  behavior"'’.  There¬ 
fore.  it  becomes  increasingly  important  to  investigate 
the  deformation  process  in  terms  of  residual  stresses. 

The  material  to  be  simulated  in  this  investigation 
consists  of  an  annealed  6061  aluminum  alloy  matrix 
reinforced  by  silicon  carbide  whiskers.  The  primary 
goal  of  this  paper  is  to  analytically  investigate  the 
generation,  arrangement  and  magnitude  of  the  resid¬ 
ual  stresses  due  to  the  difference  in  ct>efficients  of 
thermal  expansion  (JCTE)  during  cooling  process ; 
then,  to  study  the  changes  in  the  residual  stress  as  a 
result  of  external  deformation  either  in  tension  or 
compression ;  and  also,  to  investigate  the  differences 
in  the  constitutive  behavior  under  uiucrent  external 
loading,  i.  e.,  deformation  conditions.  In  the  experi 
mental  procedure  described  elsewhere"’’"”,  the  mate¬ 
rials  are  annealed  at  77.3K  and  then  tested  either  in 
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uniaxial  tension  or  compression  at  room  temperature. 
Our  commitment  in  this  investigation  is  to  adhere  to 
this  physical  process  and,  thereby,  to  simulate  the 
constitutive  behaviors  of  the  composites  and  to  moni¬ 
tor  the  change  of  the  residual  stresses.  These  residual 
stresses  were  induced  by  the  thermal  cooling  process 
due  to  the  mismatch  of  the  coefficients  of  thermal 
expansion  between  the  matrix  and  the  reinforcement, 
and  subsequently,  they  were  altered  by  the  external 
load. 

While  several  investigators  have  concentrated 
their  work  on  the  thermal  residual  stresses''"'*'*’, 
others""''"*’  have  worked  on  the  elasto  plastic  behav¬ 
iors  of  the  short  whisker  reinforced  composites  with 
emphasis  on  the  stress  strain  relationships,  i.  e.,  con¬ 
stitutive  behaviors,  of  the  composites.  Their  research 
includes  numerical  treatments  which  utilize  the  finite 
element  procedures"’’'"”.  Some  of  them  paid  special 
attention  to  the  microstructural  change  in  the  ma¬ 
trix"'’’'"’’  while  others  were  based  on  various  theoreti¬ 
cal  bases  with  different  assumptions"'"'"’*'. 

Recently,  Arsenault  and  Taya"”  and  Withers  et 
al.'“’  employed  the  Eshelby  method"" to  explain 
the  yielding  behavior  of  the  short  whisker  reinforced 
metal  matrix  composites  (MMC'  under  the  influence 
of  the  thermal  residual  stresses.  They  found  that  the 
residual  stresses  are  totally  responsible  for  the 
difference  in  the  yield  strength  of  the  whisker  rein¬ 
forced  MMC  between  tension  and  compression.  How¬ 
ever,  despite  the  effort  by  the  previous  investig- 
ators'""’''"''.  the  role  of  residual  stresses  played  on  the 
constitutive  behavior  of  the  composites  or  the  effect 
of  external  loading  on  the  formation  of  residual  stres¬ 
ses  has  never  been  systematically  investigated. 
Nevertheless,  these  experimental  or  theoretical 
results"'"'""’’  have  provided  us  with  a  clear  indication 
that  the  residual  stresses  may  be  a  dominating  factor 
w'hich  accounts  for  the  differences  of  constitutive 
behavior  as  a  funtion  of  loading  direction. 

It  is,  therefore,  of  interest  to  investigate  the 
constitutive  behavior  of  the  compxjsites  under  the 
influence  of  the  residual  stresses  due  to  the  JCTE,  and 
to  evaluate  the  changes  that  occur  in  the  magnitude  of 
residual  stresses  when  the  composite  is  subjected  to 
either  a  uniaxial  tensile  or  compressive  deformation. 
The  analytical  investigation  in  this  study  is  carried 
out  by  two  distinct  methods  :  a  finite  element  analysis, 
and  a  phenomenological  modelling  scheme  involving 
two  models.  The  main  objective  of  the  finite  element 
analysis  is  to  monitor  the  stress  strain  relationship 
and  the  changes  that  occur  in  the  value  of  the  residual 
stresses  due  to  the  external  deformation.  The  purpose 
of  the  phenomenological  models,  on  the  other  hand,  is 
to  obtain  some  physical  insight  into  the  various  fac¬ 


tors  which  control  the  processes  of  residual  stress 
generation  and  alteration  and  the  role  these  play  in 
influencing  the  composites’  constitutive  behavior. 

II.  The  Finite  Element  Model 

In  order  to  determine  the  material’s  response 
under  the  influence  of  residual  stresses,  a  simplified 
two-dimensional  finite  element  analysis  was  carried 
out  using  the  ADINA  code.  Particular  interest  was 
paid  to  determining  the  macroscopic  stress-strain 
behavior  of  the  composites  and  to  perform  a  micro¬ 
scopic  analysis  of  the  subsequent  changes  that  occur 
in  the  values  of  the  internal  residual  stresses  after  the 
external  load  was  removed. 

In  a  recent  study,  Takao  and  Taya""*’  have  shown 
that  the  use  of  the  volume  average  aspect  ratio  will 
provide  an  accurate  account  for  oveiall  composite 
properties,  and  the  employment  of  the  perfect  align¬ 
ment  will  give  rise  to  satisfactory  results,  although 
not  as  accurate  as  the  ''olume  average  aspect  ratio 
representation.  Therefore,  in  the  analyses,  the  value 
of  the  volume  fraction  of  the  whiskers  as  well  as  the 
whisker  aspect  ratio  were  held  constant  (li.  -20% 
and  /  (/--I).  The  value  of  I  d  —  \  is  typical  of  the 
average  values  observed  experimentally 
Furthermore,  it  was  assumed  that  the  composites 
consists  of  a  hexagonal  arrangement  of  arrays  of 
parallelepiped  SiC  whiskers  embedded  in  a  6061  alumi¬ 
num  alloy  matrix  (  Fig.  1  ).  Previous  investiga¬ 
tions""""  and  our  current  analysis  have  indicated  that 
within  a  particular  region,  a  small  variation  in  the 
values  of  interparticle  spacing  does  not  significantly 

affect  the  material  behavior  (e.  g..  for  the  influence  on 
yield  stress,  see  Fig.  2i.  The  value  of  T  //  in  Fig.  1, 
therefore,  was  not  selected  as  a  critical  parameter  in 
the  analysis  and  was  kept  as  a  constant  (0.29) 
throughout  the  analyses. 

Fig>jre  3  indicates  the  exact  details  of  the  finite 
element  mesh  where  six  to  eight  node  isoparameter 
elements  were  employed.  Only  a  quadrant  of  the 
composite  material  which  contains  the  SiC  whisker 
was  selected  to  represent  the  w'hole  infinite  array  of 
fibers  in  the  matrix,  and  a  multi  constraint  boundary 
condition  was  employed  accordingly  in  the  analysis 
which  corresp  )nds  to  the  bounday  condition  indicated 
as  the  following : 

where  is  the  stresses  at  the  boundary,  rij  is  the  unit 
outward  normal  and  /*  represents  four  different  sides 
of  the  boundary  of  the  unit  cell. 

Since  the  constitutive  behavior  of  the  SiC  whisk- 
er/AI  matrix  composites  has  been  extensively  inves¬ 
tigated  experimentally"’'**’'"’’,  it  will  be  appropriate 
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to  compare  our  FEM  results  with  those  obtained  from 
experimentation.  Because  the  differences  in  con¬ 
stitutive  behavior  of  the  composites  are  believed  to  be 
due  mainly  to  the  presence  of  the  residual  stresses,  by 


Fig.  1  The  whiskers  are  so  arranged  in  a  periodic  hexag¬ 
onal  packed  infinite  array  that  by  considering 
symmetry,  only  a  quadrant  of  the  particle  is  need¬ 
ed.  as  indicated  in  the  figure. 


1 


Fig.  2  FEM  results  indicate  that  the  yield  stress  of  the 
composites  is  not  sensitive  to  the  ratio  of  f',  // 
defined  in  Fig.  1, 


considering  the  initial  thermal  residual  stresses  which 
are  generated  during  the  cooling  process,  the  FEM 
model  should  reproduce  such  differences.  On  the  other 
hand,  if  such  differences  are  predicted,  then  the  resid¬ 
ual  stresses  which  are  responsible  for  producing  such 
differences  should  be  reasonably  well  predicted. 

The  FEM  simulation  is  set  up  in  such  a  way  that 
the  matrix  and  reinforcement  material  is  assumed  to 
be  stress-free  at  773K.  Then  a  temperature  change 
of  -48()K  was  applied  uniformly  on  the  FEM  mesh, 
w'hich  corresponds  to  the  cooling  process  from  anneal¬ 
ing  temperature  to  the  room  temperature  at  a  very 
slow  rate.  This  is  actually  true  for  the  furnace  cooling 
condition.  During  the  entire  cooling  process,  it  is 
assumed  that  the  SiC  whisker  reinforcement  behaves 
perfectly  elastic  in  response  to  the  incompatibility 
between  the  A1  matrix  and  the  SiC  whisker  reimforce- 
ment  due  to  JCTE.  On  the  other  hand,  the  A1  matrix 
becomes  softer  and  softer  as  the  temperature  goes  up. 
The  detailed  material  constants  are  listed  in  Table  1. 

After  the  composites  were  cooled  to  room  tem¬ 
perature,  the  FEM  mesh  was  externally  loaded  in¬ 
crementally.  Figure  4  shows  the  uniaxial  stress-strain 
curves  produced  by  the  FEM  analysis  where  the 
global  stress  and  strain  are  defined  as : 

where  u,  and  <(?,>  (in  thi'  case.  ;  is  equal  to  3)  are  the 
average  displacement  and  the  normal  stress  along  the 
boundary  whose  normal  is  parallel  to  the  uniaxial 


Fig.  .'i  FEM  mesh  where  six  to  eight -n(xie  elements 
were  employed. 
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loading  direction.  <o,>  can  also  be  defined  as  the 
applied  load  intensity,  (g),  and  V  and  /  are  defined  in 
Fig.  1. 

From  these  curves,  we  can  see  that  the  FEM 
results  are  in  reasonable  agreement  with  the  experi¬ 
mental  data"'.  For  example,  the  apparent  Young’s 
modulus  for  tension  is  larger  than  that  in  compression 
and  the  compression  yield  strength  is  larger  than  that 
of  tensile  yield  strength,  which  are  both  apparent  in 
the  experimental  curves"’. 

A  Bauschinger  test  can  also  be  carried  out  by 
FEM.  There  are  many  indicators  for  the  Bauschinger 
Effect ;  one  of  them  is  the  Bauschinger  strain  (s*), 
which  is  conveniently  defined  as  the  strain  difference 
required  for  reverse  flow  stress  to  reach  the  previous 
level  of  forward  flow  stress.  An  easy  visualization  of 
the  Bauschinger  strain  can  be  made  by  replotting  the 
complete  stress  strain  cycle  so  that  the  reverse  load¬ 
ing  portion  of  the  stress-strain  curve  was  plotted 
symmetrically  against  the  unloaded  state  form  the 
forward  straining  cycle  on  the  stress  strain  curve,  so 
that  the  complete  stress  strain  curve  (forward  and 
reverse  loading)  would  appear  on  the  same  first  quad¬ 
rant  of  the  cartisian  coordinate  system,  as  can  be  seen 
in  Fig.  5.  From  the  FEM  results,  it  can  be  easily 
revealed  that  the  direction  of  initial  loading  has  an 
effect  on  the  magnitude  of  the  Bauschinger  strain  (Sh). 
If  initially  tested  in  tension  and  then  reverse  loaded 
into  compression,  the  Bauschinger  strain  is  smaller 
compared  to  the  case  where  initial  loading  is  in  com¬ 
pression  followed  by  tension  [Fig.  5].  If  we  just  con¬ 
sider  the  reverse  loading  curves,  then  the  (7  vs.  e  for 
the  reverse  loading  curve  in  the  case  of  tension  first 
has  a  higher  yield  stress  than  that  of  tension  reverse 
loading  of  compression  first  [Fig.  5). This  has  been 


Fig.  4  Stress- strain  curves  predicted  by  FEM  where  the 
apparent  Young’s  modulus  of  the  composite  is 
higher  when  it  is  in  tension  load,  whereas  the  yield 
stress  is  greater  when  the  sample  is  under  com- 
presion  load. 


observed  experimentally"  . 

In  the  study  ot  the  internal  stresses  in  composite 
materials,  it  is  a  common  practice  to  employ  the 
concept  of  volume  average  both  experimentally'®’  '"” 
and  theoretically'"”''®''^''.  Consequently,  the  authors 
would  like  to  adopt  the  same  methodology  in  dealing 
with  this  issue.  Because  the  longitudinal  thermal  resid¬ 
ual  stress  is  much  larger  than  the  transverse  residual 
stress""’,  it  is  believed  that  the  observed  differences  in 
constitutive  behavior  are  mainly  due  to  the  longitudi¬ 
nal  residual  stress.  The  longitudinal  thermal  residual 
stress  in  the  composites  was  therefore  found  by  first 
evaluating  the  residual  stress  in  a  typical  finite  ele¬ 
ment  mesh  and  then  averaging  these  values  over  the 
whole  specimen.  The  overall  average  residual  stress  in 
the  matrix  was  calculated  through  a  simple  averaging 
scheme  given  by  the  following  equation  : 

2( 

=  (?=y  =  3)  (3) 

k 

where  is  the  stress  in  element  k.  and  .1*  is  the 
area  of  that  element. 

Using  this  simplified  averaging  scheme,  it  was 
found  that  upon  application  of  1%  plastic  strain,  if  the 
applied  loading  is  tensile,  the  average  matrix  residual 
stress  changes  from  its  initial  tensile  value  of 
47.5  MPa  to  the  compressive  value  of  -1(1.8  MPa. 
while  if  the  applied  loading  is  compressive,  the  aver¬ 
age  matrix  longitudinal  (tensile)  residual  stress  drops 
from  its  as-annealed  value  of  47.5  MPa  to  20.1  MPa  at 
^  0.3%  strain  and  then  increases  back  up  to  29.3  MPa 
when  it  reaches  1%  plastic  strain.  The  overall  change 
in  resiual  stress  as  a  function  of  applied  strain  is 
shown  Fig.  6.  It  is  clear  that  the  primary  effect  of 
compressive  loading  is  merely  to  intensify  the  residual 
stress  in  the  matrix  after  an  initial  reduction  while 
tensile  deformation  prtxluces  an  overall  reduction  in 


Fig.  5  FEM  prediction  of  Bauschinger  strain  (£^)  where 
£*  is  larger  if  the  sample  is  under  the  compression- 
first  loading  scheme. 
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the  tensile  residual  stress,  and  eventually  such  a  reduc¬ 
tion  will  lead  to  changes  in  the  residual  stress  from 
tensile  to  compressive.  That  is,  the  change  in  residual 
stress  does  not  occur  monotonically  with  increasing 
deformation.  To  recap,  in  the  case  of  compression 
deformation,  there  is  an  initial  decrease  in  the  tensile 
residual  stress  due  to  deformation,  at  a  total  strain  of 
0.3%  which  is  within  close  proximity  of  0.2%  strain 
of  yielding  and  the  tensile  residual  stress  reaches  a 
minimum  at  this  point.  Upon  further  deformation, 
there  is  an  increase  in  the  tensile  residual  stress.  In  the 
case  of  tensile  deformation,  the  tensile  residual  stress 
initially  decreases  at  a  more  rapid  rate  to  ^  0.3% 
total  strain,  and  then  the  rate  of  decrease  slows  down. 
At  a  total  strain  of  0.7%,  the  residual  stress 
becomes  compressive.  Even  though  the  total  average 
longitudinal  thermal  residual  stress  is  in  tension, 
previous  FEM  results"^'  have  indicated  that  in  whis¬ 
ker-reinforced  composite  materials,  the  residual 
stress  distribution  is  such  that  near  the  tip  of  the 
whisker,  the  matrix  is  under  compression  while  the 
rest  of  the  matrix  is  in  tension,  thus  giving  rise  to  an 
overall  tensile  average  residual  stress  in  the  matrix. 
Interestingly  enough,  our  finite  element  analysis  in¬ 
dicated  that  under  applied  tensile  load,  the  compres¬ 
sion  and  tension  zones  exchanged  their  signs  within  a 
reasonable  range  of  strain  (Fig.  7',  while  compressive 
loading  led  to  a  residual  stress  increase  in  both 
regions.  Another  persuasive  fact  to  be  noted  is  that 
while  the  global  average  residual  stress  exhibited  a 
non  monotonic  behavior,  the  residual  stress  near  the 
whisker  matrix  interface  ne%'er  showed  an  initial  drop 
after  loading  as  it  did  in  the  case  of  global  average 
residual  stress.  The  response  of  the  residual  stress  to 


Fig.  6  As  a  result  of  tension  loading,  the  arrangement 
of  residual  stress  alters  so  that,  at  the  whisker 
tip,  the  sign  of  the  residual  stress  changes  from 
tension  to  compression. 


the  external  load  near  the  interface  was  such  that  the 
longitudinal  residual  stress  increased  upon  an  initial 
compressive  plastic  prestrain  and  decreased  with  an 
initial  tensile  plastic  prestrain,  as  can  be  seen  in  Fig.  8. 

Although  the  simplified  averaging  scheme  discus¬ 
sed  above  is  easy  to  use,  it  works  best  only  with  a 
uniform  field.  It  is  evident  that  this  scheme  does  not 
reflect  the  fact  that  the  deformation  upon  loading  is 
inhomogeneous  in  the  matrix,  and  because  of  this 
deformation  inhomogeneity,  the  influence  of  the 
residual  stress  on  the  yield  point  of  that  particular 
part  of  matrix  will  vary  depending  on  the  local  defor¬ 
mation  rate  or  deformation  gradient.  The  higher  the 
deformation  rate,  the  more  important  the  role  the 
residual  stress  will  play.  Therefore,  to  best  evaluate 
the  effect  of  the  thermal  residual  stress  on  the  subse¬ 
quent  loading,  the  authors  believe  that  application  of 
a  weighted-averaging  scheme  is  more  appropriate  and 
will  lead  to  more  accurate  accounting  of  the  actual 
physical  process.  To  do  so,  the  local  loading  rate  Oo,., 
dq),  as  a  measure  of  deformation  gradient,  can  be 
directly  incorporated  into  the  average  stress  formula. 

An  averaging  formula  of  the  type 


where 

(— 1(7,.,) *  ( Oij}k\ I'M  ( ^7,^)* I  f  ( 4  *  b ) 

and 

^<7  =  4'l<.jr-(7l(  (4*c) 

was  then  used,  where  q  is  the  applied  load  intensity 
and  3,j  now  denotes  the  new  average  residual  stress  in 
the  matrix. 

From  Eq  (  4  ) ,  it  can  be  easily  seen  that  as 
compared  with  Eq.  (  3  ) .  the  effect  of  residual  stresses 
at  each  point  is  weighted  literally  by  the  magnitude  of 
the  stress  rate.  Since  the  stress  rate,  which  is  defined 
as  a  stress  increase  at  a  local  point  per  unit  external 
load  increment,  is  a  measure  of  the  local  stress 
concentration.  Eq.  (  4  )  gives  a  scalar  value  which 
indicates  an  overall  impact  of  residual  stress  on  the 
subsequent  loading.  In  other  words.  Eq.  (  4  )  puts 


t  t  ♦  t 


Fig.  7  The  residual  stress  versus  total  strain  of  a  whisker 
composite. 
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more  emphasis  on  the  residual  stress  value  at  loca¬ 
tions  where  it  will  have  higher  stress  concentration 
during  subsequent  loading,  and  at  the  same  time,  plays 
down  the  importance  of  those  with  lower  stress  con¬ 
centration.  Besides  what  is  illustrated  above,  another 
significance  of  Eq.  (  4  )  is  that  by  comparing  it  with 
Eq,  (  3 ) ,  an  indicator  of  the  degree  of  deformation 
inhomogeneity  can  be  worked  out.  Such  an  indicator 
is  important  because  in  analytical  elasto- plastic  anal¬ 
yses  involving  the  Eshelby-type  mean  field  theory,  one 
of  the  key  assumptions  is  that  the  deformation  in  the 
matrix  does  not  deviate  greatly  from  the  state  of 
homogeneous  deformation.  One  easj-  way  to  obtain 
such  an  indicator  is  to  take  the  ratio  of  the  outcome  of 
the  simple  average  [Eq.  (3)1  over  the  weighted  aver¬ 
age  [Eq.  (  4 )  ]  since  in  a  homogeneous  deformation 
process.  Eq.  (  3  )  is  equivalent  to  Eq.  (  4  ) . 

To  apply  the  above  illustration  to  the  current 
analysis,  it  was  found  that  upon  an  initial  load,  the 
ratio  between  the  two  averages  is  : 

(/=y=3)  (5) 

The  factor  1.55  is  an  indicative  of  the  degree  of 
deformation  heterogeneity  of  the  composite  material, 
which  undoubtedly  contributes  to  the  change  in  resid¬ 
ual  stress  upon  unloading.  Equation  ( 5 )  indicates 
that  the  simplified  average  thermal  residual  stress  is 
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EXTERNAL  STRESS  (MPa) 

Fig.  8  Average  residual  stress  near  the  longitudinal 
matrix  whisker  interface  as  a  function  of  applied 
stress.  The  residual  stress  is  taken  from  an  area  of 
0.3/  X  (),2d  around  the  transverse  center  line  of  the 
whisker  ear  the  interface. 


greater  than  that  produced  from  the  weighted  scheme, 
which  means  that  the  weighting  factor  is  larger  in  the 
compressive  residual  stress  region  where  the  residual 
stress  has  a  negative  value.  That  simply  means  that 
the  compression  zone  at  the  tip  of  the  whiskers  is 
responding  faster  to  the  external  loading  than  the  rest 
of  the  matrix.  In  order  to  better  understand  the  heter¬ 
ogeneous  deformation  field  in  the  matrix  due  to  the 
existence  of  the  whiskers,  a  stress  rate  (do/ d?)  dia¬ 
gram  was  constructed  for  the  composite  in  different 
parts  of  the  matrix  when  the  composite  is  initially 
loaded  in  tension,  as  shown  in  Fig.  9.  In  particular, 
three  regions  within  the  matrix  were  of  interest.  The 
analysis  indicated  that  the  average  stress  rate  within 
the  regions  A,  B  and  C  were,  respectively,  6  x  10  ^ 
(MPa/MPai,  9x10  “^(MPa/MPa)  and  13x10  '^ 
(MPa/MPa).  These  results  support  the  notion  that 
the  compressive  zone  at  the  tip  of  the  whisker 
responds  faster  to  external  loading  than  any  other 
region  in  the  matrix.  On  the  other  hand,  in  the  vicinity 
of  the  longitudinal  whisker-matrix  interface,  the 
response  of  the  internal  stress  to  applied  loading  is  the 
slowest.  Finally,  in  region  B,  where  some  interaction 
with  the  whisker  tip  is  expected,  the  response  to 
applied  loading  falls  within  a  somewhat  intermediate 
state.  The  finding  that  there  is  a  stress  concentration 
at  the  whiskers  tips  is  consistent  with  w'hat  was 
observed  experimentally.  N.  Sato  et  al.‘“'  have  report¬ 
ed  extensive  damages  at  short  fiber  tips  due  to  excess 
plastic  deformation. 

In  consideration  of  the  deformation  process  in  the 


Fig.  9  FEM  results  indicated  that  the  residual  stress  field 
can  be  characterized  into  three  regions  in  terms  of 
stress  rate,  where  at  the  whisker  tip.  the  stress 
rate  is  highest. 
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matrix,  it  is  believed  that  the  differences  in  micro 
strain  rate  in  different  parts  of  the  matrix  are  solely 
due  to  the  presence  of  the  whiskers.  That  is,  near  the 
longitudinal  whisker-matrix  interface,  deformation 
will  be  restricted  due  to  the  high  stiffness  of  the  SiC 
whisker,  and  at  the  same  time,  the  matrix  material  at 
the  whisker  tip  will  have  to  deform  more  to  accommo¬ 
date  the  global  strain.  The  microstress  rate  in  the 
initial  compression  zone  at  the  end  of  the  whisker  is 
always  larger  than  the  microstress  rate  in  the  initial 
tension  zone,  irrespective  of  the  global  strain.  It  is  also 
important  to  note  that  because  of  the  differences  in 
the  micro-strain  rates  in  different  regions  of  the 
matrix  and  because  the  changes  made  by  plastic 
deformation  cannot  be  recovered,  the  compatibility 
condition  between  different  regions  will  be  permanent¬ 
ly  altered  as  a  result  of  external  loading.  Therefore, 
the  corresponding  residual  stress  will  change  after 
unloading. 

The  finite  element  analysis  presented  in  this  paper 
has  been  performed  with  the  sole  intention  of  obtain¬ 
ing  a  better  understanding  of  the  effect  of  the  residual 
stresses  on  the  constitutive  behavior  of  dis- 
continuously  reinforced  metal  matrix  composites.  It  is 
evident  that  the  most  prominent  characteristics  of  the 
composite  constitutive  behavior  have  been  adequately 
predicted  using  such  an  approach.  .Although  the  2  D 
FEM  model  provides  a  good  economical  tool  for 
understanding  how  the  system  operates,  it  was  not  the 
intention  of  the  authors  in  the  present  investigation  to 
quantitatively  compare  the  analytical  results  with 
those  obtained  experimentally.  It  is  evident  that  the 
plane  strain  representation  of  an  otherwise  much 
more  complex  three-dimensional  array  of  whiskers  in 
a  matrix  is  an  o'.'er  simplification. 

III.  Phenomenological  Modelling 

The  results  of  the  finite  element  analysis  discus¬ 
sed  in  the  previous  section  indicate  that  the  change  of 
the  magnitude  of  the  residual  stress  induced  by  JCTE 
and  its  subsequent  change  under  the  influence  of 
mechanical  deformation  are  primarily  governed  by 
the  plastic  strain  gradients  in  the  matrix.  It  is  our 
belief  that  the  existence  of  the  initial  compression 
zone  is  a  direct  consequence  of  volume  mismatch 
between  the  whisker  and  the  matrix,  which  occurs  as 
a  result  of  JCTE  in  the  composite.  On  the  other  hand, 
the  fact  that  the  average  matrix  residual  stress  is  in 
tension  can  also  be  obtained  through  interface  shear 
load  transfer.  It  is.  therefore,  clear  that  the  processes 
of  residual  stress  generation  and  residual  stress 
change  are  basically  controlled  by  two  factors ;  the 
nature  of  the  interfacial  shear  load  transfer  between 
the  SiC  whiskers  and  the  A1  matrix  and  the  volume 


mismatch  between  the  two  constituents.  The  mismat¬ 
ch  of  the  volume  is  introduced  by  either  the  CTE 
mismatch  when  under  thermal  loading  or  the 
differences  in  Young’s  modulus  when  under  external 
loading, 

In  order  to  maintain  overall  compatibility,  the 
two  factors  must  work  together.  The  volume 
difference  is  matched  through  interfacial  normal 
stress  and  the  interfacial  compatibility  is  satisfied  by 
means  of  interfacial  shear  stress.  In  order  to  compre¬ 
hend  the  separate  roles  of  "interface  effect”  and  "vol¬ 
ume  mismatch"  on  the  residual  stress  response  after 
the  mechanical  loading,  two  distinct  analytical  models 
were  developed.  The  first  model,  subsequently  refer¬ 
red  to  as  the  "interface  model”,  is  shown  in  Fig.  10, 
which  assumes  that  the  material  consists  of  a  discon¬ 
tinuous  filament  embedded  in  the  matrix.  It  is  essen¬ 
tially  the  same  as  the  shear  lag  model  with  an  addi¬ 
tional  consideration  of  thermal  stress  and  average 
plastic  deformation.  It  is  also  assumed  that  there  is 
full  elastic  continuity  across  the  whisker-matrix  inter¬ 
face.  The  whiskers,  therefore,  remain  fulh'  bonded  to 
the  matrix,  causing  complete  load  transfer  between 
the  two  constituents.  The  primary  purpose  of  this  type 
of  modelling  is  to  understand  the  effect  of  load  trans¬ 
fer  at  the  interface  on  the  residual  stress  due  to  the 
applied  loading. 

In  the  second  model,  it  is  assumed  that  a 
discontinuous  whisker  is  embedded  in  the  matrix. 
Interfacial  sliding  is  allowed  across  the  whisker 
matrix  interface  along  the  direction  of  the  .vhisker 
axis,  but  there  is  full  elastic  continuity  at  the  tip  of  the 
whisker  between  the  whisker  and  the  surrounding 
matrix  provided  by  normal  interfacial  stress.  This 
model,  which  is  referred  to  as  the  "volume  mismatch 
model",  is  shown  in  Fig.  11.  The  main  purpose  of  this 
model  is  to  explain  how  the  volume  mismatch  affects 
the  residual  stress  arrangement  as  a  result  of  tempera¬ 
ture  change  and  external  loading. 

As  in  the  case  of  FEM  analysis,  a  bilinear  stress 
strain  materials  model  is  assumed  for  the  matrix  and 
a  perfectly  elastic  model  is  assumd  for  the  whisker  in 
both  of  the  analyses.  The  two  analytical  models  are 
discussed  in  more  detail  below. 

The  Shear  Lag  .Model 

The  simplest  and  yet  most  effective  approach  to 
account  for  the  mechanism  of  shear  load  transfer  at 
the  matrix-reinforcement  interface  is  the  "shear  lag" 
theory  introduced  by  Cox'^”.  Ever  since  the  advent  of 
this  theorj-,  it  has  been  widely  used  due  to  its  simplic¬ 
ity,  Some  investigators'^''*^'  have  claimed  it  to  be 
quite  accurate  in  predicting  mechanical  properties  of 
the  composites  ;  others'**'-'*'  have  found  it  less  attrac¬ 
tive.  In  this  section,  a  modified  "shear  lag"  method  is 
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employed  to  try  to  account  for  the  generation  of  the 
residual  stresses. 

In  developing  a  modified  shear  lag  theory,  as  in 
Ref.'“',  the  change  of  the  normal  fiber  stress  can  be 
equated  to  the  interface  shear  stress  through  an  equi¬ 
librium  condition.  Furthermore,  because  the  matrix 
material  is  assumed  to  follow'  the  bilinear  stress  strain 
model,  an  assumption  was  made  that  the  amount  of 
unrecoverable  interfacial  shear  strain  (plastic)  will 
induce  a  proportional  amount  of  the  interfacial  shear 
stress  in  the  matrix  because  of  a  constant  work 
hardening  rate.  It  is  also  assumed  that  the  overall 
constitutive  response  of  the  matrix  remains  elastic 
until  the  average  matrix  stress  state  satisfies  the 
Tresca  yielding  criterion. 

V/hen  the  composite  is  cooled  down  from  the 
annealing  temperature,  longitudinal  thermal  residual 
stresses  are  developed  in  the  matrix  and  the  whisker. 
These  will  be  denoted  as  <i7m>  and  <(Ti>,  respectively. 
Using  a  modified  shear  lag  model  which  incorporates 
temperature  effect  and  plasticity,  appropriate  expres¬ 
sions  may  be  readily  derived  for  <t7m>  and  <(Ti>.  The 
governing  equation  for  the  modified  shear  lag  model  is 
(further  details  are  given  elsewhere'^”)  : 

=//(iJr/|  \JTU  +  Jh’-Jw,)-  = 

(6) 

where  Ja  is  JCTE,  JT  is  the  magnitude  of  tempera 
ture  change,  Ju  =u  -t'  where  u  and  v  are  the  dis¬ 
placements  along  whisker  axes  at  and  away  from  the 
interface  in  the  matrix,  respectively,  Jw^  represents 
Jic  when  interface  shear  stress  ro  approaches  yielding 
( Tp).  Hw  is  the  whisker  stress  and  h  is  a  constant  to  be 
determined.  As  shown  in  Fig.  10,  D  and  d  are  the 
diameters  of  the  unit  cell  and  the  whisker,  respective¬ 


where  is  the  volume  fraction  of  the  whiskers,  £„ 
and  Euj  represent  the  Young’s  moduli  of  the  matrix 
and  whisker,  respectively,  Emp  is  the  work  hardening 
rate  of  the  matrix,  (T,vm  is  the  matrix  yield  stress,  /  is 
the  total  length  of  the  whisker  and  h  and  d  can  be 
determined  as  follows ; 

where  Gmp  is  the  tangential  shear  modulus  of  the 
matrix. 

When  a  uniaxial  loading  (ffca)  corresponding  to 
the  total  strain  (ewt)  is  applied  along  the  whisker’s 
direction  in  either  tension  or  compression,  and  the 
composite  is  then  unloaded,  it  may  be  shown  that  the 
matrix  residual  stress  (<(7^'»  as  a  result  of  applied 
loading  now  becomes : 


VwEjAE  m  —  Emp) 


tanh  ^-j 


<aS‘>=-  (10-b) 

V  U- 

s— 2y2  Gmp^^w  nn.r'i 

where  s?w  represents  the  applied  plastic  strain  and  Ec 
is  the  Young’s  modulus  of  the  composites  \Ec  =  Em{\ 
-  Vu)  + EwVu.].  Since  Emp<Em.  clearly,  the  matrix 
residual  stress  decreases  upon  tensile  loading  ( efap  >0) 
and  increases  on  compressive  loading  (£cbp<  0).  The 
overall  changes  in  the  residual  stress  remain 


ly.  If  the  temperature  differential  (IJTI)  is  such  that 
the  strains  in  noth  the  matrix  and  the  whisker  remain 
elastic,  then  Jwy  and  Vp  may  be  taken  away  from  Eq. 
'  fi  i .  and  Eq,  I  6  )  will  be  reduced  to  the  form  given 
Cox  "  with  a  new  thermal  mismatch  term  (h\Ja\ 


|J7'i).  For  20%  SiC  A1  composite,  it  is  likely  that  the 
matrix  will  be  plastically  deformed  as  a  result  of  the 
cooling  process.  Following  the  procedure  summarizerl 
by  Kelly  *"',  the  general  form  of  the  solution  is 


obtained  as  follows  : 

<f7;> 

T%  eX: 

'  M  -  Vp.  E„pV 


tanh(i7/,/2) 

G3//2) 


(7*a) 


(7-b) 


(7-c) 


Fig.  10  Schematic  diagram  for  the  interface  model. 
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monotonic. 

The  set  of  equations  stated  above  applies  equally 
to  both  cases  of  compressive  and  tensile  loading 
provided  we  consider  the  transformation  : 

ecap  ~  ifirapl 

This  means  that  <ff™>  is  symmetric  about  <(7™> 
assuming  only  that  plastic  deformation  is  considered. 
It  also  can  be  seen  that  Eq.  ( 10)  is  only  a  function  of 
applied  plastic  strain  which  means  that  elastic  strain 
will  not  alter  the  value  of  the  residual  stress  which 
contracicts  the  fact  that  there  is  an  initial  residual 
stress  drop  regardless  of  the  loading  directions 
predicted  by  FEM.  It  is,  therefore,  clear  that  the 
modified  shear  lag  model  cannot  adequately  predict 
the  effect  of  the  loading  directions  on  the  asymmetric 
behavior  of  the  average  residual  stress.  In  other 
words,  the  mechanism  of  load  transfer  through  the 
matrix -whisker  interface  cannot  have  a  predominant 
effect  on  the  average  value  of  the  residual  stress  and 
the  residual  stress  change  and,  consequently,  on  the 
constitutive  behavior.  In  fact,  the  outcome  of  this 
model  matches  the  results  of  the  finite  element  analy¬ 
sis  given  earlier,  in  that,  within  the  immediate  proxim¬ 
ity  of  the  longitudinal  matrix  reinforcement  interface, 
the  change  of  the  residual  stress  is  monotonic  (Fig.  8). 
This  signifies  that  the  influence  of  the  interface  effect 
is  primarily  restricted  to  a  very  small  region  in  the 
close  vicinity  of  the  whisker-matrix  interface  and 
does  not  significantly  influence  the  bulk  of  the  compos¬ 
ite. 

The  Volume  Mismatch  Model 

In  this  model  (Fig.  11).  it  will  be  assumed  that 
when  a  whisker  is  embedded  in  the  matrix,  no  bonding 
exists  between  the  matrix  and  the  whiskers  along  the 
axis  of  the  whisker.  Therefore,  no  load  transfer  can 
occur  between  the  whisker  and  the  matrix  at  their 
longitudinal  interface,  and  hence  the  interface  effect 
can  be  separated  from  current  consideration.  On  the 
other  hand,  full  continuity  of  stress  from  the  matrix  to 
the  whisker  is  assured  at  the  interface  of  the  whisker 
and  the  matrix  at  the  tip  of  the  whisker  so  that  the 
nitrmal  stress  continuity  would  be  provided.  It  is 
further  assumed  that  a  boundary  is  made  along  the 
extension  of  the  SiC  A1  interface  along  the  whisker’s 
axis  to  separate  the  compressive  and  tensile  residually 
stressed  region.  Considering  that  the  thermal  residual 
stress  is  in  compression  at  the  tip  of  the  whisker,  the 
material,  therefore,  may  be  regarded  as  being  com¬ 
posed  of  three  distinct  zones  ;  the  region  of  the  matrix 
which  is  under  tension  (denoted  by  subscripts  "mt”), 
the  region  of  the  matrix  which  is  under  compression 
(denoted  as  "me")  and  the  whisker  zone  (denoted  by 
“w”).  The  combination  of  the  two  regions  in  the 
matrix  (mt  f  me)  makes  up  the  entire  matrix  region 


in  the  unit  cell. 

Although  the  local  field  in  each  zone  is  not  uni¬ 
form,  since  our  interest  is  the  average  residual  stress, 
it  is  often  convenient  to  work  with  volume  averages  of 
the  non-uniform  field  in  each  zone  defined  by  Eq  (11) , 
which  is  equivalent  to  Eq.  (  3  )  defined  in  a  previous 
section  : 

<  •  >A.=  ■  )dv  (11) 

where  An  is  the  area  of  each  zone.  i.  e.,  me  and 

w. 

Considering  virtual  work  and  Gauss’  theorem,  it 
can  be  shown  that  the  following  relationship  exists  : 

jj^  c!ijdA  =  j^^  a.ntik-ijdl  (12) 

i.  e., 

AnAn  —  ^Oikyik^iy{An\\An\  (  13) 

where  is  the  normal  to  the  boundary  |.T,|. 

Furthermore,  it  is  a  common  practice  to  assume 
that  the  average  strain  and  the  average  stress 

(<<t„>/i,)  are  related  in  the  same  way  as  in  a  local  field 
by  a  set  of  effective  constants'®’. 

.\n~  CukASkl^  An  I  )  .tni  —  ft.Vm 

An~  ^\^(ykty  An~  dkt\  K  ttw)  (  Ctj  | 

(14) 

where  CT,ki  and  C,%  are  the  stiffness  and  the  tan¬ 
gential  stiffness  of  the  matrix,  and  eZ  and  are  the 
average  strains  and  stresses  on  the  yield  surfaces  in 
the  strain  and  the  stress  space,  respectively. 

T.  Yano  et  al."”'  stated  that  the  lateral  stress 
generated  by  the  difference  in  Poisson’s  ratio  is  small 
with  a  A1  All  Ni  fiber  reinforced  eutectic  composite 
system  w'hich  has  a  Poisson’s  ratio  of  !/,  =0.16  and 
=  0.33.  Since  a  SiC  'A1  system  has  a  Poisson’s  ratio  of 


Fig.  11  Schematic  diagram  for  the  volume  mismatch 
model. 
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2/sic  ^0.19  vs.  i/ai^0.33.  we  shall  adapt  their  statement. 
Therefore,  if  only  the  longitudinal  normal  stress  is 
considered,  Eq.  (14)  will  reduce  to  the  following: 


5^33^ /In 

Em 

Emp 


(Ta'-l  X,.  —  fJ.vm 
^  fTvm 


(15) 


An  =  mt.  me 

From  the  above  simplifications,  it  is  easily  revealed 
that  the  concept  of  volume  mismatch  has  been 
reduced  to  the  mismatch  between  the  whisker  and  the 
matrix  at  the  whisker  end.  Because  all  the  quantities 
we  are  interested  in  are  along  the  whisker  axis,  w’e 
may  substitute  a  simpler  notation  for  the  element  to 
be  investigated,  i.  e., 

(16) 

It  is  clear  that  Eqs.  12  and  13  define  the  relation¬ 
ship  between  average  stress  over  the  area  of  one 
particular  zone  and  the  line  average  of  the  stress  on 
the  boundary  |.4,1,  Therefore,  the  global  equilibrium 
condition  (equilibrium  between  zonesl  can  be  expres¬ 
sed  in  terms  of  volume  average  stress,  and  the  equilib¬ 
rium  is  automatically  satisfied  when  adapting  Eq.  ( 13) 
Furthermore,  during  the  entire  analysis,  the  continuity- 
condition  at  the  interfaces  between  zones  is  satisfied 
on  an  average  basis. 

From  the  requirement  of  the  condition  of  equil- 
brium  of  stresses,  it  may  be  shown  that  the  average 
residual  stre.ss  in  the  whisker  is  related  to  the  average 
residual  stre.ss  in  the  matrix  by  the  following  equa¬ 
tion  : 


^  “  yVtm)  (17) 

where  again,  I  „  refers  to  the  volume  fraction  of  the 
whiskers  and  \  I  'u  1  m  ~  I  --  I 'm  denotes  the  total 
volume  fraction  of  the  matrix  portion.  The  overall 
average  residual  stress  in  the  matrix  is.  in  turn, 
related  to  the  residual  stress  in  the  tensile  part  of  the 
matrix.  <iin:>  : 

where  Fu  is  defined  by  : 

“  IV  1  IV  ■ 

It  may  also  be  shown  that  the  residual  stress  in 
the  compressive  part  of  the  matrix  is  given  by  ; 

<n;,>-=  ,  v;,>-  (20) 

1  V  mt 

Considering  the  cooling  process,  it  is  known  that 
the  mi,smatch  of  thermal  expansion  coefficients  is 
large  enough  to  induce  plastic  deformation  around  the 
reinforcement'"”.  By  considering  equilibrium  and  com¬ 
patibility.  in  an  average  sense,  it  can  be  seen  that  the 
following  equations  can  be  obtained  when 
<  |z)7’,|  (The  details  of  the  derivation  is  shown  else¬ 


(18) 

(19) 


where”'^’)  where : 
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(21) 


Therefore,  the  thermal  residual  stress  along  the  longi¬ 
tudinal  direction  is  ; 


/  ._r\ 14i,F u- 

\\~-F.)  +  Eu\Aa\\AT\ 

L  t^mp _ J  _ _ _ _ 

A.  .  Tiki 

Em  1  1  mr  L  Emp  Eu  J 

This  means  that  within  a  certain  range  of  temperature 
variation  (|z]T-l£(.|/IT|<  |zJT,|).  plastic  deformation 
will  be  introduced  in  the  matrix.  However,  until  \AT\ 
>\ATt\.  the  matrix  will  remain  partially  plastically- 
deformed.  For  [AT[  =  772K .  it  is  both  predicted  by  Eq. 

(21)  and  experimentally  obseved  by  Vogelsang  et 
al.”"  in  terms  of  generation  of  dislocations  that  the 
matrix  is  partially  plastically  deformed.  The  predict¬ 
ed  average  residual  stress  in  the  matrix  is  given  by  Eq. 

(22) . 


It  is  of  particular  interest  to  evaluate  the  new 
residual  stresses  in  the  composite  when  a  uniaxial 
loading  (  ffca)  corresponding  to  the  total  applied  strain 
(£„>,)  is  applied  along  the  axis  of  the  whisker  (.n 
direction)  and  is  subsequently  removed.  Clearly,  the 
applied  stress  changes  the  residual  stress  pattern 
when  plastic  deformation  is  involved.  If  the  condi¬ 
tions.  as  stated  before,  are  such  that  only  portions  of 
the  matrix  have  initially  yielded  due  to  the  thermal 
residual  stresses,  it  may  be  shown  that  when  the 
applied  load  is  tensile,  the  overall  average  residual 
stress  in  the  matrix  after  unloading  is  given  by  (for 
derivation  details,  refer  to'  ’’')  ; 
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and 
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and  Ece  represents  the  all-elasic  modulus  of  the  com¬ 
posite  in  the  absence  of  any  plastic  thermal  residual 
stresses. 

When  the  applied  load  is  compressive : 


<aS‘y  =  <a'my- 


Vw  r  Ew  Emp^  tot _ 

1-  Vu,  L  (1  -  Fw)Ew+  FuEmp 


_ E  wEm _  Gca 

{\  —  Fw)Eui  +  FuiEm 


(23-d) 


Equation  (23)  clearly  indicates  that  <<j™>  is  no 
longer  a  monotonic  function.  These  equations,  there¬ 
fore.  predict  a  non-monotonic  change  in  the  residual 
stresses  in  the  composite  depending  on  whether  the 
system  is  subjected  to  a  compressive  or  tensile  load¬ 
ing. 

By  considering  applied  composite  stress  (<Tm)  and 
the  corresponding  composite  strain  (fw,  =2  £<■«,) 


Gca  2  Eci 


(24) 


where  Ec,  is  the  composite  Young’s  modulus  or  com¬ 
posite  tangential  Young's  modulus  depending  on  the 
nature  of  the  applied  strain  (sca,).  i.  e.,  elastic  or 
plastic,  as  listed  in  Table  1.  The  result  is  shown  in  Fig. 
12,  It  is  clear,  by  comparing  with  Fig.  5.  that  the  most 
prominent  features  predicted  by  the  FEM  are  re¬ 
produced  by  the  volume  mismatch  model. 

Using  Eq  (23),  the  average  values  of  the  residual 
stress  in  the  matrix  were  plotted  in  Fig.  13  after  the 
application  of  tensile  and  compressive  loadings.  For 
the  purpose  of  comparison,  the  results  of  the  finite 
element  analysis  described  before  have  also  been 
superimposed  on  the  same  diagram. 


Discussion 


The  results  of  the  FEM  modelling,  which  have 
been  duplicated  by  the  phenomenological  modelling, 
indicate  that  the  residual  stress  produced  by  the 


Table  1  Tensile  material  properties  for  different  con¬ 
stituents  of  the  composites 
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JCTE  and  cooling  is  affected  even  at  initial  levels  of 
external  deformation.  If  the  deformation  is  the  result 
of  compressive  loading,  the  residual  stress  at  first 
decreases  and  then  increases  with  increasing  deforma¬ 
tion,  whereas,  if  the  deformation  is  the  result  of  ten¬ 
sile  loading,  the  residual  stress  decreases  initially  at  a 
rapid  rate  and  the  rate  slows  down  upon  further 


(‘10’’) 


Fig.  12  Stress-strain  curves  predicted  by  volume  mis¬ 
match  model  which  are  in  general  agreement 
with  the  FEM  prediction  shown  in  Fig.  5. 


Fig.  13  Comparison  of  average  residual  stress  predicted 
by  FEM  and  ‘‘volume  mismatch"  model.  Note  the 
similarities  of  slope  change  between  the  two 
models. 
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deformation.  If  the  tensile  deformation  is  continued, 
the  residual  stress  changes  sign  from  tensile  to  com¬ 
pressive. 

From  the  phenomenological  analyses,  it  is  clear 
that  the  volume  mismatch  model  gives  qualitative 
agreement  with  the  finite  element  analysis  on  the 
general  trend  of  the  residual  stress  change.  Not  only 
the  immediate  average  residual  stress  drop  upon  load¬ 
ing  is  apparent,  but  the  slope  change  can  also  be 
predicted.  Note  that  in  the  case  of  compressive  load¬ 
ing,  there  is  a  sign  change  in  the  slope  after  the 
composite  has  started  yielding,  which  could  not  be 
predicted  by  the  interface  model.  Furthermore,  by 
comparing  Figs.  4  and  12,  the  constitutive  behaviors  of 
the  composite  predicted  by  FEM  and  the  volume 
mismatch  model  are  generally  in  agreement. 

To  understand  the  underlying  mechanism  for  the 
change  in  residual  stress,  it  is  helpful  to  retrieve  the 
argument  of  weighted  residual  stress  discussed  in 
section  II .  Since  the  ratio  of  the  simple  average  stress 
over  the  weighted  average  stress  is  greater  than  1  (1. 
55),  this  means  that  the  magnitude  of  the  weighting 
factor  (stress  rate  da3-\  'aq)  is  greater  in  the  coi  'res- 
sive  region.  This  is  equivalent  to  stating  that  in  order 
to  maintain  compatibility  between  the  whisker  tip 
region  of  the  matix  and  the  rest  of  the  matrix,  the 
compression  zone  at  the  whisker  tip  will  respond 
faster  than  the  rest  of  the  matrix,  i.  e.,  higher  stress 
(strain)  rate,  and  result  in  localized  plastic  deforma¬ 
tion  before  global  yield.  This  local  stress  (strain)  rate 
gradient  will,  therefore,  explain  the  pre-global  yield 
residual  stress  rearrangement  which  is  represented  by 
the  initial  residual  stress  drop  (Fig.  1.3). 

It  is  also  interesting  to  note  that  the  current 
Eshelby  method  does  not  predict  such  a  residual  stress 
reduction  after  an  initial  prestrain.  Withers  et  al.'“’  in 
their  recent  mean  field  investigation  implied  that  the 
residual  stress  stays  as  a  constant  within  a  pre-global 
yield  region  and  that  there  is  an  overall  monotonic 
change  in  residual  stress  as  a  result  of  prestraining. 
The  discrepancy  lies  in  that,  first,  in  the  mean  field 
theory,  the  yielding  criteria  is  satisfied  in  an  average 
sense.  However,  well  before  the  average  global  yield¬ 
ing  criteria  is  satisfied,  localized  micro -yielding  has 
already  occurred,  e.  g.,  at  the  whisker  end,  as  pointed 
out  previously,  and  incompatibilities  between  different 
parts  of  the  matrix  which  are  responsible  for  the 
generation  of  the  matrix  residual  stress  have  already 
developed.  Therefore,  residual  stress  alteration  is 
induced  prior  to  the  global  yielding.  However,  by 
considering  only  a  global  average  field,  such  localized 
microplastic  deformation  can  not  be  taken  into 
account.  Secondly,  the  ratio  of  simple  average  stress 
over  average  stress  weighted  by  the  initial  stress  rate 
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gives  a  value  of  1.55,  a  55%  deviation  from  unity.  This 
indicates  that  there  is  a  high  gradient  in  the  initial 
deformation  field  which  means  that,  in  this  particular 
case  (with  a  large  deformation  gradient  within  the 
matrix  phase)  the  assumption  of  uniform  plastic 
strain  before  and  at  the  onset  of  the  global  plastic 
deformation  is  probably  unsubstantiated.  Such  an 
assumption  may  be  employed  at  the  post -yield  region. 
In  fact,  if  the  initial  non  monotonic  behavior  of  the 
change  of  the  residual  stress  is  set  aside,  the  latter 
portion  of  the  residual  stress  as  a  function  of  applied 
strain  at  higher  strain  value  is  largely  monotonic. 

In  the  post  yield  region,  it  is  worthwhile  to  point 
out  that  in  the  case  of  compressive  loading,  there  is  a 
sign  change  in  the  slope  d<.am'>  / after  the  yielding 
of  the  composite  yielding,  which  is  predicted  by  both 
FEM  and  the  volume  mismatch  composite  model. 
This  is  an  indication  that  due  to  plastic  deformation, 
the  deformation  gradient  across  the  tension  and  com¬ 
pression  zones  in  the  matrix  becomes  totally  different 
1  order  to  maintain  the  deformation  compatibility. 
The  initial  compressive  load  will  produce  a  geometric 
compromise  between  tension  and  the  compression 
zone,  and  therefore,  there  will  be  a  reduction  in  resid¬ 
ual  stress.  After  the  matrix  material  has  begun  defor¬ 
ming  plasiicrdly  throughout,  the  two  zones  become 
more  and  more  imcompatible  as  the  compressive  load 
increases  followd  by  unloading  and.  as  a  result,  the 
residual  stress  will  increase  again.  It  is  also  important 
to  point  out  that  while  the  tendency  of  residual  stress 
change  is  in  agreement  with  both  FEM  and  volume 
mismatch  models,  the  actual  magnitude  of  the  resid¬ 
ual  stresses  predicted  by  the  analytical  model  is  much 
smaller  than  that  predicted  by  the  FEM  and  the 
relative  positions  between  tension  and  compression 
curves  predicted  by  FEM.  are  not  in  agreement  with 
the  volume  mismatch  model.  It  is  speculated  that  the 
absence  of  interfacial  bonding  in  the  volume  mis¬ 
match  model  is  the  key  reason  for  the  discrepancy. 
For  example,  the  discrepancies  of  relative  positions  of 
the  curves  of  residual  stress  after  tension  or  compres¬ 
sion  between  the  FEM  and  volume  mismatch  can  be 
explained  by  the  fact  that  the  effect  of  the  interface 
will  enhance  the  initial  residual  stress  drop  when 
under  an  external  tension  load,  and  oppose  such  a 
drop  when  under  compression.  Such  an  effect  can  be 
verified  easily  through  Eq  (10),  which  is  directly  der¬ 
ived  from  shear  load  transfer.  However,  the  pattern  of 
the  residual  stress  changes  is  solely  determined  by  the 
volume  mismatch  mechanism. 

Conclusion 

In  summary : 

( 1 )  The  predicted  change  of  the  average  residual 
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stress  is  not  monotonic  with  respect  to  the  applied  pre 
-stress  (strain),  i.  e.,  regardless  of  the  sign  of  the 
loading.  The  average  residual  stress  drops  upon  the 
initial  elastic  deformation,  and  it  will  increase  or 
decrease  depending  on  whether  the  loading  history  is 
compressive  or  tensile. 

(  2 )  The  results  predicted  by  FEM  are  generally 
comparable  to  that  predicted  by  the  theoretical  model 
which  is  based  on  the  volume  mismatch  between  the 
reinforcement  and  the  matrix.  This  indicates  that  the 
volume  incompatibility  is  the  controlling  factor 
regarding  the  change  in  residual  stress. 
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Transient  and  Residual  Stress  in  Porcelain/ 
Alloy  Strip  for  Dental  Use  as 
Affected  by  Tempering* 

Kenzo  ASAOKA‘* 

Porcelain-fused-to-metal  (PFM)  restorations  may  develop  cracks  during  process¬ 
ing  or  in-mouth  service  if  the  relative  physico-mechanical  properties  of  the  porcelain 
and  alloy  are  highly  mismatched.  The  precise  conditions  when  it  might  occur  are  not 
known.  Many  processing  and  property  variations  can  affect  the  stresses  developed 
throughout  a  porcelain-alloy  system.  To  understand  this,  a  computer  simulation  of 
stress  in  a  PFM  strip  was  conducted.  The  simulation  considers  cooling  from  tempera¬ 
tures  higher  than  the  porcelain  sag  point.  The  following  temperature-dependent 
factors  were  incorporated  :  the  elastic  modulus,  shear  viscosity  (porcelain),  and 
coefficient  of  thermal  expansion.  Cooling  rates  during  tempering  of  porcelain  were 
measured.  The  cooling  rate  dependencies  of  the  glass  transition  temperature  and 
temperature  distribution  during  cooling  w'ere  also  included  in  the  simulation.  The 
results  suggest  that  tempering  with  a  high  cooling  rate  make  up  for  compressive 
residual  stress  in  body  porcelain  and  the  high  curvature  of  the  composite  beam. 

Key  Wort/s ;  Biomatcrials.  Laminated  Construction.  Residual  Stress.  Thermal 
Stress.  Viscoelasticity.  Transient  Stress 


1.  '  atroduction 

Porcelain  fu.sed  t  metal  (PFM)  restorations 
w'ere  presented  to  the  .'ental  profession  about  three 
decades  ago.  This  tech  lologica!  advance  precipitated 
a  change  to  the  era  ot  ull  coverage  restorations  that 
serve  both  the  functio’  al  and  cosmetic  requirements 
of  patients  with  carious,  dishgured.  or  missing  teeth. 
PFM  restorations  con^  st  of  a  cast  metal  coping  onto 
which  a  porcelain  veneer  is  fired.  With  this  technique, 
opaque  porcelain  is  fu>  d  against  the  casting  to  mask 
the  metal.  The  opaque  s  heavily  loaded  with  metallic 
oxide  opacifiers  and  n  ly  be  applied  with  a  thickness 
of  0.1  to  0.2  mm.  Sirce  the  opaque  porcelain  is 
reflective,  it  must  be  t  wered  with  at  least  1  mm  of 
body  porcelain  if  rea  onable  esthetics  are  to  be 
achieved. 

Currently,  the  use  of  PFM  re.st orations  represents 
a  significant  component  of  dental  restorative  practice. 
However,  the  dental  practitioner  has  been  faced  with 
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the  clinical  problem  of  occasional  failure  of  the  porce- 
lain  alloy  system  during  function  ;  a  problem  of  small 
magnitude  in  terms  of  the  frequency  of  occurrence, 
but  of  serious  consequence  when  it  does  occur.  Fail¬ 
ures  can  be  attributed  to  a  variety  of  clinical,  process¬ 
ing  and  material  system  deficiencies.  With  the  continu¬ 
ous  presentation  of  more  and  different  types  of  alloys 
and  porcelains  for  use  by  the  dental  profession,  it 
becomes  of  increasing  clinical  importance  to  assess 
their  compatibility  with  existing  problems. 

For  this  study,  a  computer  simulation  using  a 
viscoelastic  stress  analysis  was  develop)ed  to  clarify 
the  effect  of  tempering  and  thermal  expansion  of  alloy 
on  transient  and  residual  stresses  in  the  PFM  strip. 

2.  Materials  and  Methods 
2. 1  Cooling  rate 

The  powder  of  the  porcelain  was  formed  into  a 
cylindrical  shape  (5.0  mm  in  diameter  and  10.0  mm  in 
length)  in  a  steel  die  by  dry-pressing.  Thermocouples 
of  0.13  mm  in  diameter  were  embedded  in  the  porce¬ 
lain  compacts.  The  porcelain  was  fired  according  to 
the  manufacturer’s  recommended  procedure.  After 
firing  once  in  a  furnace,  the  porcelain  was  reheated  to 
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